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A vacuum ionization spectrometer has been constructed 
in which an externally controlled, motor-driven slide brings 
either of two powder specimens into the path of the 
primary x-ray beam. Comparison of the diffracted power in 
the Debye-Scherrer rings from NaCl, NaF, LiF, and KF 
powders have been made, with a F P-54 Pliotron circuit to 
detect the ionization currents. In the studies on KF a range 
of wave-lengths from 1.538A to 4.146A has been used. On 
the short wave-length side, the expected decrease of the 
crystal structure factor was observed in KF. Beyond the K 


INTRODUCTION 
HE crystal structure factor F is defined by 


the relation 


F=f; exp [27 (1) 


i 


Ni(hx;+ky;+lz;) ], 


where the summation is to be carried out over 
all the atoms in the unit cell. f; is the atomic 
structure factor of the atom j, whose coordinates 
in the unit cell are x;y,z;. N is the order computed 
from a simple translatory lattice corresponding 
to the classification of the crystal under investi- 
gation and having one atom per unit cell. #&/ are 
the Miller indices of the reflecting plane. Relative 
values of F may be obtained by comparing the 
powers of the diffracted Debye-Scherrer rings 
from powdered samples, provided suitable pre- 
cautions are taken. If F is a known function of 
the f;’s, i.e., if the crystal structure is completely 
known, it is possible, for simple crystals, to 
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edge of potassium, the crystal structure factor rose again, 
but considerably more than predicted. Conclusions: Hénl's 
computations based on the wave-mechanical theory of 
x-ray dispersion do not correctly give the observed struc- 
ture factor differences for the (200) diffraction halo from 
powdered potassium fluoride on the long and short wave- 
length sides of the potassium K edge. The predicted 
due to K 


elements 


the structure factor decrement 


for light 


invariance of 


electrons with scattering angle is 


confirmed in our experiments. 


compute f values of the constituents of the 
crystal from a sufficiently complete set of F 
values. 

The object of the experiments reported here 
was to investigate the variation of f with A in 
the region of the K critical absorption frequency. 
We chose to study a light atom, because if the 
total number of extranuclear electrons is small, 
the structure factor due to the 
anomalous K electron behavior are a relatively 
larger fraction of the structure factor. Other 
considerations, such as the necessity of finding an 
atom having its K absorption edge very near the 
wave-length of a strong x-ray line led to the 
decision to study the behavior of the potassium 
atom in KF. The power diffracted from powdered 
KF was compared with that diffracted by NaF 
and LiF, which should show only very minute 
anomalous behavior in the region around the 
potassium K_ limit, in turn were 
calibrated against NaCl at wave-lengths far from 


decrements 


and these 
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Fic. 1. Photograph of vacuum ionization spectrometer 
with the cover lifted. 


the chlorine critical absorption edge. The potas- 
sium K absorption limit in KF has been measured 
at 3.4277A by Stelling,' and this is a different 
value than that observed in metallic potassium 
by Lindh? (3.4310). Wave-lengths in this region 
are so readily absorbed in air that a vacuum 


spectrometer is needed. 


A. EXPERIMENTAL PART 


The vacuum ionization spectrometer 

The vacuum ionization spectrometer was 
similar to the vacuum double spectrometer 
designed and constructed in this laboratory, and 
used by Parratt.* Fig. 1 shows a photograph of 
the spectrometer with the lid raised. The metal 
x-ray tube, of the type designed in this laboratory 
by Dershem,‘ is attached to the raised cover, 
upon which are mounted the diffusion pumps 
used in its evacuation. Thus the vacuum in the 
instrument itself, which is merely the lowest 
pressure attainable from a Megavac pump, acts 
as a fore-vacuum for the exhaustion of the x-ray 
tube. At the right of Fig. 1 are seen the slits, 
which were carefully machined and ground steel 
bars, placed side by side, and actuated by a 
mechanism such that the distance between them 
may be changed, but they are constrained to 


1 QO. Stelling, Siegbahn, Spektroskopie der Rintgenstrahlen 
(J. Springer, Berlin, 1931), p. 290. 
>A. E. Lindh, Siegbahn, Spektroskopie der Réntgen- 
strahlen (J. Springer, Berlin, 1931), p. 289. 
L. G. Parratt, Phys. Rev. 41, 553 (1932). 
4E. Dershem, Rev. Sci. Inst. 7, 86 (1936). 
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move parallel to themselves. Another slow 
motion device turns the entire slit assembly 
about a vertical axis, so that the central ray may 
be made to pass through the axis of rotation of 
the sleeves carrying the powder frames and the 
ionization chamber. 

At the left of Fig. 1 is seen an electric motor, 
with speed reduction gears built in, which moves 
the frame carrying the powders forward and 
backward, so that first one sample and then the 
other can{be brought into the path of the x-ray 
beam. The position of the framework carrying 
the powders can be read through a glass window 
in the top of the lid, and settings can be repeated 
with ease to 0.0025 cm. 

The front slit of the cylindrical ionization 
chamber can be seen in the central background. 
The chamber was filled with argon gas at one 
atmosphere. The outer wall of the chamber is 
charged to 90 volts; the insulated collecting 
electrode leaves the chamber at the rear, and the 
ionization current is conducted along insulated 
rods which enter the central shaft of the instru- 
ment and pass down the axis of rotation. Some 
of the ionization current readings were taken ona 
Cambridge Compton electrometer, later a F P-54 
Pliotron amplifying circuit was installed. 

The angular positions of the ionization chamber 
and powder frames can be changed by rotation 
of sleeves which extend down the central axis, 
and the rotations can be read to about 2 minutes 
of arc on a divided circular scale. A sine screw 
operating on the sleeve carrying the powder 
frames allows angular rotations of 30 seconds of 
arc to be easily made. Numerical data on the 


spectrometer are given in Table I. 


TABLE I. Dimensions of the vacuum ionization spectrometer 
From x-ray tube target to central axis 27cm 
Length, height, width, of slits 15, 1.55, 0.041 
From central axis to ionization chamber 


window 9.0 
Diameter, length, of cylindrical ionization 

chamber 3.5, 6.5 
Width, height, of ionization chamber slit 0.148, 1.0 


Experimental technique 

The substances which were used, namely, 
NaCl, NaF, LiF, and KF were obtained from 
various chemical supply houses which sell an- 
alyzed chemicals. All were stated to contain only 
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traces of foreign substances. Because of the fact 
that the presence of chlorine as an impurity in 
the fluorides would have introduced serious error, 
tests for Cl were made on the actual samples 
used, insuring its absence in appreciable amounts. 
The powders were ground for several hours in a 
porcelain mortar, and in some cases finally in an 
agate mortar; then compressed into briquets in 
steel frames. These frames were rectangular 
pieces of steel 3 mm thick and 3.2 X4.0 cm, with 
a rectangular aperture 1.62.3 cm to hold the 
powder. Enough powder was used so that after 
compression the aperture was over full. The face 
of the briquet from which diffraction was to 
occur was then eroded down to the level of the 
frame by rubbing on sandpaper. In this way we 
avoided using layers which had been immediately 
adjacent to a metal surface during compression 
and which might thus have been oriented. The 
normal way in which the powers diffracted into 
the (200), (220), and (222) rings from KF 
decreased, indicated lack of orientation, at least 
at 1.538A wave-length. 

Potassium fluoride is hygroscopic, forming a 
hydrate, KF 2H,O, which melts at 41°C. Since 
the hydrate was supplied, dehydration was 
necessary. The hydrate was heated overnight on 
asteam-bath, thus removing all but traces of the 
water without danger of hydrolysis. The dried 
material was then heated to 600°C in a platinum 
crucible. The grinding operations were carried 
out in a mortar placed on an electric heater 
which maintained a temperature of over 100°C 
throughout the operation, and the compressing 
was done in a heated press. The KF was then 
stored in a desiccator over P2O; pending use, and 
a large pan of P,O; was always kept in the 
spectrometer. Air was admitted to the instrument 
through a drying train. 

The Pliotron circuit® used in the latter stages 
of the experiment had a sensitivity of 4X10™ 
amp. per mm scale deflection. The weakest 
powder lines, such as the diffraction of Ti Ka in 
the order (222) from KF gave deflections of 1.0 
cm above the background, with the x-ray tube 
running at 30 kv and 3 ma. Strong lines, like 
Cu Ka, from (200) NaF would have given 60 cm 


ee 


*L. A. DuBridge and H. Brown Rey. Sci. Inst. 4, 532 
1933), 
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deflection at 30 kv and 10 ma if the sensitivity of 
the circuit had not been reduced. Thus the 
measured ionization currents varied from 4 
X10- to 2.4107 amp. 

In the theory of diffraction by a powder 
consisting of randomly oriented particles, one 
calculates the power diffracted into an arc of a 
Debye-Scherrer ring. Two experimental methods 
were used to measure relative values of this 
power. 

(1) The ionization chamber was moved in 
angular steps equal to the angular width of its 
window through the region in which a diffracted 
maximum appeared. The background readings 
were determined at positions adjacent to the 
diffracted beam. The sum of the differences 
between the readings in the beam and the 
background readings is proportional to the power 
which it is desired to measure. The method is 
easily seen to be equivalent to using a single 
position of an ionization chamber with a slit 
wide enough to take in the entire beam of 
diffracted x-rays. The angular width of the 
ionization chamber window (40 min.) was 
measured by moving the chamber past the 
diffracted beam from a single calcite crystal, 
which diffracts a beam of the same angular 
divergence as that leaving the slits. For the 
powder specimen three or four adjacent positions 
were always sufficient to cover the diffracted 
beam. 

(2) The ionization chamber was moved in 
small steps, 5 or 10 minutes of arc, through the 
region of the diffracted beam, and adjacent to it 
until a level background was obtained. With 
Pliotron readings as ordinates and angles as 
abscissae, a smooth curve was drawn through the 
points obtained. The area under the peak and 
above the base-line is proportional to the power 
whose measurement is desired. 

These two methods were compared many 
times and were shown to give the same power 
ratio for two diffraction rings. The technique of 
making a measurement was as follows: Let us 
suppose the relative powers of LiF (111) and KF 
(200) for a certain wave-length are required. 
The powder frame support is moved by means of 
the motor until LiF is in the x-ray beam. The 
sleeve carrying the powder frames is then rotated 
until the incident and diffracted beams enter and 
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PABLE II. Sets of readings to indicate the accuracy of repetition. 


F values based on | Fnac)(200) = 


20.80, | Fnar(200) = 14.55. 


POWDERS ACA ORDERS Dati METHOD [P. P RESULT 
LiF; NaCl 1.934 111); (200) Dec. 24, '35 2 1.66 iFycir(111)= 4.98 
53 - 5 ge 2 1.72 ” 5.08 
™ « Jan. 25, '36 rs 1.63 4.94 
Average 5.00 
KF; NaF 4.146 200); (200) Nov. 28, '35 1 0.984 1 Fx p(200) = 18.50 
7 - - Nov. 29, '35 | 0.917 a 17.86 
” ” Nov. 30, '35 1 0.939 18.08 
ss " Dec ‘. ae l 0.918 a 17.88 
3 - Dec. 3, °35 2 0.914 " 17.84 
Average 18.03 


* Fresh sample of NaCl 


leave the briquet at the same angle, the so-called 
“reflection method” being used. The ionization 
chamber is moved according to method (1) or (2) 
above, and the Pliotron circuit readings noted. 
This constitutes run A. Run B is then immedi- 
ately made on KF by bringing KF into the beam 
by means of the motor-driven slide. Going back 
to LiF, run C is made, and then a second run, D, 
on KF completes the set. The four runs could be 
made in about 90 minutes. In computing the 
power ratios KF /LiF, B/}(A+C) and }(B+D)/C 
are used. This is a rough method for correcting 
for a steady variation of power with time, such as 
might come from slow deposition of tungsten on 
the target, or on the cellophane window of the 
x-ray tube. 

As an example of the accuracy attained in 


relative power measurement, Table II shows 
some determinations made with Fe Ka and 
Ag La wave-lengths. We estimate that our 


structure factors on the short wave-length side 
of the potassium K absorption edge are accurate 
to about 2 percent; on the long wave-length side 
to about 3 percent. 


Difficulties and sources of error 

In the longer wave-length region, from 3 to 
4A, the deposition of tungsten on the target gave 
some difficulty. Two effects may result from this 
deposition, (1) a decrease in intensity of the soft 
line radiation with time, (2) an appearance of 
higher order diffractions of the tungsten L-series 
lines superposed on truly soft x-ray diffraction 
haloes. It was found desirable to use a tube with a 


shielded filament,‘ which was very satisfactory in 
removing this source of error. 

A curious behavior of NaF in the spectrometer 
was discovered in some experiments undertaken 
to detect an effect of particle size, if such existed. 
Two samples of NaF, which under a microscope 
showed different average sizes of particles, were 
placed in the spectrometer and compared. One 
of the specimens had been used for several weeks; 
the other was fresh. Contrary to expectation, the 
fresh specimen, of relatively small particle size, 
showed a smaller diffracted power in the (200 
ring for the Ag La wave-length. After several 
days standing in the spectrometer, which was 
essentially a vacuum desiccator, the two speci- 
mens showed equal power in the (200) diffraction 
haloes. This ‘‘aging effect’’ caused us to question 
all previous comparisons in which fresh specimens 
of NaF had been used. We repeated many of 
these earlier measurements with LiF the 
comparison powder. No such anomalous effects 


as 


were found for LiF. 

Orientation of the powders was tested for by 
brushing the briquet with a camel's hair brush 
between determinations. No good evidence for 
such an effect was found. 

Where the x-rays fell on the powder a dis- 
coloration was produced after some time. To test 
whether this effect was connected with a change 
in diffracted power, the samples were shifted 
slightly so that a fresh spot was selected. No 
change in diffracted power could be detected. 

Schafer® has described an effect of particle size 
which is equivalent to the statement that with 





°K. Schafer, Zeits. f. Physik 86, 738 (1933). 
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very high absorption in the particles the effective 
absorption coefficient of the material in the 
briquet is no longer directly proportional to its 
average density, which relation is tacitly assumed 
in the expression for the computation of the 
atomic structure factors, Eq. (2). Our experi- 
ments with various particle sizes in NaF at 
4.146A could not establish the presence of such 
an effect. 


Experimental results 
It may be shown that? 


Fy Me) jem. sin = (Se ~) 1 
F, M,\ juss sin A> 


1+cos? 262\ [Pa i)! 
x - — | » 4% 
sin 262 [Pale 
In this expression the subscripts refer to two 


different diffraction haloes, either from the same 
powder or from different powders. 


sin 26; 


Nm 


>] 


is the crystal structure factor, Eq. (1). 


— 
— 


is the number of unit cells per unit volume. 
is a factor showing how many orientations of the unit 


cell can contribute to the diffraction halo in question. 


= 


is the linear absorption coefficient of the particles 
composing the briquet. 


~ 
=> 


is the angle between the forward directions of the 
incident and diffracted beams. 

P, is the power in the arc of the diffraction halo entering 
the ionization chamber window. 


All of the crystals used were of the NaCl 
structure, and it can be shown® that F can only 
be 0, or 4(f.+/,), or 4(f.—f,) where f, and f, are 
the atomic structure factors of the atoms com- 
posing the crystal. In cubic crystals the inverse 
ratios of the cubes of the length of the sides of the 
unit cells (a9) can be used for M2/M;,. The 
trigonometric factors of Eq. (2) may be grouped 
together to give 


(6) = (1+ cos? 26)/sin? 6 cos 6 (3) 


for each reflection. Pertinent data on the crystal 
powders used are given in Table III. The 





"Compton and Allison, X-Rays in Theory and Experi- 
ment (D. Van Nostrand, 1935). The expression is a com- 
bination of Eqs. (6.78), p. 414 and (6.85), p. 419, re- 
membering that (effective absorption coeff.) = (coeff. for 
particles) Eidensits of powder)/(density of particles) ). 


*Compton and Allison, reference 7, pp. 354-355. 
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TABLE III. Properties of the crystals used. 


CRYSTAL lo DENSITY Vf 
NaCl 5.628A 2.164 5.610 x 10" 
KF 5.328 2.481 6.612 

NaF 4.62 2.809 10.14 

LiF 4.01 2.65 15.51 


densities are those computed from x-ray data. 

Table IV gives the results obtained for atomic 
structure factors of KF at various wave-lengths 
and various values of (sin @) \. All the structure 
factors in this table are based fundamentally on 
the value’ 


, Fy ac) (200) =frathfi |= 20.80 
(1.538A, room temperature). 


Based on this value are the two subsidiary 
standards, 


t Fwar(200) =fratfr = 14.55 
(1.538A, room temperature) 
and 


tFriv(111) =f —fij = 5.00 
(1.538A, room temperature). 


Although our experiments on the comparison of 
LiF and NaCl were actually carried out for 
Fe Ka radiation, we have assumed that the value 
found is valid for Cu Ka radiation. A correction 
for the variation of the f’s of Na and F with 
wave-length will be discussed later. 

7 
for }Fywaci(220) at room temperature, ob- 
tained }Fyaci(200)=20.80, } Fy ar(200) = 14.25, 
i Frir(111) =4.61. His LiF value does not agree 
well with that determined by us, and we therefore 
made several determinations of the LiF : NaC] 
ratio, as shown in Table II, using different NaC] 


Havighurst,'® assuming a value of 15 


samples. There is no doubt that our results are 


*R. J. Havighurst, Phys. Rev. 28, 869 (1926). 

1” R, J. Havighurst, reference 9, Table I, lists a value of 
the effective linear absorption coefficient in LiF for Mo 
Ka which is altogether too high according to recent tables 
of absorption coefficients, for instance, that due to S. J. M. 
Allen and quoted in Compton and Allison’s X-Rays in 
Theory and Experiment. These tables give uw for solid 
LiF =3.82 at 0.710A; Havighurst used 4 =4.59. The cor- 
rection to his results seems, however, to increase rather 
than to decrease the discrepancy between our structure 
factor and his. In our experiments the iron K radiation 
passed through a MnO; filter which removed almost all the 
Fe KB line. 
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] ABLI I\ = served, une rrected, crystal Structure factors. 


Based on 1 Fy 


Pow h ‘ I OF 
Nal NaCl 1.538 ( Ka 00 00 19° 25 15° 53 
Lik NaCl 1.934 Fe Ka 111 00 24° 41 0° 06 
KI NaCl 1.538 Cu Ka 200 00 16° 48 15° 53 
KI Nal 1.538 Cu Ka 00 00 6° 48 19° 25 
1.934 Fe Aa  o 24° 44 
?.499 Va Aa 7a 32° 44 
».744 Ti Ae 31° 00 36° 26 
4.353 ( Ka 39° OO 416° 32 
3.432 Sh La 10° 06 } 58 
3.593 Sn La $2 Se" 
4.146 \c Lae 51° 06° : 63° 49 
Kl Lil 3.593 Sn La 00 1 42 ; 0° 54 
3.927 Ag Lp 47 ’- 58° OF 
4.146 Av Le 31° 06’ : 63° 34’ 
KI 1.538 Cu Ka 20 200 24° 09’ : 16° 48 
Y 00 30° 00’ : 16° 48 
744 Ti Ke 20 00 46° 45 31° 00 
222 200 63° 08 31° 00 
3.353 Ca Ka 220 200 62° 49 39° 00 


TABLE V. Correction for variation of the structure factors 


LIN! ACA Mina My tf Nap(200 
Cu Ke 1.538 —(0.13 —().06 14.55 
Fe Ka 1.934 -—0.18 —().09 14.63 
Va Ka 2.499 —0.25 —0.13 14.74 
Ti Ke 2.744 —().27 —0.15 14.78 
Ca Ka 3.353 — (0.32 —0.19 14.87 
Sb La 3.432 — 0.33 —().19 14.88 
Sn Le 3.593 —0).34 —().21 14.91 
Ag LS 3.927 —().35 «23 14.94 
Ag La 4.146 —().36 —().25 14.97 

LiF use is i i 


uniformly higher than his. James and Firth" 
found a value of } Fyc\(200) = 20.65 which differs 
from 20.80 by approximately the correction for 
anomalous dispersion of the K electrons in 
chlorine at the Cu Ka wave-length. 

In spite of the fact that the K-absorption 
limits of the sodium and fluorine in our standard 
crystals are far from the wave-length region of 
the observations, the variation in their f’s with 
wave-length is appreciable and must be corrected 
for. We have used Hénl’s theory in calculating 
this correction, as given in Eq. (10). The value of 
A for sodium and fluorine has been taken as 0.30. 
The critical absorption wave-lengths of Na and F 
were taken as 11 and 18A respectively. Eq. (10) 
leads to a set of Af values for Na and F. We have 


t James and Firth, Proc. Roy. Soc. Al17, 62 (1917). 


-}(200) = 20.80. 


b(0 rm R 
15.42 : 23.92 6:6 1.0: 162.6 2.36 i} Nal x 
9.00 : 14.28 8:6 69.4: 312 1.67 tFLip( 00 
21.22 : 23.92 6:6 43: 162.6 0.688 iF 00 ” 
1 15.4 6:6 43 1.0 0.296 if’ KR (200 

4 8.95 483 39 0.299 
6 + 900 8 0.295 
5.3 RS 1 } 476 0 Ss ) 
3.39 199 65 0.18 } 

4 2.74 358 : 705 0.547 0.80 

0 5 106 5 & 

5 3.86 602 1213 0.936 8 05 
3.02 4 } 6 S 106 116 1.805 il KI oO 
74 3.13 520 : 542 1.44 0 
> 75 3.8) 602 : 638 1.25 17.46 

1 \ 

9.44 21.22 12:6 265 188 0.528 Fy 0 4 
S7 > 8-6 325 188 0.136 ) 0 
».75 5.37 12:6 65 188 0.597 0 13.86 
3.75: 5.37 8:6 325: .188 0.260 22 06 
3.70 3.39 12:6 65 188 1.25 20 


of the standard crystals NaF and LiF with wave-length 


ey ipa iF KF(200 iP KR(zZ20) LF KY 
5.00 19.14 14.75 11.70 
5.03 19.33 
5.07 18.79 
5.09 18.49 14.09 9.7 
5.13 15.18 11.52 
3.33 11.04 
5.15 16.96 iz.ee" 
5.17 17.64* 
5.19 18.57 18.11* 


computed the corrected values of 1 Fy (200) and 
| Fyir(111) as follows: 
(3 Frar(200) ], = 14.55 — [(Afwa)a— (Afsa)i. 

—[ (Afv),— (Afr) iss 
(4 Friv (111) J, =5.00—[(Afe— (Afr): 


Thus the crystal structure factors of our standards 
slowly increased with wave-length through the 
region investigated. The KF structure factors in 
the last column of Table IV were computed on 
the assumption of no variation of the NaF and 
LiF factors with wave-length; the corrected 
values are shown in the latter columns of Table V. 

The values of fx+fr thus obtained at room 
temperature were reduced to values for atoms at 


rest in the lattice by the general relation f=/ce 








00 


1.00 


11.70 


9.7 


) and 


dards 
h the 
ors in 
ed on 
F and 


‘ected 


ble V 


room 


yms at 


= foe 





v 





ATOMIC STRUCTURE 
where fy is the corrected structure factor and f 


is the factor observed at absolute temperature 7. 


Waller's 


The value of J is determined by 


6h? s(x 1\ sin? @ 
M- ( + ) | 
m ke £ 4 ” 


where n is Planck's constant, 


formula™ 


m, is the mass of the atom considered, 
k is Boltzmann's constant, 


: is the ratio O/T, where © is the characteristic 


» 


temperature of the crystal. 


The function (x) is the result of the evaluation 
of a definite integral and has been tabulated by 
Debye.'* The value of the characteristic tempera- 
ture for a crystalline substance may be deter- 
mined Lindemann 
melting point formula.'* In this manner the 


approximately from the 
characteristic temperature of KF was found to 
be 297°K. In reducing the f’s to fy’s, a tempera- 
ture correction was computed separately for each 
of the two ions in KF, by using the proper mass 
of each ion for m,. The values of e~™ for each ion 
are shown in columns 2 and 3 of Table VI. Since 
in the experiments no diffracted beams from 
planes in which } Fxr=fx —fr were measured, no 
values of either fx or fr could be obtained 
directly from experiment. A value of fy? was 
therefore taken from the Hartree'® curve for the 
corresponding value of (sin @) A. The fr was then 
computed from this Hartree value of fi? by 
multiplication with the factor e-“r. The re- 
sulting value of fy was then subtracted from the 
observed value of fx+fr, leaving a value of fx. 
This value of fx was then reduced by multi- 
plication with e”’r to fy* and the above value of 
fo® added to it, giving }FyS®=fp¥ +f,". The fifth 
column of Table VI therefore represents the final 
values of } F)"*® from our experiments, and thus 
those used to compare with theory. 

It should be noted that although a theoretically 
calculated value of fy” is introduced into the 
temperature correction, the final result, }F)*?, 
has but a negligible dependence upon the value 





* [. Waller, Zeits. f. Physik 17, 398 (1923), Diss. Uppsala 
1925). 

'’ P. Debye, see Compton and Allison, reference 7, p. 437. 
“Lindemann, Zeits. f. Electrochemie 17, 817 (1911). 


ye D. R. Hartree, see Compton and Allison, reference 7, 
p. 781. 
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raBLe VI. Temperature corrections applied to the values 
iF of Table J 


| } 
‘Ki Vv 


TableV e«° A ans Ilartree j/ _ I : ACA 
19.14 0.964 0.930 6.90 20.10 200) 1.538 
19.33 . . ' 20.30 1.934 
18.79 7 19.74 2.499 
18.49 19.43 2.744 
15.18 “s a " 15.98 3.353 
11.04 we o op 11.69 3.432 
16.96 ng " - 17.84 3.593 
17.73 “ “ ” 18.64* 3.593 
17.64 ” ' . 18.54* 3.927 
18.57 ‘ “ “ 19.51 4.146 
18.11 2s 35 19.03* 4.146 
14.75 0.933 0.867 5.50 16.21 220) 1.538 
11.70 0.900 0.807 4.43 13.45 (222) 1.538 
14.09 0.938 0.867 5.50 15.49 220) 2.744 
9.7 0.900 0.807 4.43 11.2 222) 2.744 
11.52 0.938 0.867 5.50 12.74 220) 3.353 


of fo” which is introduced. Such a value need be 
known only approximately, since e~”« and e~™ 
differ but little for any of the planes considered. 
Should another value of f)* differing by one full 
electron be chosen, the effect on the final value of 
1Fo¥® would be 0.1 
diffraction, and about 0.04 electron for the (200). 


electron for the (222) 


B. THEORETICAL PART 
Statement of the theories 


The theory of the variation of the atomic 


structure factor with wave-length has _ been 
developed by Coster, Knol and Prins,'* Glocker 
and Schafer,'? and Hénl.'* As the frequency of 
the primary radiation decreases and approaches 
the critical frequency of the K electrons, the 
amplitude of the wave scattered by them 
increases, becoming much larger than that of the 
J. J. Thomson scattering. At the same time the 
phase of the scattered wave changes with respect 
to that of the incident wave. As a result the 
scattered wave from the K electrons interferes 
destructively with the waves scattered from the 
remaining electrons in the atom, and reduces the 
resultant amplitude observed at any distant 
point. The atomic structure factor for an atom 


is the ratio of the wave scattered by the atom at 


6 Coster, Knol and Prins, Zeits. f. Physik 63, 345 (1930). 

’ Glocker and Schafer, Zeits. f. Physik 73, 289 (1931) 

's Hénl, Zeits. f. Physik 84, 1 
18, 625 (1933). 


1933); Ann. d. Physik [5} 
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O to a distant point P, to the amplitude which 
would be scattered to P by a single electron at O, 
the single electron assumed to be scattering in 
the J. J. Thomson manner. In order to account 
for the type of absorption observed in the x-ray 
region, it has been necessary to replace the K 
electrons by a distribution of virtual oscillators 
=x to 


whose natural frequencies extend from v= 
the frequency of the K limit. If we empirically 
distribute these oscillators in such a manner that 
the absorption coefficient increases as \*, and 
retain the rest of the framework of the Lorentz 


dispersion theory, we find 


2amc? 
J(O= 
e*h? 


In this equation, f(0) is the atomic structure 
factor in the forward direction (scattering angle 
zero). gx is the oscillator strength of the K 
electrons of the atom of atomic number Z. 5x is 
the unit decrement of the refractive index per K 
electron and x is the index of absorption per K 
electron. The equation assumes that only the 


e?r? 2 


KOx+(Z—Sx) + 


2am 2| 


K electrons are resonating in the atom. 
On the short wave-length side of the K limit 


e*d?2 
6x= {1-+x? In (x-?-—1)—2aryx*! (5) 
2amc* 
and on the long wave-length side 
ed? 
oxk>= f1+x? In (1—x-*)!. (6) 
2rmc? 


In these equations, x=\/Ax, and y is a damping 
constant which we have estimated from the 
width of the K series lines of potassium. By 
extrapolation of Parratt’s'® values of the Ka line 
widths of V (23), Ti (22), and Ca (20), we estimate 
that the Ka lines of K (19) have a full width at 
half-maximum of 1.67X10-" cm. The value of 
y is this width divided by the wave-length (3.73 
X 10-8 cm) or 4.51074. 

We have computed values of 8x according t 


) 


19L. G. Parratt, Phys. Rev. 44, 695 (1933). 
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the following formula: 


hA r-!l 


OK = 
4arNox r 


Here A is the atomic weight of potassium, NV jis 
the Avogadro number, (6.023 1075), r the 
absorption jump ratio of the potassium K edge, 
(10.2), and wu,» is the mass absorption coefficient 


Is 


for potassium for wave-length \<Ax. 8x =0 for 
A>Ak. 

In the earliest form of x-ray dispersion theory, 
due to Kallmann and Mark,”’ the factor zx which 
enters the preceding equations was set equal to 2, 
the electron population of the K shell, according 
to the Pauli principle. More recent development 
of the dispersion theory has indicated that for the 
K electrons at least, an oscillator strength less 
than 2 should be observed. The reason for this 
has been clearly stated by Kronig and Kramers.2 
In an atom having only one extranuclear 
electron, for instance, hydrogen, in its normal 
state, the sum of the oscillator strengths for all 
possible transitions, including those to states of 
positive energy in the continuous region, is one.” 
For a two electron atom however, the second 
electron may remove one of the possibilities of 
transition for the first one, by occupying one of 
the levels to which it could possibly jump. This 
subtracts a certain oscillator strength from the 
original electron, and also subtracts a corre- 
sponding negative oscillator strength from the 
second electron, because for it a possibility of 
emission removed. Thus the decrement in 
oscillator strength is transferred from the first to 


is 


the second electron. In the case of the K elec- 
trons, practically all the important discrete 
levels to which they might go in absorption are 
already full, and their oscillator strength corre- 
sponds to that calculated for transitions into the 
continuous region alone. 

An important quantity in these considerations 
is Af, which may be defined as the structure 
factor decrement. Since the normal value of 
f(0) is Z, the decrement Af is Z—f(0). By a 


20 Kallmann and Mark, Naturwiss. 14, 648 (1926); Ann. 
d. Physik 82, 385 (1927). 

1 R. de L. Kronig and H. A. Kramers, Zeits. f. Physik 48, 
174 (1928). 

22? W. Kuhn, Zeits. f. Physik 33, 408 
Naturwiss. 13, 627 (1925). 


1925): W. Thomas, 
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ATOMIC STRUCTURE 


combination of the preceding equations, it may 
be shown that on the short wave-length side of 
the limit 


Af=Z—| [sx(x* In (x-*—1) —2ryx*)+Z}? 


m*c!A? (r—1)? : 
+ Ibm* (8) 
$r2e4N? or? 


and on the long wave-length side 
Af=—2«x7 In (1—x-?). (9) 


These equations are often referred to as those 
coming from the Kramers-Kallmann- Mark treat- 
ment of dispersion theory. 

Hon! has made a more fundamental theoretical 
investigation of the problem than those which led 
to Eqs. (8) and (9). He computed from theoretical 
grounds the variation of the absorption coeffi- 
cient with wave-length, instead of arbitrarily 
inserting a A* variation law, and computed the 
oscillator strengths zx from hydrogen wave 
functions, making approximate corrections for 
screening. Hénl’s equation for Af is 


2’e 4x7 In | 1—x-? 
Af= —_- 
9 (1—A)? 


1 | 1-—x | 
— \x* In +2x? | (10) 
(1—A)3| 1+x 


where A=0.246 and represents a correction for 
screening in the potassium atom. The e appearing 
in Eq. (10) is the basis of natural logarithms. 
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For large values of x, Eq. (10) approaches the 
value Af=1.43, which is the oscillator strength 
of the K electrons of potassium. The values of Af 
computed from Eq. (10) are very close to those 
computed from Eqs. (8) and (9) with sz, =1.43. 

Hénl has also investigated the possibility that 
the Jf values may depend on scattering angle, or 
on (sin @) A. He finds that although such an 
effect is theoretically indicated, the magnitude 
is so small as to be undetectable experimentally 
except possibly for very heavy atoms. This 
makes it possible to use the structure factor 
decrements of Eqs. (8) and (9), computed with 
Z=19, at the actual scattering angles which can 
be experimentally realized. 


Comparison of theory and experiment 


We shall first discuss the comparison of theory 
and experiment for the variation of }F,(200) 
with wave-length in KF. In making this com- 
parison, we assume that our measured }F)(200) 
contains the A[} F)(200) ] calculated from Hénl's 
theory for potassium and fluorine at the wave- 
length of Cu Ka (1.538A). We have 


A} Fo(200) J = AfoX + Afv’. (11) 


The left-hand member of the preceding equation 
represents the decrement in the reduced crystal 
structure factor; the right-hand members repre- 
sent the decrements in the reduced atomic 
structure factors of potassium and fluorine, and 
are to be computed directly from Eqs. (8), (9), or 
(10). The computed value of Af}F,(200)) at 
1.538A is —0.38, according to Hénl's theory. Our 


TABLE VII. Comparison of observed and theoretical values of \(\ Fy) in the (200) diffraction of K F. 


Eqs. (8 

Lint ACA) v=A Ay Mm ZK =! 
Cu Ka 1.538 0.4486 137 —0.38 
Fe Ka 1.934 0.5643 272 —0.31 
Va Ka 2.499 0.7289 540 —0.03 
Ti Ke 2.744 0.8006 656 0.36 
Ca Ka 3.353 0.9785 1128 3.87 
KK abs 3.4277 1.0000 1183 11.02 
Sb La 3.432 1.0013 (125) 8.93 
Sn La 3.593 1.0482 141) 3.80 
Ag L8 3.927 1.1146 (179) 2.71 
Ag La 4.146 1.2096 (204) 2.42 
— * 20 1.43 


7 and 8 from col. 4, Table \ 


+4 Values of Afo” used in computing cols 
LiF used as intermediate standard. 


9 AfoX MAF o)# LF y)# 
43 Honl, E 10 Eqs. (8), (9 Eq. (10 EXPERIMENT 
—().32 — 0.44 —0.38 —().38) 
—),20 —0.40 —().29 —(0.58 
0.09 —0.16 — 0.04 —0.02 
0.58 0.21 0.43 0.29 
4.15 3.68 3.96 3.74 
10.83 
7.78 8.74 7.59 8.03 
3.59 3.59 3.38 1.88, 1.08* 
2.56 2.48 2.33 1.18* 
2.33 2.17 2.08 0.21, 0.69* 
1.43 
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Fic. 2. The wave-length variation of the crvstal structure 
factor for the (200) diffraction of powdered KF. Hénl’s 
theory is the solid curve. The variation of the atomic 
structure factor of fluorine in this region is so small that 
practically the entire effect may be ascribed to the potas- 
sium K electrons. 


measured value of 4)(200) here is 20.10, 


according to Table VI; hence the value after 


“correction” for resonance of the K electrons 
of potassium and fluorine would be 19.72. 
Thus having fitted theoretical and observed 


Af} F)(200)] values together at this one wave- 
length, we can compare theory and experiment 
at other wave-lengths. Thisis shown in Table VII. 

The experimental values of A(}/)) shown in 
the last column of the table were obtained in the 


following way: 
Al} Fy(200) ] = 19.72 — [1 Fy(200) }y. 


The crystal structure factors were taken from 
column 5 Table VI. The extent of the agreement 
between theory and experiment is shown in Fig. 
2. It is seen that on the short wave-length side of 
the limit the agreement is within experimental 
error, but that large discrepancies are found on 
the long wave-length side. The observed decre- 
ments for Sn La, Ag Lg, and Ag La are con- 
siderably less than half those predicted. Taking 
only the long wave-length side of the curve into 
consideration, our results agree fairly well with an 
oscillator strength of the K electrons of about 0.6. 
Inserting such a result into the equations, how- 
ever, completely disrupts the observed agreement 
between theory and experiment at Ca Ka on the 
short wave-length side, and also fails to predict 
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Fic. 3. Variation of the observed structure factor decre 


ments with scattering angle. The point for Ti Ka (222) 
contains a very large experimental error. 


the experimental value for Sb La, which is very 
close to the limit. 

It is also interesting to note the extent of the 
agreement of our results on }F)(200) with the 
computed structure factors of Pauling and 
Sherman” and of Hartree.*4 From Table VIII it 
is seen that our measured factors are somewhat 
lower than 1.538A. In 
computing the theoretical factors in Table VIII, 
0.38 has been added to the tabulated values in 
the references to allow for the resonance of the K 


the theoretical ones at 


electrons of potassium and fluorine at this wave- 
length. In computing the structure factor for 
atoms according to Hartree it has been assumed 
that the K and K* curves are indistinguishable 
at our values of (sin @) A. Our (200) factor is 
within 0.2 electron of those computed for atoms 
according to Hartree, but our values fall off more 
rapidly with (sin @)/A than do his. 

Fig. test of the 
dependence of A(}Fo) on (sin @) ». The solid 
curves represent results taken from column 5 of 


Table VI. Hartree's 


3 shows an experimental 


The dashed curve shows 


observed and theoretica 


TaBLe VIII. Comparison of 
atomic structure factor sums at Cu Ka, }=1.538A 


PAULING AND SHERMAN HARTRI EX 
INDI lons Atoms lons Atoms MEN 
200) 22.0 21.1 20.9 20.3 20.10 
220) 18.0 17.5 17.4 17.2 16.21 
222) 15.5 15.2 15.1 14.9 13.45 


Pauling and Sherman, Zeits. f. Krist. 81, | 


1932) 


4D. R. Hartree, see Compton and Allison, reterence / 


p. 781. 
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Fic. 4. Results of other observers on the wave-length 
variation of atomic structure factors in the region of the K 
absorption edge. The solid curve is Hénl’s predicted decre- 
ments for iron. The points are experiments on chromium, 
iron, nickel, and copper. 


values for ions at the Cu Ka wave-length. If the 
structure factor decrements are independent of 
(sin @) A, the solid curves of Fig. 3 should be 
parallel. This is seen to be very nearly the case 
with the exception of the curve at the wave- 
length of Ti Ka between (220) and (222). The 
(222) diffraction halo of Ti Ka was the weakest 
measured in the entire experiment. At the 
maximum sensitivity of the Pliotron circuit it 
showed a deflection of about 1 cm above the 
rather high base-line. Our experimental error 
for this point is so great that we believe the 
apparent deviation from parallelism for this 
point is of no real significance. Our conclusion 
is that the structure factor decrements are 
independent of (sin @) A’ as predicted for light 
atoms by Hénl. 


C. DiscussION OF RESULTS 


Résumé of previous work on this problem 

Since no previous determination has been 
made of the variation of the atomic structure 
factor with wave-length for any atom of as low 
an atomic number as potassium, no direct 
comparison of our results with those of other 
experiments is possible. However, some experi- 
mental data are available for atoms of somewhat 
higher atomic number and the agreement be- 
tween such data and Hdénl’s calculations is 
shown in Fig. 4. 

From Eq. (10) it may be seen that Jf is 
dependent upon the variable \/Ax and also on 
the parameter A which is a correction for 


screening. Since the latter is a very slowly 
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changing function of atomic number, one might 
consider it almost constant throughout a very 
limited range of elements of adjacent atomic 
number. Then Af would depend to an extremely 
good approximation upon the variation in \, Ax 
alone. The closeness of this approximation is 
shown in the table in Hénl’s*® paper. Making 
use of this fact, the previous experimental data 
for 24Cr, 26Fe, 28Ni, 29 Cu, have been con- 
densed into a single plot in Fig. 4. The solid 
curve, which is actually Hénl’s curve for iron, 
represents a very close approximation to the 
predictions of the theory for all these elements. 

As representative of previous work covering a 
large range of wave-lengths, we have shown the 
results of Rusterholz** on copper (+), of Schafer?’ 
on iron (X) and chromium (O), and of Bradley 
and Hope* on the intermetallic compound FeAl 
(4). All of these measurements were made by the 
photographic method. The values of Rusterholz 
are from measurements of the (220) diffraction of 
metallic copper for various wave-lengths, all, 
however, on the long wave-length side of the 
Cu K absorption edge. They show good agree- 
ment with the calculations of Hénl in the vicinity 
of the limit, if we assume a value of f,“"(220) = 16, 
as indicated by Hartree for the limiting value 
approached by the atomic structure factor of 
copper as x approaches zero. At longer wave- 
lengths than the copper K limit region, Ruster- 
holz’s experimental results fall below the Hénl 
curve. Schafer’s results for iron and chromium 
are in fairly good agreement also, but the 
experimental points seem a little higher than 
theory demands in the immediate neighborhood 
of the absorption edge on the long wave-length 
side. This discrepancy is probably within the 
experimental error. The values of Bradley and 
Hope are distinctly higher on the long wave- 
length side and rise to a value equal to the 
maximum value on the short wave-length side. 

Scattered points are also given on the plot 
showing results for a less complete range of 
wave-lengths. Both Miss Armstrong” and 


Hénl, Ann. d. Physik [5] 18, 625 (1933), p. 651 
* A. A. Rusterholz, Zeits. f. Physik 82, 538 (1933) 
K. Schafer, Zeits. f. Physik 86, 738 (1933). 
‘A. J. Bradley and R. A. H. Hope, Proc. Roy. Soc 
36, 272 (1932). 
* Alice H. Armstrong, Phys. Rev. 34, 931 (1929) 
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Wyckoff*® have measured the atomic structure 
factor for metallic copper ( a) for the wave- 
length 1.539A using the ionization spectrometer, 
and Wyckoff has also made measurements with 
the wave-length 1.655A. Since no temperature 
correction was applied in either case, proper 
temperature corrections were made using the 
values of e” given by Rusterholz. Since determi- 
nations were also made in these two papers of the 
value of the atomic structure factor of copper 
using Mo Ka (0.710A) radiation, the Af values 
plotted above are based upon the experimental 
value obtained with Mo Ka. This latter value 
was made to coincide with the Hénl Af curve 
just as in the treatment of our own experiments. 
A NaCl sample was used as a standard by Miss 
Armstrong and by Wyckoff throughout the 
range of wave-lengths which they investigated. 
No correction has been applied for the small 
change of the F value of this standard due to 
resonance of the K electrons of chlorine and 
sodium. A recent determination of the atomic 
structure factor of copper for 1.539A radiation 
has been made by Brindley and Spiers*' by a 
photographic method. Their value, which is not 
plotted, would almost coincide with that of 
Wyckoff. 

Finally there is shown on the plot a series of 
values of fy (e) taken from the dispersion curve 
obtained by Lameris and Prins® using the 
method of total reflection of the general x- 
radiation from a mirror of metallic nickel. These 
readings indicate a steeper rise of the structure 
factor on the long wave-length side of the 
absorption edge than the Hénl theory predicts. 
Direct measurements (not plotted) of the atomic 
structure factor of nickel by Wyckoff are in good 
agreement with the points of Lameris and Prins 
in this region. 

Thus it may be said that taken as a whole, the 
experiments undoubtedly show that the be- 


"R. W. G. Wyckoff, Phys. Rev. 36, 1116 (1930). 

1G. W. Brindley and F. W. Spiers, Phil. Mag. 20, 865 

(1935). 
** A. J. Lameris and J. A. Prins, Physica 1, 881.(1934). 


AND W. P. 


JESSE 


havior of an atomic structure factor at the K 
edge is qualitatively that predicted by the 
modified dispersion theory that has been applied 
to x-rays. As to quantitative agreement, a 
general statement is difficult to make. Of the 
extended surveys, Rusterholz 
Schafer seem to confirm the theory in the region 


those of and 
of atomic numbers 24 to 29, whereas those of 
Bradley and Hope and Lameris and Prins lie 
definitely above the predicted values. 


CONCLUSION 


Before making a definite statement of our 
conclusions, we wish to call attention to the 
manner in which we have treated our experi- 
mental data. We have endeavored to keep ou 
values of A(}F)) on as purely an experimental 
basis as possible. Thus we have not obtained 
them by subtracting our measured values from a 
theoretical atomic structure factor sum computed 
by the method of Hartree or of Pauling and 
Sherman. Our A(}F)) values of column 9, Table 
VIII represent differences between an experi- 
mental value of }F) at 1.538A and our experi- 
mental values at other wave-lengths. A constant 
has been added in so that H6nl’s computed 
A(iFo) values are fitted to ours at one wave- 
length (Cu Ka). 

In conclusion we may state the following: 


(1) Honl’s calculations do not correctly give the 
observed structure factor differences for the 
(200) diffraction from powdered potassium 
fluoride on the long and short wave-length sides 
of the potassium K edge. 

The predicted invariance of the structure factor 
decrement due to K electrons with scattering 
angle for light elements is confirmed in our 
experiments. 


i) 
— 


The type of discrepancy between theory and 
experiment is such as to indicate that a revision 
both of the oscillator strengths and the oscillator 
distribution functions of the K electrons of the 
potassium atom in solid KF is necessary. 








the 


ied 


the 
ind 
ion 
of 


our 
the 
eri- 
our 
ital 
ned 
ma 
ted 
and 
ble 
eri- 
eri- 
ant 
ted 


ive- 


the 
the 
um 
des 


ctor 
ing 
our 


and 
ion 
tor 


the 








APRIL 1, 1936 PHYSICAL 


REVIEW VOLUME 49 


K X-Ray Absorption Spectra of Some Compounds of Bromine and Rubidium 


S. T. STEPHENSON, State College of Washington 


(Received January 17, 1936) 


The K x-ray absorption spectrum of Br in KBr, RbBr, CdBre, and NaBrO; and of Rb in 


RbBr is obtained by using powder absorbing screens with a double crystal ionization spec- 


trometer. Methyl iodide is used in the ionization chamber to eliminate the effect of Br in methyl 


bromide. The width of the Br edge is about 9 volts and of the Rb edge about 11 volts. The 


positions of secondary fine structure discontinuities are given and discussed. No discontinuity 


due to a possible simultaneous ejection of K and L electrons is observed. 


INTRODUCTION 


+ lap fine structure and width of x-ray 
absorption edges and the secondary fine 
structure often found on the short wave-length 
side of these edges have been investigated 
recently by several observers. According to 
Kronig' the observed secondary fine structure is 
dependent upon the crystal lattice of the ab- 
sorber. It was thought that it would be of 
interest to investigate the AK x-ray absorption 
spectrum of bromine in a few bromine com- 
pounds and of both bromine and rubidium in 
RbBr for two reasons. First: The universal use 
of AgBr photographic plates when the photo- 
graphic method is employed or the almost 
universal use of methyl bromide in ionization 
chambers when the ionization method is em- 
ployed raises some question as to whether 
observed effects in the bromine absorption 
spectrum are due to bromine in the absorber or 
to bromine in the detecting plate or chamber. 
Second: In certain cases*: *»* two elements in 
the same compound do not exhibit the same 
secondary fine structure, making a study of both 
the Rb and Br absorption spectra in RbBr of 
importance since both edges are comparatively 
close together and, therefore, easy to investigate 
and compare. 


APPARATUS AND METHOD 


The double crystal ionization spectrometer is 
of the type described by Ross,* constructed from 


'R. de L. Kronig, Zeits. f. Physik 70, 317 (1931); 75, 
191 (1932) 
?V. P. Barton and G. A. Lindsay, Phys. Rev. 46, 362 
(1934). 

*D. Coster and G. H. Klamer, Physica 1, 145 (1934). 

*G. P. Brewington, Phys. Rev. 46, 861 (1934). 

*P. A. Ross, Rev. Sci. Inst. 3, 253 (1932). 


a good Societé Génevoise optical spectroscope 
Crystals A and B were split from the same piece 
of optical quality calcite and were mounted with 
faces parallel to each other by the optical 
telescope method using a Gauss eyepiece. Crystal 
B, mounted on the prism table, is moved by a 
rigid lever arm run by Brown and Sharpe 
micrometer calipers, and the angular motions 
are further checked by viewing with a good 
telescope a millimeter scale as reflected from a 
mirror waxed to the top of the crystal. The 
angular readings can be relied upon to within 3 
seconds of arc or about 0.05 X.U. in the region 
investigated. The ways supporting crystal A 
were taken from a micrometer microscope and 
fastened rigidly to the spectrometer in place of 
the collimator. A quadrant electrometer working 
at 3000 mm/volt sensitivity on a scale 1} meters 
distant measures the ionization currents in a 
small steel ionization chamber. It should be 
emphasized that here, as well as in a previous 
investigation of KBr,®: 7 methyl iodide is used 
in the ionization chamber. Fig. 2 shows the 
Mo Kajae doublet in the (1,1) position, run to 
insure that the apparatus was working properly. 

The metal x-ray tube* with Mo or W target 
was run on the pumps at about 40,000 volts 
peak and 30 to 40 milliamperes tube current. 
In view of Bearden’s® experience with oil 
diffusion pumps and of their increasing popu- 
larity, it may be of interest to state that a 
regular General Electric X-Ray Universal broad 
focus filament has been used in the tube for over 


°S. T. Stephenson, Phys. Rev. 44, 349 (1933). 

*C, D. Cooksey and S. T. Stephenson, Phys. Rev. 43, 
670 (1933). 

* Constructed, with some minor changes, following the 
design of N. S. Gingrich, Rev. Sci. Inst. 1, 691 (1930). 

* J. A. Bearden, Rev. Sci. Inst. 6, 276 (1935). 
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200 hours and shows no signs of failure. The tube 
is pumped with one single stage metal diffusion 
pump!’ backed by a Megavac with one charcoal 
trap between fore pump and diffusion pump, no 
traps of any kind being used between diffusion 
pump and tube. Apiezon oil B is used and all 
wax joints are of Apiezon wax. Doubtless a CO. 
or liquid-air trap would make the vacuum 
conditions still better but no trouble has been 
experienced and faster pumping is obtained 
without such traps. 

The absorbing screens were made by saturating 
cigarette papers with a solution of the salt and 
drying over a hot plate.‘ Seven or more papers 
were used in each screen to secure uniformity. 
Various precautions were taken to eliminate any 
apparent structure due to (1) nonuniform screens 
or (2) to crystal imperfections or (3) to variations 
in primary intensity by (1) moving the screens 
between runs, (2) taking blank runs without the 
absorber present ; also sliding the crystals so that 
the same portions were not always used, and 
(3) removing the absorber after each reading to 
check the intensity without absorber on at least 
one run for every substance except CdBre. All 
readings were made with the crystals in the 
(1,1) position, the absorbing screens being placed 
between crystals A and B. Six readings were 
usually taken to obtain each point in a given run. 


RESULTS AND DISCUSSION 


The curves presented in Figs. 1 and 2 show the 
results obtained. The ionization currents in 
mm/sec. due to energy transmitted through the 
absorber are plotted as ordinates and X.U. are 
plotted as abscissae starting from an arbitrary 
zero, no attempt having been made to measure 


TABLE |. Widths of absorption edges and positions of 
secondary structure in volts. 


EpGE COMPOUND WiIpTH \ a B 8 ( 
Br KBr 9+1.5 63 12 22 3352 
Br RbBr 94+1.5 63 10.5 17 

Br CdBry 9+1.5 63 I11 22 

Br NaBrO; 9.5+1.5 6.5 11 19 48 77 
Rb RbBr 1141.5 5.8 13(?) 27(?) 


10 The pump is of the type described in detail on page 80, 
D. H. Sloan, R. L. Thornton and F. A. Jenkins, Rev. Sci. 
Inst. 6, 75 (1935). 
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Fic. 1. The K absorption spectrum of Br in KBr, RbBr, 
CdBre and NaBrQs. 
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Fic. 2. The K absorption spectrum of Rb in RbBr. Mo 
Kaa, doublet in inset. 


the actual wave-length position of the edge itself 
or any shift of the edge due to chemical combi- 
nation, though this latter could not have ex- 
ceeded 3 volts between the compounds studied. 
In Fig. 1 each of the curves following KBr is 
displaced an arbitrary amount downward and 
approximately 1 X.U. to the right of the one 
preceding it. A summary of the quantitative 
results obtained is presented in Table I, the 
average values from several runs being given. 
The widths of the Br and Rb edges in the 
compounds investigated, given in volts in column 
3 of Table I, were obtained from the experimental 
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raBLe I]. Physical constants of compounds used. 


CRYSTAL GRATING MELTING 
CoMPOUND STRUCTURI SPACE i Point 
KBr Cubic 6.57 730 
RbBr Cubic 6.86 683 
CdBr: Hexagonal 6.63 568 
NaBrQOs; Cubic 6.71 381 


curves by the triangle method! * and are not 
corrected for the effect of the crystals. This is 
not the true width of the edge’ but serves as a 
convenient means of recording and comparing 
results. The widths obtained, though slightly 
narrower, agree within experimental error with 
those quoted by Zinn." No pronounced fine 
structure in the edges themselves of the type 
found by Barnes" for elements of lower atomic 
number was recorded. 

The last five columns in Table I give the 
displacements of the secondary fine structure 
from the center of the main edge in volts, A 
being the first minimum in the transmission 
curve, a the first maximum, etc. Table II gives 
the physical characteristics of the substances 
used which are pertinent to the discussion. 

The most pronounced structure is shown by 
NaBrQOs, the substance of lowest melting point, 
in violation of the general rule that more 
pronounced structure is obtained for higher 
melting point substances. Others’ have noted 
such violations. RbBr with a comparatively 
high melting point showed the least pronounced 
structure though several different screen thick- 
nesses were tried. 

The structure for KBr checks satisfactorily 
with the first two maxima found in a previous 
investigation of a KBr single crystal.*:? RbBr 
should yield a similar structure slightly closer 
to the main edge, because of its larger grating 
space according to Kronig’s theory. Because of 
the unaccountably low intensity of the secondary 
structure for RbBr only the first maximum was 
obtained. The points a and B for RbBr are 
closer to the edge by an amount in accord with 
the theory. The fact that CdBr, though hex- 
agonal shows a close-in structure almost identical 





"P. A. Ross, Phys. Rev. 44, 977 (1933). 

" W.H. Zinn, Phys. Rev. 46, 659 (1934). 

* F. K. Richtmyer, S. K. Barnes and E. Ramberg, Phys. 
Rev. 46, 843 (1934). 

A. H. Barnes, Phys. Rev. 44, 141 (1933). 


with that of KBr and RbBr is interesting and 
perhaps explainable on the same basis as that 
used by Coster'® to account for similar structures 
exhibited by close-packed cubic and hexagonal 
lattices. NaBrQO, is cubic and of about the same 
grating space as RbBr and KBr; but as 3 Oatoms 
are associated with each Br atom in the lattice 
the close-in fine structure should not be similar 
to that of KBr and RbBr. Comparing we see 
that 8 and C differ for KBr and NaBrQOy but A, 
a, and B are quite closely the same. Hanawalt!® 
has obtained photographically pronounced sec- 
ondary fine structure for NaBrOs. A, a, and B 
obtained here agree with his data. 8 and C do 
not agree although it is possible that the large 
peak 8 is complex, being made up of two discon- 
tinuities close together as was evidenced by some 
of the runs, in which case the data are in better 
agreement. In an ionic crystal the outer electronic 
shells of Br are complete and after the ejected 
K electron reaches the outer shells the Br 
nucleus and incomplete A shell will act on the 
electron much as the element of atomic number 
one greater, namely, krypton. Now krypton has 
a large ionization potential of about 14 volts and 
it is probable that structure within 15 volts or 
so of the Br edge should be governed more by 
the atom!'? than by the crystal lattice, particu- 
larly since Kronig’s assumptions do not hold 
accurately close to the edge. This may be the 
main factor resulting in the similarity of A, a, 
and B in the cases studied. 

The comparison of the secondary fine structure 
of Br and Rb in RbBr is inconclusive because of 
the meager structure found for Br and the almost 
complete lack of structure for Rb, a and B here 
being doubtful. The structures might be expected 
to differ at least close to the edge.'* The two 
here are similar in that both exhibit compara- 
tively weak structure. a and B are not the same. 
No very extended structure was found which 
would prove whether the two are in better 
agreement 100 or more volts from the edge. 

It is of interest that a discontinuity due to 
the simultaneous ejection of a K and an L 
electron from a Br atom would have appeared 
at about 815 X.U. whereas the wave-length of 

'® D. Coster, Physica 2, 606 (1935). 

6 7. D. Hanawalt, Phys. Rev. 37, 715 (1931). 


‘7 W. Kossel, Zeits. f. Physik 1, 119 (1920). 
‘8 PD. Coster and R. Smoluchowski, Physica 2, 1 (1935) 
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the Rb edge is 814 X.U. It was necessary to 
search for such a discontinuity due to Br to 
make sure that it did not affect the runs on the 
Rb edge in RbBr. Two runs over the region 
810-820 X.U. with NaBrQ; as absorber failed to 
show any pronounced evidence for simultaneous 
K and L ejection large enough to affect the re- 
sults for the Rb edge. 

fine 
polycrystalline absorbers of 


The secondary structure obtained by 
others!" using 
atomic number greater than 34 and the double 


crystal spectrometer is not very pronounced and 
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is obtained only comparatively close to the main 
edge as in this investigation. More secondary 
fine structure of the Kronig type further from 
the edge than that obtained here undoubtedly 
exists for all the compounds but does not appear 
prominently enough to be recorded above experi- 
mental errors. It would be more pronounced if 
the absorption took place in a single crystal, 
the case of KBr®: 7 being one example. A study 
of Rb and Br in a single crystal of RbBr would 
prove interesting if a perfect enough crystal 
could be prepared. 
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A Simple Method for Demonstrating and Measuring Approximately the Index of 
Refraction of Crystalline Substances for X-Rays 


H. C. Hoyt anp Geo. A. Linpsay, University of Michigan 


Received January 21, 1936) 


An x-ray spectrum line may be formed by-reflection at 
the surface of a crystal even if that surface is not parallel to 
the atomic planes concerned in the reflection. If two faces 
are good very close to the line of their intersection, then, 
using a slit that is not too narrow, the same line may be 
reflected simultaneously from these two faces of a stationary 
crystal by turning it so that for the proper reflecting angle 
the rays must fall at the intersection. If the wave-length of 
the line is such that the reflected beam from the beveled 


N 1922 Davis and Terrill':? suggested a 

method of increasing the refraction of a 
beam of x-rays by grinding a crystal surface at 
an angle with the atomic planes so that the 
incident or emergent beam might make a very 
small angle with the reflecting surface. 

The present investigation retains the principle 
of the Davis method of unsymmetrical reflection, 
but introduces two important modifications. 
Firstly, two natural faces were used instead of a 
natural face and a ground surface; secondly, the 
same spectrum line was reflected simultaneously 
from two faces of the crystal by turning the 
crystal so that for the proper reflecting angle the 
rays fell at the intersection of the two faces. 
The crystal was stationary during the exposure. 


! Davis and Terrill, Proc. Nat. Acad. Sci. 8, 537 (1922) 
* Nardoff, Phys. Rev. 24, 143 (1924). 


face almost grazes the surface, then that portion will be 
refracted much more than the portion of the line reflected 
from the other face, and the two parts of the line are 
plainly separated. A very simple calculation gives the 
index of refraction from the width included between the 
two portions of the line. In this way 6/X* for a cerussite 
well with the 


crystal has been found to agree fairly 


calculated value. 


The wave-length of the line was such that the 
rays reflected from the almost 
grazed the surface, and were therefore refracted 
much more than the portion of the line reflected 
from the Bragg face, so that the two parts of 
the line were plainly separated. No reference 
line is needed but the two reflecting surfaces 
selected must be good very close to their line 


beveled face, 


of intersection. 


THE CRYSTAL STRUCTURE OF PbCO; 


The crystal specimens chosen were cerussite 
(PbCQO;) which in form and 
resembles aragonite in its structure and sym- 
metry.* The density is 6.5 g/cc. Sharp x-ray 
emission lines were obtained by reflection from 
the faces used. In a number of investigations 


is orthorhombic 


W. L. Bragg, Proc. Roy. Soc. A105, 16 (1924). 
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the refracting medium has been used as a grating 
to determine the index of refraction with crystals 
of low density. Thus Hjalmar* used gypsum 
specific gravity 2.3), and Larsson® used mica 
(sp. gr. 3.0), calcite (sp. gr. 2.7), and quartz 
sp. gr. 2.65). Cerussite was selected in the 
present work because its density is more than 
twice the density of the heaviest of these, and 
hence it should give a much greater refraction 
Since the method employed requires the use 
of two faces of the crystal for each exposure, a 
knowledge of the grating spaces and of the angle 
between the faces used is necessary. 
Measurements of the grating spaces for several 
faces of the crystal were made, together with 
two of the angles between faces. The results were 


as follows : 


Crystal plane (001) (O11) (110) (111) (102) 
d (x.u.) 3065 2483 4416 3583 2637 


The angle between the (100) and (110) faces 
was found by an optical spectrometer to be 
31° 21’ 30”, and that between the faces (111) 
and (102), 31° 6’ 10’’. From these values of d 
and one of the angles, the axes of the cry stal 
were calculated to be as follows: 


? 


a=5173+4 x.u.; b=8478+5; c=6129+33. 


The axial ratios from the above values are 
a:6:c=0.6102 :1 : 0.7229. 

These ratios agree well with those given in 

crystallographic tables, the latter values varying 

according to the specimen used for the measure- 


ments. For example, the following are three 
samples of previous values. 


a rhs 6S Source 
0.6102 : 1 : 0.7232 Kokscharow 
0.6100 : 1 : 0.7230 Goldschmidt's Atlas 
0.6101 : 1 : 0.7230 Artini 


Our probable errors are based on the assump- 
tion that the above values of d are accurate to 
about 1 x.u. 


lhe axial values given by Zachariasen® are: 


a=5.14A, b=8.45A, c=6.10A. 





*Hjalmar, Dissertation, Lund (1923). 

a Larsson, Dissertation, Uppsala (1929). 

W. H. Zachariasen, Norske Videnskaps-Akad., Oslo. 
Mat.-Naturv. Klasse (1928). 
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l—u o=tan (0—a) sine G— a je 


Fic. 1. Geometrical relations involved in the derivation for 
the expression of the index of refraction 


These latter values were determined by the 
powder method. 


APPARATUS 


All spectra in this investigation were obtained 
with a Siegbahn vacuum spectrograph. The 
radius of the spectrograph was 183 mm, and the 
width of the slit was 0.07 mm. Cellophane was 
used to seal the slit and visible light from the 
filament or other sources was intercepted by red 
colored gold beater’'s skin. The materials, for the 
several characteristic lines studied, were either 
soldered to the target or rubbed on as a powder. 
The crystal was mounted with the line of 
intersection between the two faces in the axis 
of rotation of the crystal table. 


(GEOMETRICAL RELATIONS AND DERIVATION OI 
THE EXPRESSION FOR THE INDEX 
OF REFRACTION 


In Fig. 1, AC is the slit width very much 
magnified and AEDJC is the outline of the beam 
directed towards the intersection of the two 
faces. FG is the portion of the spectrum line 
from the Bragg face and J// the portion from 
the beveled face. ¢ represents the angle of 
separation, measured from the inside edges, 
between the two portions of the line. 

The path of the incident and refracted rays 
ACED is shown in Fig. 2. 

From the definition of the index of refraction 


sini cosp 
= = =1—4, (1) 
sinr cosy 


HOYT AND 














Fic. 2. Path of the incident and refracted ravs 


from which 


cosy—cosB8 2 sin }(y+B) sin }(B—y) 
= = . (2) 
cos ¥ COs ¥ 
8—y is the deviation of the ray ED due to 


refraction. 8 and y are nearly equal, therefore 


(y+8)/2 equals 8 or y approximately (until 8 


becomes very small). The sines and tangents 
are approximately equal to the angles expressed 
in radian measure. 


Therefore 


2 sin 8 sin 3(8—y) 
= = 2(B—y) 


cos ¥ 


very approximately. (3) 


From Fig. 1, 8B=@—a and (B—y)=e, 


so that 6=€(@—a). (5) 


Thus 6 is simply the product of two angles 
expressed in radian measure. 

Since the separation which was measured is 
the difference between the pe sitions of the inside 
edges of the two portions of the separated 
spectrum line, a correction is necessary for the 
refraction of the portion of the line coming from 
the Bragg face. The deviation due to refraction 
of the ray reflected at the Bragg face is given 
by 0’—@=4cot 6 where @ is the glancing angle 
within the crystal and 6’ is the corresponding 
angle between the emergent ray and the Bragg 


face. Hence the amount to be added to e is 

6 cot 6 

or 6=€(0—a)+6(0—a) cot 8, (6) 
€(06— a) 

whence 6= (7) 


1—(#@—a) cot 0 
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Fic. 3. Plot of 6/X° against A 


The quantum theory of dispersion as developed 
by Kronig, Kramers,’ Kallman and Mark,’ and 
Prins® gives the following relation between 6 and 
the frequency 


” Nj vii? v 

5= » i 1+{—) Inj1i- 
2rc?2m v? | y v; 
where N; is the number of electrons of frequency 
v; per cubic centimeter, and » is the frequency of 
the incident radiation. Damping is neglected. 

Rewriting the equation for 6, after substituting 
1/ for v and dividing by \* we have 


where N is the total number of electrons of the 
medium per cubic centimeter. 

5/d*, calculated from this equation is plotted 
with A giving the dispersion curve shown in 
Fig. 3. 

A few comments may be made on this method 
of determining the index of refraction for x-rays 


a) The slit must not be too narrow. 

(b) For good lines an intense radiation should strike 
both faces. This is especially true of the beveled face since 
the radiation reflected from the beveled face traverses a 
considerable amount of the crystal and is therefore strong] 
absorbed. 

c) The direction of radiation cannot be reversed, that is, 
it must take the direction indicated in Fig. 1. If it were 
sent towards the crystal in the reverse direction, there 
would be no separation of the line, since, for a given wave- 
length, the reflected rays would emerge practically parallel 
from both faces. 

d) The limiting wave-length for a given crystal can be 
determined by substituting a, Fig. 1, for @ in the Bragg 


H. A. Kramers, Nature 113, 673 (1924). 
8’ Kallman and Mark, Ann. d. Physik (4) 82, 585 
J. A. Prins, Zeits. f. Physik 47, 479 (1928). 
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Fic. 4. Separation of a monochromatic beam due to 
refraction. 


relation m\=2d sin @. The minimum wave-length which 
may practically be employed is somewhat larger than this, 
since the refracted beam is nearly all absorbed if @—a is 
very small. 

e) As stated previously an imperfect edge between the 
two planes is likely to cause serious error in the readings 
For one crystal the character of the intersection between 
the (111) and (102) planes was studied by taking photo- 
graphs at shorter distances. The regular distance from the 
axis of the spectrograph to the photographic plate was 
183 mm. One photograph was taken at 103.5 mm and 
another at 41.4 mm. This procedure was repeated for 
different lines. Since the photograph is a record of two 
reflected beams of rays having the same wave-length, any 
separation due to an imperfect edge would be the same for 
all distances of the plate. With these photographs the 
separation was proportional to the distance, indicating that 
the separation was due to an actual difference in direction 
between the two beams and that any irregularities in the 
surfaces close to the edge were too small to affect the 
readings. 

A microscopic examination of the second crystal gave no 


indication of imperfections near the edge. 


RESULTS 


The edge of the line to be measured at G 
(Fig. 1) was in all cases quite sharp, while that 
at H was less so, especially when the angle 
between the face and the emergent ray was 
small. Since the incident rays are concentrated 
into a narrower beam after reflection from the 
beveled face, as shown in Fig. 1, it might be 
expected that the refracted line would be rather 
intense but, as explained in (d), when (@—a) is 
small the refracted line is strongly absorbed. 


Observations were made over a range of 


raB_e I. Summary of computed values of 5/d* for two crystals 
and the various wave-lengths 


LIN! A(x i(x.u mit Simm € SR 2-10 


Crystal No. 1, a=31°7 


Cu KB;(2) 1389 2637 41 0.195 0.001065 6.71 
W La;(2) 1473 2637 171 .053 000289 7.17 
Fe Ka,(2) 1932 3583 91 161 000781 6.52 
Sn Ly; 2995 2637 210 153 000836 6.25 
Sb Lg 3218 2637 390 .0S2 .000503 6.45 
Sb La; 3432 2637 570 .067 000366 6.39 
Sn La; 3592 2637 711 .057 .000311 6.42 
Ag LA, 3927 3583 127 393 (002147 5.45 
Ag La, 4146 3583 254 217 001186 5.69 
Rh LB, 4364 3583 384 .137 000748 5.14 
Rh La; 4588 2637 1764 .029 000158 5.44 
Crystal No. 2, a=31° 23.5 
Ru Lae 4836 4416 108.5 0.469 0.002562 3.63 
SAB, 5021 4416 195.5 308 001667 4.14 
Mo Lg 5166 4416 264.3 301 001644 5.30 


wave-lengths extending from 1389 x.u. to 5166 
x.u. Several photographs of the separated lines 
are shown in Fig. 4. Table I contains a summary 
of the lines with their wave-lengths, the grating 
constants used, the angle (@—a) between the 
reflected ray and the beveled face, the angle « 
due to the difference in refraction between the 
two portions of the beam, and the final value of 
6/* as calculated from Eq. (11). The readings 
for the last three wave-lengths appearing in 
Table I were taken with the second crystal. All 
observations were taken in the first order with 
the exception of the lines Fe Ka,;, W La, and 
Cu K8,, which were taken in the second order. 

For each line the difference S is the mean of 
the measured separations on two or more photo- 
graphic plates. Ten settings of the comparator 
were made on each position. 

Fig. 3 contains the theoretical curve for 6/X 
and also the experimental values represented by 
small circles. 

The experimental points obtained do not show 
conclusively the influence of the absorption 
edges. However, in the region of the smallest 
value of 6/X*, that is, near the My absorption 
edge, the two points for Ru La; and S KB, fell to 
low values, and they lie nearly on the curve. 

The results, though approximate (probably 
accurate to about ten percent), are simply and 
directly obtained without special apparatus. 
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Excitation Potential, Relative Intensities and Wave-Lengths of the Ka’’ X-Ray 
Satellite Line 


LyMAN G. PARRATT,* Cornell University 


Received January 25, 1936) 


two-crystal vacuum spectrometer, ionization 


With a 
curves of the Ka” 
elements S(16) to V(23). This satellite line, on the short 


satellite line have been recorded for 


wave-length side of and very close to the Ka line, was 
observed in 1922 by DolejSek and has been until now a 
forgotten line. The excitation potential of this satellite has 


been determined: The initial atomic state for its emission 


I. INTRODUCTION 


ATELLITES comprise a large fraction of the 
total number of characteristic x-ray emission 
lines; but as yet we do not have a generally 
acceptable explanation of satellite origins. The 
greatest need in resolving this enigma is that of a 
more complete and accurate experimental de- 
scription of these nondiagram lines. 

In another paper’ the writer has recently 
reported the results of experiments on the 
excitation potential of the satellite lines of the 
Kas, 4 group. This group of lines consists (for 
elements near titanium (Z=22)) of 
‘and ay. For these lines? 


five com- 
ponents: a’, a3, as’, a3’ 
the excitation potential was found to be in agree- 
ment with the voltage required to ionize doubly 
the emitting atom by ejection of a K and an L 
electron by a single cathode-ray impact. This 
result is in support of the Wentzel-Druyvesteyn 
theory which supposes that the atomic (or ionic) 
transition KL—LL gives rise to the Kaz, , lines. 
The present paper is concerned with determina- 
tions of the excitation potential, relative in- 
tensities and wave-lengths of the Ka” satellite 
line. We shall conclude that this line is not of 
the a3, 4 group but, in terms of the Wentzel- 
Druyvesteyn theory, is produced by a K /—-LM 
transition. 


* The author is indebted to the Carnegie Foundation for 
a grant-in-aid of this research. 

‘L. G. Parratt, Phys. Rev. 49, 132 (1936). 

* The intensity of the two strongest lines, a; and ay, 
dominates the intensity of the group, and, as discussed in 
reference 1, whether or not all five components of the group 
have the same excitation potential could not be deter- 
mined in these experiments. 


appears to be one of KM ionization. The Ka” intensity 
relative to the Ka intensity varies with atomic number 
reaching a maximum of 2.3 percent at Ca(20). The Ka 
wave-length positions have also been determined for each 
element. The observed data are compatible with both the 
Wentzel-Druyvesteyn and the Richtmyer theories for the 


origin of this satellite line. 


Il. EXPERIMENTAL PART 


The 


which were mounted calcite crystals A-By, was 


two-crystal vacuum spectrometer, on 
used in the present as in the previous experi- 
ments.' The advantages of this type of instru- 
ment over the single-crystal spectrograph are (1 
in the greater practical resolving power and 
dispersion, and (2) in the greater accuracy with 
which intensities can be measured. The relative 
the method of 


disadvantages of ionization 


intensities as compared with the 
photographic method the 
difficulties encountered in preparing satisfactory 


2) in 


registering 
are (1) in greater 
targets of the materials to be studied, and 
maintaining the higher vacuum in the x-ray 
tube necessary for accurate intensity measure- 
ments. These advantages and disadvantages 
have been mentioned elsewhere.!: * The reader is 
referred to previous papers for discussions of the 
which the 
and of 


present ionization 


the 


accuracy! with 


curves' were recorded instrument 
perturbations® which distort the true x-ray line 
contours into the observed curves. 

The targets used in these experiments were 
prepared as follows: CaCle, KCI, SrCle, K2PQ,, 
CaO and Sc,O; in the powdered form were 
pounded into a roughened copper surface. CuC] 
was formed by passing chlorine gas over copper. 
A metallic calcium made in the 
manner previously described.* The pure metallic 


target was 


titanium and vanadium targets have also been 


L. G. Parratt, Rev. Sci. Inst. 6, 372 (1935). 
‘The experimental set-up in the present work was 
identical with that in the previous work, reference 1. 
>L. G. Parratt, Rev. Sci. Inst. 6, 387 (1935). 
‘LL. G. Parratt, Phys. Rev. 44, 695 (1933). 
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Fic. 1. lonization curve of the Ca Ka region. This curve was recorded with a CaO target and with a tube potential of 
20 kv. The as, « group of satellites is plotted to an intensity scale 68 times the scale used in plotting the a’’, a and ay 
lines. The satellite background (the side of the a, line) has been sketched in and the satellite contours, corrected for this 


background, are also drawn in the figure. 


described,*: * and likewise the sulphur targets.’ 
All of these targets were satisfactory in studying 
the Ka;, » doublet and the Ka” lines, but only 
FeS, K2ePO,;, CaO, ScsO3, questionably CaCle, 
and metallic Ca, Ti and V served satisfactorily in 
obtaining data on the Kaz , lines. This latter 
group of targets would dissipate, with the focal 
spot approximately 2 mm in diameter, from 400 
to 700 watts of power. 


Notation 


The term Ka’’ in the present paper is used to 
refer to the ‘“‘bump” which appears on the short 
wave-length side, near the base, of the Ka, line 
with elements S(16) to V(23). That the reader 
may gain perspective, an ionization curve of the 
Ka region, Ka;, », a’ and the as. , group, for 
calcium (20) is shown in Fig. 1. Ca Ka’’ appears 
at the wave-length 3349.1 X. U. In Fig. 1 of the 
previous paper' the same Ka region is reproduced 
for titanium (22). 

The notation used to designate x-ray satellites 
is confused; this confusion becomes more ap- 
parent as our knowledge of these lines becomes 
more complete. Many faint x-ray lines known in 
the literature as satellite or nondiagram lines can 
now be accounted for easily as diagram lines for 


"L. G. Parratt, Phys. Rev. 49, 14 (1936). 


which the quadrupole, rather than the dipole, 
selection rules are obeyed. Another type of con- 
fusion is illustrated by the following: Hulubei 
recently reported® a triplet nature of the Kaz, ; 
satellites in the region of elements Cu(29) to 
As(33). The third component of this group, the 
one between a; and a;, Hulubei calls a;’. This is 
an extrapolation of the notation of Backlin®: '® 
from the range of elements Si(14) to Cl(17). The 
writer, using the two-crystal spectrometer, has 
been able to trace" the component a;’ from S(16) 
to Ge(32). Relative to the a; or to the a, line, a;’ 
moves progressively as the atomic number 
changes: In the region of low atomic numbers the 
a;’ is in the position reported by Backlin® and by 
Ford,'® that is, on the long” wave-length side of 





§ M. H. Hulubei, Comptes rendus 201, 544 (1935). 
*E. Backlin, Zeits. f. Physik 33, 547 (1925); 38, 215 
(1926). 

1°Q. R. Ford, Phys. Rev. 41, 577 (1932). 

't Discussed in more detail in another paper in process 
of publication. 

2 1f these trends are extrapolated to atomic numbers 
lower than Si(14) the a; and a;’ lines should be coincident 
in wave-length at Al(13) or Mg(12): In fact, with these 
elements, Al and Mg, no a;’ has been found while the rela- 
tive intensity of a; has been reported as being anomalously 
large (reference 10). From ionization curves at S(16) the 
writer estimates the intensity of the resolved a;’ compo- 
nent to be approximately 80 percent of the intensity of the 
a; component. With elements of atomic number less than 
12 the a;’ line may be expected to appear, if one continues 
the extrapolation, on the long wave-length side of az. 
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a,;; at Ti(22) the a; and a, are coincident in 
wave-length; and with elements of atomic num- 
ber greater than 22 a;’ is on the short wave-length 
side of as. These observations, in agreement 
with those of Hulubei, show that in the range of 
higher atomic numbers az, ; is complex, actually 
a quadruplet since the a’ line, although very 
faint, is still present, and also show that the 
terms a;’ and a, as assigned by Hulubei should be 
interchanged. 

A third type of confusion in notation is in re- 
gards to the present line Ka’. In 1922 DolejSek"™ 
observed an “emission band”’ on the short wave- 
length side of the Ka, line with elements Cl, Kk, 
Ca and Sc. The short wave-length edge of this 
band he called a’; the known as a’ 
Dolejsek, following Siegbahn’s lead, called az. 
Dolejsek’s band” and the writer’s 
“bump,” a’ 
It speaks well of DolejSek’s care in his observa- 


tions and of the resolving-power of his spectro- 


line now 


“emission 
’, are the same line (or group of lines). 


graph that he was able to detect this emission 
band: Of the numerous investigators following 
him, none, so far as the writer is aware, has re- 
ported finding it. 

Druyvesteyn"™ refers to DolejSek’s emission 
edge a’ as the a” line, and reverts to Hjalmar’s'® 
use of the term a’. In Fig. 1 the a’ line is at 3340.1 
X. U. Deodhar" has reported a possible doublet 
nature of the Ka’ line with elements Si(14) and 
P(15). The terms a’ and a”’ were used by Deodhar 
to refer to this doublet. So far as the writer is 
doublet 


observer has seen a 


or other elements, and 


aware no other 
structure in a’ for Si, P 
the writer has found no indication of such struc- 
ture for S(16). Deodhar states that his a’’ may 
have been due to chemical combination or may 
have been spurious. We shall therefore adhere to 
the notation used by Druyvesteyn in assigning 
the term e@” to the “bump.” 


Relative intensities 


lonization curves of Ka’’ with elements S(16) 


to Cr(24) are reproduced in Fig. 2. The a”’ line 
is most prominent at Ca(20). This may or may 
not mean that the intensity of a” relative to the 
intensity of a; is a maximum at Ca: We do not 


M. V. DolejSek, Comptes rendus 174, 441 (1922). 
'*M. J. Druyvesteyn, Doctor's Dissertation, Groningen, 
1928. 
'* E. Hjalmar, Zeits. f. Physik 1, 439 (1920). 


‘6G. B. Deodhar, Proc. Roy. Soc. A131, 633 (1931). 


YMAN G. 


PARRATT 


know the shape of the a, line and consequently 
we have considerable latitude in drawing in the 
background of the a” satellite. In Fig. 2 the side 
of the a) line for each element has been draw nin 
such that the width at half-maximum intensity 
of the a” line is approximately the same as the 
width of the a, line. There is no justification for 
choosing this resolution; in fact, 
satellite lines in general are wider than the associ- 
The component lines of the 


particular 


ated diagram lines: 
Ca Kaz; , satellites are 1.5 to 3 times as wide as 
the a, line. Furthermore if a@”’ is in reality a group 
of lines the width should certainly be greater 
than is here assumed. This question bears not 
only on the shape and possible multiplet struc- 
ture of a” but also on the asymmetry of the a 
lines, a subject beyond the scope of the present 
paper." In Fig. 3 is plotted the intensity of a” 
relative to the intensity of a; when the resolution 
is arbitrarily made as indicated in Fig. 2. These 
intensity ratios are listed in Table I. The un- 
certainty in the relative intensity due to the 
indefinite background is greatest for S(16) and, 
with this element, may be as great as a factor of 
two or even three. The values given in Table | 
are probably a minimum. A slight trace of the 
existence of a”’ for Cr(24) was observed but this 
evidence was within the observational error. The 
reality of the presence of a’ for V(23) may be 
questioned since the ‘“‘bump”’ with this element is 
close to the observational error. 

TABLE I. Relative intensitves and wave-lengths of the Ka 
satellite line. Four different chemical compounds of Cl(17) 
CaCle, KCI, CuCl and SrCle, and two chemical compounds 
of Ca(20), CaCle and CaO, in addition to metallic Ca 
were used as targets. No difference in the a” line due to 
chemical binding was observed greater than the observa- 
tional error. In the case of CaO and metallic Ca, however 
a marked difference is observed in the widths of the a 
lines. See Figs. 1 and 4. The wave-lengths given by 
DolejSek"* for the short wave-length edge of his “emission 


band” are listed in the last column. DolejSek did not give 
wave-lengths of the Ka, line. 


RELATIVI Wave-LENGTHS 
INTENSITY a 
TARGET ELEMENT (a’’/ay) X 100 ai* a’ Dol k 
FeS S(16) 1.0 5361.28 5358.6 
CaCl. Cl(17) 1.4 4718.2 4714.8 4712 
KPO, K(19 2.1 3733.68 3730.9 3730 
CaO Ca(20) 2.25 3351.69 3349.1 3349 
Sc.O;  ‘Sc(21) 1.4 3025.03 3022.8 3023 
metal Ti(22) 0.9 2742.87 2740.9 
metal V(23) 0.3 2498.42 2496.5 


* The wave-lengths of the Aa line of vanadium and of titamun 
ken from Bearden and Shaw, Phys. Rev. 48, 18 (1935); of s 
from Siegbahn, Spektroskopie der Réntgenstrahien 


43. 612 (1933 


ta 





to chlorine 
and of sulphur from Valasek, Phys. Rev 
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Fic. 2. lonization curves of the Ke” line for elements $(16) to Cr(24). The tube potential used in 
recording these curves was 15 to 20 kv. The abscissae and ordinate scales are arbitrary in the case of 


each curve. 
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Fic. 3. Intensity of the Ka” satellite line relative to the intensity of the Ka, line for elements 
S(16) to Cr(24). As mentioned in the text, these values are very sensitive to slight changes in the 
arbitrarily assumed shapes of the a lines: Whether or not this curve reaches a maximum point 
depends upon the shape one chooses to draw in each of the curves of Fig. 2. 





Wave-lengths 


Measured from the peak of the Ka, line as a 
standard, the wave-length positions of the peaks 
of the a” lines, resolved as discussed above, are 
as given in Table I. Only differences in wave- 
length can be accurately measured with the 
writer's spectrometer, so the wave-lengths of the 
a, lines used as references are included in the 
table. Surprisingly good agreement exists be- 
tween the present a’’ wave-lengths and those 
given by DolejSek for the edge of his ‘emission 
band."’ 


Excitation potential 


As discussed in the previous paper! a knowl- 
edge of the minimum voltage required to produce 
a given satellite line is especially requisite in 


developing or in testing any theory of satellite 
origin. In the case of the Ka” satellite, however, 
the two principal contemporary theories, the 
Wentzel-Druyvesteyn theory" and Richtmyer's 
“double-jump” theory,'? predict the same excita- 
tion potential: The initial state of the atom 
emitting the Ka” line is supposed in each case to 
be that of KM ionization. 

Fig. 4 shows ionization curves of Ca Ka” 
recorded at various voltages. For these curves a 
consistency in the background is managed by 
superimposing the observed curves, matched in 
intensity at the a; peak, on an assumed but fixed 
a, shape. Consequently, the differences in the 
relative intensity which may be present with 


F. K. Richtmyer, J. Frank. Inst. 208, 325 (1929). 
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Fic. 4. lonization curves of Ca Ka” recorded at various tube voltages. A metallic calcium target was used. No relative 
uncertainty is present in the background of these curves as discussed in the text. 

Although there is a very pronounced difference in the widths of the Ca Ka, ¢ lines observed from targets of metallic 
calcium and of CaO (Fig. 1), no difference in the a” satellite was observed greater than the observational error. 


’ is a satellite line. There can be no doubt that 





different tube voltages or tube currents can be a 


measured without the background uncertainty 
mentioned above. In testing for a possible effect 
of chemical binding by using targets of calcium 
oxide and of metallic calcium the background 
uncertainty is again large because of the different 
a, shapes. See the widths of a; in Figs. 1 and 4. 


the excitation voltage is less than that which is 
required to produce a state of KL ionization. This 


the conclusions of the 


indicates, in view of 
previous experiments,' that @”’ is not of the a, 
group of satellites. 

If one assumes that the curve of Fig. 5 should 


be similar in shape to the curve of the intensity of 


The a” relative intensity was observed to be 
constant with tube current from 10 to 30 ma at az, 4 vs. voltage,' one obtains, by the consequent 
15 kv. extrapolation to zero intensity, the excitation 


The a” relative intensity as a function of tube 
voltage is plotted in Fig. 5 and the data are 
listed in Table II. As the voltage decreases the 


voltage of a’’ as 40704250 volts. The energy 
required to remove a K electron of atom Z=20 
(calcium) is 4030 volts.'* The energy required to 














curve of this figure definitely breaks at about 2 
twice the excitaton potential and exhibits the 2 : 
° ° - . ” . . 
same general features as the voltage function for 20h / 
3 : > . . - , 9 
the intensity of the Kaz, 4 satellites, Fig. 2 of / 
reference 1. This may be taken as evidence that a 7% 
<x | | 
| 
TABLE II. The intens‘ty of the Ca Ka’’ satellite line rela 
tive to the intensity of the Ka, line as a function of the x-ray & 
tube voltage. =o. | 
RELATIVE i — : 
CURRENT INTENSITY i. i i s ae) 35 5) "26 2S 
VOLTAGI i? ma a’ ai) X 100 K KL Voltage (KV) Ye ws : si 
500 115 0.80 . ~ , ° ‘ = , 
eo 89 1.50 Fic. 5. Intensity of the Ca Ka” satellite line relative to 
3500 88 190 the intensity of the Ka; line as a function of the tube 
6500 65 15 voltage. The voltage for which the intensity ratio a” /a 
8000 57 30 falls to zero is the satellite excitation potential. The 
10 000 rT, 218 voltage indicated by K, 4030 volts, is the excitation po- 
; ) 2.18 yrege “dates “lbs KL. 4440 
15 000 7 > 20 tential of the Ka; 2 lines; that indicated by KZ, 
0000 15 7 0) volts, is the excitation potential of the Kaz, , satellites. 
? 725 -ae Cc. 
30,000 12 2.25 * M. Siegbahn, Spektroskopie der Réntgenstrahlen, second 


edition (Julius Springer, 1931), p. 348. 
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remove an Jy, yp electron of atom Z=21 


(scandium) is approximately 37 volts. Then the 


energy requisite to producing a state of K.V/ 
ionization in the calcium atom is 4067 volts, in 
excellent agreement with the experimental value. 


III. CONCLUSIONS 


From the excitation potential measurement we 
conclude that the initial state of the atom 
emitting the Ka” satellite is one of K.V ioniza- 
tion, a conclusion which is in agreement with the 
Wentzel-Druyvesteyn theory and also with the 
Richtmyer theory for the origin of this satellite. 

From the wave-length positions there 1s 
no choice between the theories. Druyvesteyn, 
using DolejSek’s wave-length measurements, has 
shown" that the Av interval between a” and a, 
is of the same order of magnitude as the interval 
between the frequencies of the WZ levels of atom 
Z+1and of atom Z. On the basis of this evidence, 
Druyvesteyn ascribed, in accord with the 


APRIL 1, 1936 PHYSICAL 


Wentzel-Druyvesteyn theory, the origin of the a” 
line as due to the transition K.1J—LA1. In this 
region of atomic numbers the energy values of the 
M shells are rather indefinite’ and by “‘order of 
magnitude’ must be meant a factor of possibly 
several-fold. Agreement of about the same order 
of magnitude is also obtained with Richtmyer’s 
theory. 

From the relative intensities, so far as theo- 
retical predictions can be or have been made,'*: °° 
there is likewise, as yet, no choice between the 
two theories. 
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Effect of Intense Illumination on Time Lag in Static Spark Breakdown 


Harry J. Wuite,* University of California, Berkeley 


(Received February 5, 1936) 


The time lag of spark breakdown in static uniform fields 
and with intense illumination of the cathode has been 
studied in air, helium and carbon dioxide. An electro- 
optical shutter was used to observe the time lag, while an 
auxiliary spark gap supplied the intense cathode illumi- 
nation. In air, the time lag was about 10-7 second for 
overvoltages of a few percent and increased very rapidly 
with decreasing overvoltage. Increasing the overvoltage 
to above 30 or 40 percent reduced the time lag to a more or 
less constant value of 2 or 3 10~* second. This, in part, is 


INTRODUCTION 
HE time lag of spark breakdown in static 
fields has been studied in a thorough and 
systematic manner by Zuber' and Tilles? for lags 
from about 1 second to below 10-* second. The 


* Now with Research Corporation, New York. 
'K. Zuber, Ann. d. Physik 76, 231 (1925). 
* A. Tilles, Phys. Rev. 46, 1015 (1934). 


shown to be due to the nature of the initiatory spark. The 
position of the midgap streamer observed in previous 
experiments has been found to depend on both the over 

voltage and intensity of illumination. An explanation of 
these observations is advanced on the basis of space charge 
effects. The results in carbon dioxide were quite similar to 
those in air. In helium a much higher overvoltage was 
found necessary to produce a given time lag than in air. 
This is explained by considerations of the relative rates of 
gain of energy by electrons in the two gases. 


effects of overvoltage and intensity of illumi- 
nation of the cathode on the time lag were 
investigated. Time lags down to about 210-7 
second with low overvoltage and high illumi- 
nation intensity have been observed by Snoddy.* 
Very short time lags of the order of magnitude of 


*L. Snoddy, Phys. Rev. 40, 409 (1932). 
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10-8 second 
Tam, 
voltages. None of the investigations of the very 


have been observed by 


Beams* and others, 


short lags has been quantitative, correlating the 
and intensity of illumi- 


effects of overvoltage 
nation on the time lag. To this end the Kerr cell 
electro-optical shutter has been used to study 
these short lags and the results are being reported 


in this paper. 


METHOD AND APPARATUS 


In the method used, the time lag of breakdown 
of a given gap is determined by measuring the 
time interval between the impulse which sets off 
the spark and the appearance of luminosity in 
the gap. The starting impulse is provided by a 
flash of the 
cathode of the gap and coming from a spark in an 


ultraviolet radiation focused on 
auxiliary gap. The elapsed time between the 
starting impulse and the appearance of lumi- 
nosity is determined by using an electro-optical 
shutter which is focused on the gap being studied, 
but controlled by the potential across the 
auxiliary gap. Thus, the shutter closes at a 
definite time after a spark occurs in the auxiliary 
gap, this time until 


luminosity is just observed to appear, the time 


and, by varying open 


lag is determined.’: *: *: '° 


A schematic diagram of the apparatus is shown 
in Fig. 1. Gz is the gap for which the time lag is to 
be determined and G, is the auxiliary gap which 
illuminates Gs and controls the closing time of the 
electro-optical shutter. The high voltage direct 


potentials are supplied by transformers with 


thermionic tube rectifiers and smoothing con- 
densers of sufficient capacity to reduce the a.c. 


ripple voltages to negligible values. The electro- 


*P. O. Pedersen, Ann. d. Physik 71, 317 (1923). 

®>R. Tam, Archiv. f. Elektrotechnik 19, 235 (1928). 

6]. W. Beams, J. Frank. Inst. 206, 809 (1928). 

7 This assumes that the photoelectrons from the cathode 
are released instantaneously, and that the appearance of 
luminosity is simultaneous with the breakdown of the gap 
as marked by the rapid fall of the potential difference 
across it. The photoelectric effect has been shown to be 
instantaneous to within 3X10~* second by the work of 
Lawrence and Beams,* while observations by Dunnington 
and the writer'® show that luminosity first appears within 
about 2X10~° second after the beginning of breakdown. 

*E. O. Lawrence and J. W. Beams, Phys. Rev. 32, 478 


’ 


(1928). 
*F. G. Dunnington, Phys. Rev. 38, 1535 
'H. J. White, Phys. Rev. 46, 99 


1931). 
1934). 


Peders« mn, 


by using large over- 
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Fic. 1. Schematic diagram of apparatus. 


optical shutter is designated by E.O.S. and 
involves no unusual features. 

The voltage supply for gap G2 is divided into 
two sections for the purpose of applying an 
overvoltage of a definite and easily variable 
value. The upper section supplies a_ voltage 
to 


voltage of Ge, while the lower section supplies the 


approximately equal the static sparking 
overvoltage. The magnitude of the overvoltage is 
measured by means of a voltmeter consisting of 
an accurate 20 megohm-resistor and a _ milli- 
ammeter. The overvoltage is applied to Gz by 
moving the switch SW2 from the right to the left 
contact. Gap G, is fired by closing switch SW1 
which raises the voltage to a few percent above 
the sparking value. A quartz mercury arc is used 
to illuminate the cathode of G, in order to reduce 
the time lag to a small value. The switches are 
synchronized so that SW2 closes a few hundredths 
of a second ahead of SW1, thus enabling an 
essentially static voltage to be applied to gap 
Gz without its sparking before being set off by 
gap G,. To be sure that the switching operation 
did not affect the voltage supplied to Go, tests 
were made for its sparking with zero overvoltage 

The procedure followed in making the measure- 
ments is as follows. The length of the leads to the 
Kerr cell is adjusted to a definite value and the 
voltage of the upper section of the supply for G: 
set at just the minimum sparking potential. Then 
gap G, is fired at intervals of about 10 seconds 
and the overvoltage on G2 raised until luminosit) 
is just seen to appear. Using this procedure the 
time lag was measured as a function of the 
overvoltage and intensity of illumination for 
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several gap lengths in air at atmospheric pressure 


and also for helium and carbon dioxide. 


DATA 


The results obtained for air at atmospheric 
pressure and for gap lengths of 1 mm, 3 mm and 
5 mm are shown in Fig. 2. The results for helium 
at a gap of 10 mm and for carbon dioxide at a 
gap of 5 mm are shown in Figs. 3 and 4, re- 
spectively. The curves for the three gap lengths 
in air are very similar. At overvoltages of several 
percent the time lags are of the order of magnitude 
of 10-7 second. Increasing the overvoltage causes 
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Fic. 2. Time lag curves for air. 


the time lags to decrease rapidly at first, but to 
finally turn up sharply and approach a nearly 
constant value of 2 or 3 10> second. 

The effect of varying the cathode illumination 
intensity was noted by interposing wire screens 
between the two gaps which reduces the effective 
aperture of the quartz lens L (Fig. 1). Decreasing 
the intensity resulted in increasing the time lag, 
the effect being much larger at low overvoltages. 
Thus, referring to Fig. 2 and the 3-mm gap, 
decreasing the intensity from 1 to 0.002 at 10 
percent overvoltage increased the time lag from 
68 X 10~* second to 130 10~-" second, while at 30 
percent overvoltage the corresponding increase 
was from 37X10-* second to 5310~-° second. 
Experimental difficulties prevented measure- 
ments being taken at higher values of overvoltage. 
However, from the trend of the curves it appears 
that decreasing the intensity at high overvoltages 
increases the time lag by a more or less constant 
amount. 





uy 100 i 
© o| 
\ 
s 
~) 80 + X 
Ss 
e \ 
~ ' 

60 | ans eal 
3 ’ r jee * 7 

| a | n $ 

K ~ 
> 40 | — } 
w 100°le 7 me 
UO | 
x 20 + + + 4 + + + 
Q 

0 20 40 #460 80 100 —f20—Ss«sf40 





TIME LAG-SECONDS«/09 


Fic. 3. Time lag curves for helium, 5-mm gap 
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Fic. 4. Time lag curves for carbon dioxide, 5-mm gap 
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The measurements show that the time lag is a 
function of overvoltage and intensity of illumi- 
nation ; that the order of 
magnitude of 10-7 second occur in gaps of a few 
mm at low overvoltages but only provided very 


also time lags of 


intense cathode illumination is used. Limitations 
of the apparatus prevented time lags of longer 
than 1.510~-7 second being measured, so that 
the curves could not be traced down to very 
small overvoltages. However, from the trends of 
the curves at low overvoltages, the time iag for 
voltages just above the minimum sparking 
potential'' must be orders of magnitude larger 
than 10-7 second, and may possibly run con- 
tinuously into the formative lags observed by 
Tilles.2 It is of interest to extend the present 
measurements into the time lags 
observed by Tilles, and work in this direction is 


R. R. Wilson at the 


region of 


being undertaken by 
University of California. 

From Fig. 3 it is seen that much higher over- 
voltages are required in helium than in air or 
carbon dioxide to produce approximately the 
same time lags. The significance of the curves for 
helium discussed in the theoretical 
interpretation of the No. particular 
significance can be attached to the curves for 
carbon dioxide which are more or less analogous 


will be 
results. 


to those for air. 

In the course of the present study it was found 
that the position of the midgap streamer’: '” 
which appears for long gaps is affected by both 
the overvoltage and the intensity of illumination 
of the cathode. Fig. 5 shows approximately how 
the location of this streamer varies with over- 
voltage for a 5-mm gap. As the gap voltage is 
raised above the minimum sparking potential, 
the streamer moves toward the anode and reaches 
the latter at an overvoltage of about 10 percent. 
Further increases in overvoltage cause the 
streamer to move back toward the midgap 
region. The position of the streamer for this 
latter range of overvoltage tended to be quite 
variable and at times even two parallel streamers 


“ The overvoltages are calculated relative to the static 
sparking potential obtained with low intensity illumination, 
i.e., in terms of the ratio of the observed potential to the 
so-called normal sparking potential of the gap. Actually 
the real minimum sparking potential for a gap with this 
illumination is several percent lower than the normal 
sparking potential with weak light sources. 
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Fic. 5. Variation of the midgap streamer position with 


overvoltage. Percent overvoltage indicated at the right. 


would appear in the same breakdown. Increasing 
the intensity of illumination had the same effect 
on the position of the streamer as would be 
caused by increasing the overvoltage and vice 


versd. 
DISCUSSION OF RESULTS 


The short time lags observed can be explained 
on the basis of the theories proposed by Loeb" 
and by von Hippel and Franck" which were also 
developed by Schumann."': !* According to these 
ideas the breakdown is precipitated by high 
fields being built up by space charge formation in 
the gap. The space charge is generated by 
electrons crossing the gap and producing by 
impact electrons, positive ions and also radiating 
molecules. The their 
relatively large drift velocity in the electric field, 
are removed from the gap before the positive 
ions move an appreciable distance. Thus a net 
positive charge is left in the gap, which, under 
favorable conditions, may give rise to fields large 
enough for positive ions to produce a sufficient 
number of secondary electrons at the cathode to 
cause breakdown to occur. The radiation from 
molecules excited by electron impact may also be 
important in releasing electrons from the cathode. 
In any event the elapsed time is about that 
perhaps several 


electrons, because of 


required for an electron or 


electrons to cross the gap, and is of the order of 
magnitude of 10-* to 10-7 second for gaps of a 
few mm in air at atmospheric pressure. 


2. B. Loeb, Science 69, 509 (1929). 

‘A. von Hippel and J. Franck, Zeits. f. Physik 57, 696 
1929). 

“4 W.O. Schumann, Zeits. f. tech. Physik 11, 131 (1930) 
6 J. J. Simmer, Zeits. f. Physik 81, 490 (1933). 
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In interpreting the results obtained, account 
must be taken of the nature of the flash of 
illumination which was used to illuminate the 
cathode. The distribution of its intensity with 
time is not definitely known. Photographs of 
sparks taken by Snoddy* with a rotating mirror 
under somewhat similar conditions show that the 
intensity rises to a large value in a few times 10-7 
second and then dies out after several micro- 
seconds. Observations by Dunnington® and the 
writer’? show that an appreciable fraction of the 
maximum intensity is attained in a few times 
10-8 second. The photo-charge released from the 
cathode by a single flash of the illuminating 
spark at maximum intensity was measured and 
found to be about 10’ electrons. This is suffi- 
ciently large to make it probable that at least 
several photoelectrons were emitted during the 
first 10-* second, a time equal to about 1/500 of 
the total duration of the flash. 

The time lags measured lie in the period during 
which the illuminating spark intensity was 
rapidly increasing. Hence one must consider the 
rapid variation of photoemission with time in 
discussing the significance of the time lag- 
overvoltage curves. The rapidly increasing 
intensity of the illumination probably had a large 
effect on the shapes of the time lag curves. Thus, 
if the breakdown is due to the formation of space 
charge in the gap, and if the distribution of the 
space charge does not change appreciably, the 
condition for a spark to occur in a given gap is 
that a certain amount of charge, say gy, be 
generated. Then 


q e=No. of positive ions or electrons produced 


N 
= > x eSi ‘adzA Np, 


where a= Townsend coefficient for impact ionization by 
electrons, 
d=Gap length, 
No=No. of photoelectrons effective in producing the 
spark. 


The integral /(‘adx is a function of the amount of 
space charge already generated and hence must 
be evaluated for each photoelectron, the latter 
being designated as AN, in the above sum. 

The space charge generated in a given gap thus 
depends essentially on a and No. Decreasing the 
time lag also decreases N,) because of the smaller 


integrated light intensity of the initiatory spark. 
To compensate for this decrease, the value of a, 
and therefore the gap voltage, must be increased 
somewhat. At least a part of the gap voltage 
increment required to reduce the time lag by a 
given amount is needed to compensate for the 
correspondingly smaller number of photoelec- 
trons. Hence the rapidly increasing intensity of 
the illuminating spark doubtless accounts in 
a large measure for the rapid decrease in over- 
voltage at the longer intervals. The shortest time 
lags observed of 2 or 3X10~* second may be due 
to entirely inadequate illumination intensities 
produced by shortening the initiating spark. 
They may also be caused by the fact that at these 
time intervals the electrons are unable to cross 
the whole gap and therefore, in view of the 
shortened distance taken, the overvoltage must 
be increased to make the ionization more 
effective. The mobility of electrons in these high 
fields is uncertain by a factor of 3 or 4, but at 
least in order of magnitude the transit time of 
electrons across gaps of a few mm corresponds to 
the shortest time lags observed of 2 or 3X10-° 
second. 

The space charge generated by a single 
electron crossing a gap under the high fields 
corresponding to the shortest time lags observed 
is significant in this connection for if, under these 
gap conditions, the breakdown is initiated by the 
first few photoelectrons, the space charge field 
generated by them must be large enough to cause 
the spark to occur. As an example consider the 
3-mm gap at 50 percent overvoltage for which 
X/p=73. The value of a for X/p=73 is from 
Sanders’'® data 200, and neglecting space charge 
distortion, /o¢adx = ad = 200 X 0.3 = 60. From this, 
the number of positive ions generated by one 
photoelectron in crossing the gap would be 
e*’=e*, and the corresponding charge expressed 


in @.S.u., is 
g=4.8X10 oy e—=4.810-'™ x 10° =5 k 10", 


Such an enormous space charge could not be 
generated under the existing conditions and this 
must mean that the photoelectron traverses only 
a fraction of the gap before building up a large 
space charge. Thus, assuming that the spark is 
initiated by space charge effects, the breakdown 


°F. Sanders, Phys. Rev. 44, 1020 (1933). 
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under these conditions probably will start before 
the first photoelectron has crossed the gap. 
Increasing the gap voltage will reduce the time 


lag somewhat, but the limiting lag will be 
determined by the initiatory photoelectric 
emission. 

The space charge hypothesis provides a 


plausible explanation for the appearance and 
changes in position of the central streamer as 
observed in Fig. 5. The point or points at which 
the breakdown begins in a gap will be determined 
by the location of the space charge in the gap. 
For gaps with relatively weak fields it is possible 
that the space charge is set up through the 
migration of positive ions or the passage of 
successive electron avalanches in such a way as to 
concentrate the field at the cathode for short gaps 
and out in the gap itself for longer gaps. As the 
field is increased, fewer photoelectrons will be 
needed to produce the intense field and the field 
produced becomes more a field due to the 
ionization by a single electron avalanche. As this 
process takes place with increasing field strength, 
the enormously enhanced ionization along a 
single electron path as the electron approaches 
the anode will tend to concentrate the positive 
space charge at first near the anode, and, with 
increased field strength, back toward the midgap 
region. With an intense positive space charge in 
the midgap region the gradient which will be 
most effective will be located there and while 
some gradient exists at the cathode, the gradient 
in the midgap region is sufficient to produce a 
beginning spark there. When visible streamers 
appear, they appear simultaneously at the 
cathode and midgap for long gaps, and within 
2 or 3X10~° second. In short gaps the midgap 
streamer probably cannot be formed because the 
breakdown gradient will always occur first at the 
cathode." 

From the variation of overvoltage and time lag 
with illumination over the whole range investi- 
gated the following may be said. For time lags in 
the neighborhood of 10-7 second, the illuminating 
spark has had time to become rather bright, so 
that the number of photoelectrons released is 


7 The field in the neighborhood of the cathode surface is 
enhanced by the image of the positive space charge in the 
cathode. This will have a large effect in increasing the field 
there only if the space charge is relatively close to the 
cathode, as it must be in the shorter gaps. 


J 


WHITE 


large. Referring to our previous equation a very 
approximate condition for breakdown in a given 
gap that g Decreasing the 
illumination intensity means that the value of a, 


is be constant. 
and therefore the gap voltage, must be increased 
in order to keep the time lag constant. On the 
other hand, for very large overvoltages, the 
time lag is very small and only a few photo- 
electrons or even one are necessary to produce a 
breakdown. The illuminating spark intensity 
increases from zero at zero time lag, and an 
average time of several billionths of a second 
elapses before the several photoelectrons neces- 
sary to produce breakdown are emitted. Now 
decreasing the intensity will result in increasing 
this average time, and the time lag will be 
increased by just this amount. Increasing the 
overvoltage cannot compensate for decreasing 
the intensity in this region, for the condition for 
breakdown is that only a few, possibly one, 
photoelectrons appear in the gap. This is in 
contrast to the situation existing for the longer 
time lags where a decrease in the intensity can be 
compensated for by increasing the overvoltage. 
The curves for helium are similar to those for 
air except that a much higher overvoltage is 
necessary to produce a given time lag of break- 
down. The relatively large overvoltage necessary 
can be explained on the basis that ionization by 
electron impact occurs at low field strengths in 
helium, and at these field strengths is so ineffi- 
cient that it is necessary in relatively shorter gaps 
to increase a materially in order to produce the 
ionization necessary for a spark. Thus, for 
example, Werner'® finds that the ionizing zone in 
the positive corona discharge in helium comprises 
at least 500 free paths with a potential drop of 
about 200 volts, while in nitrogen about 460 
volts are required and only 70 free paths. This 
means that the rate of gain of energy is relatively 
small in the low fields in helium, while in the 
relatively large fields in nitrogen it is large. 
The writer is very much indebted to Professor 
L. B. Loeb for his interest in the writing of this 
paper and for his help in interpreting the 
experimental results. He also wishes to thank 
Dr. R. N. Varney for valuable criticism of the 


manuscript. 


‘SS. Werner, Zeits. f. Physik 90, 384 (1934). 
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On the Mean Lifetime of Metastable Neon Atoms* 


EUGENE W. PIKE, Scott Paper Co., Chester, Pennsylvania 


(Received February 1, 1936) 


Penning has shown that illumination of the Townsend 
discharge in neon with red light from an auxiliary neon arc 
reduces the number of excited atoms in the former. The 
concentration of excited atoms in the Townsend discharge 
as a function of illumination intensity can be computed 
from the electronic properties of neon; a comparison with 
Penning’s direct measurements evaluates the mean life- 
times of both the metastable and the resonance states in 


HE problem of the mean lifetimes of 

metastable atoms in neon, while less im- 
portant fundamentally than it once was thought 
to be, is still of great moment in the detailed 
mechanism of electrical conduction in gases. 
Up to the present all the experimental work has 
been done with the absorption method,' most of 
it before the true complexity of the problem was 
known. For this reason it seems valuable to 
point out a second independent method con- 
firming the previous results approximately, 
especially since the argument contributes ap- 
preciably to the mechanism of discharges. 

In a Townsend discharge in neon, excited 
atoms in both the metastable and the resonance 
states are generated in numbers comparable to 
the number of ions formed.? The atoms in the 
resonance state rapidly radiate their energy out 
of the gas, but the longer-lived metastable atoms 
accumulate to a high concentration.* If the 
discharge is illuminated with red light from an 
auxiliary discharge in neon, some of these 
metastable atoms are raised to the 18.5-volt 
levels (Paschen’s p levels) by absorption, and a 
certain fraction of these latter go down to the 
resonance states, instead of returning to the 
metastable states. The resulting decrease in 
metastable atom concentration can be measured 
directly,'* or by using a trace of some second 





* Presented at the Washington Meeting of the American 
Physical Society, April 27th, 1935. This work was done at 
Palmer Laboratory, Princeton, N. J. in 1934. 

‘ (a) K. Meissner and W. Graffunder, Ann. d. Physik 84, 
1038 (1927); (b) H. B. Dorgelo and T. P. K. Washington, 
Proc. Roy. Acad. Sci. Amst. 30, 37 (1927). 

*Cf. J. S. Townsend and S. P. McCallum, Proc. Roy. 
Soc. A124, 543 (1929); F. M. Penning and M. C. Teves, 
Physica 9, 97 (1929). 

*Cf. W. de Groot, Naturwiss. 14, 104 (1926). 


the gas. The numerical results for 36 mm pressure are: 
mean lifetime in sec. of the *P. metastable state = 2 x 10~*/ p; 
of the 2°P; resonance state=4X10~°/p, where p is the 
pressure in mm Hg. Values for the mean lifetime of the 
24P, state derived by Zemansky’s method from absorption 
measurements at about 3 mm Hg pressure are: 9X107* p 
from Meissner and Graffunder’s measurements, and 
8X 10-* p from Dorgelo and Washington. 


‘it can also 


gas such as argon as an indicator; 
be computed from the known electronic proper- 
ties of neon, using the mean lifetimes as ad- 
justable parameters. A comparison provides a 
numerical value for the mean lifetimes of both 
the metastable and the resonance states.° 

The equilibrium conditions for the metastable 
and resonance atom concentrations are simply 
formulated. Let the subscripts , and , refer to 
the resonance and metastable states, respec- 
tively; let E be the total number of excited 
atoms generated by the discharge per cc per sec., 
and m be the concentration of excited atoms;' 
let 8 be the reciprocal of the mean lifetime a, L 
be the illumination intensity in arbitrary units, 
2q be the probability that an excited atom will 
absorb a quantum from radiation of unit in- 
tensity in 1 sec., and } the probability that this 
atom will return to the opposite type of lower 
state, i.e., from a metastable and fo a resonance 
state, or vice versa. Writing the rate of generation 


on the left 
aE +2b.q,Ln,=n,(8,+2b,q,L), 
(1-—a)E+2b,9q,Ln,=n,(8,+20,.q.L). (1) 


‘F. M. Penning, (a) Zeits. f. Physik 46, 335 (1928); (b) 
ibid. 57, 723 (1929). 

°> By mean lifetime is meant here the concentration di- 
vided by the rate of generation; if energy of excitation is 
transferred from atom to atom in the gas by radiation or by 
collisions of the second kind, then the mean lifetime is the 
lifetime of the excitation energy, equal to the sum of the 
lifetimes of all the atoms excited by this energy as it 
traverses the gas. 

® Both £ and » are functions of position in the discharge. 
At pressures sufficiently high that diffusion is small, the 
relation between £ and nm will be true for any volume ele- 
ment in the discharge; since m is proportional to E (vide 
Eq. (2)), the relations developed below may be summed 
over the discharge as a whole. 
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It has been shown’ that slow electrons excite 
neon to the 18.5-volt levels only, rather than to 
the 16.5-volt levels. Because of this, a can be 
determined from Dorgelo and de Groot’s® pho- 
tometry of the light from neon excited by a 
22-volt electron beam (a sufficiently close ap- 
proximation to the Townsend discharge) by a 
simple summation of the line intensities to the 
various lower states, giving a=0.7. The excita- 
tion was sufficiently weak that the error arising 
from the re-absorption of the lines is second 
order. Because the *Py) state disappears from 
the gas so rapidly, the line intensity to it was 
divided equally between the metastable and 
resonance states. 

The “branching ratios’? can be determined, 
less accurately, from Dorgelo and de Groot’s* 
photometry of the light (‘‘Rezonanzlicht’’) re- 
emitted from a Townsend discharge irradiated 
by a strong beam of red light from an auxiliary 
neon arc. To within the error of measurement, 
b,=b,= 3. This ratio is probably dependent on 
the pressure and the illumination intensity, but 
experimentally the a’s are insensitive to varia- 
tions in 0 or (q, qs). 

The quantity qg is proportional to 2,B,,L»ps,° 
where the B's are the Einstein coefficients, and 
the L’s are the line intensities in emission from 
a positive column in neon; the sum is over all 
lines from the lower state in question, with the 
exception of the lines 5852, 6074, and 6402, 
whose upper states connect with only a single 
lower state.'® The Z values are given by Dorgelo 
and de Groot,* and the B’s are related to the A 
coefficients measured by Ladenburg, Kopfer- 
mann, and Levy" as B,,= (2 »/8s) (Ap. 82h)A ps. 
The sums for the two resonance states 'P; and 
3P, are 1655 and 1310, respectively; for the 
metastable *P) and *P» states, 790 and 680. 
Within the accuracy required here, g,=2g,. 
Implicit in this calculation is the assumption 
that the relation between the emission and 


*F. M. Penning, Zeits. f. Physik 41, 769 (1926); M. J. 
Druyvesteyn, ibid. 64, 781 (1930); 73, 33 (1931). 

* H. B. Dorgelo and W. de Groot, Zeits. f. Physik 36, 897 
(1926). 

® The subscripts , and , refer to the 18.5-volt and the 
16.5-volt levels, respectively. 

‘0 These lines were also omitted in the determination 
of b. 

“Cf, R. Ladenburg, Rev. Mod. Phys. 5, 248 (1933); 
Zeits. f. Physik 65, 189 (1930). 
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absorption contours is the same on the average 
for both types of states, probably a safe assump- 


tion here. 


Introducing these values into Eq. (1) and 
solving 
n,= E(q,.L+0.78,) [8.8,+9.L(8,+28 


Nm 


n, = E(2q,.+0.38,) /[8,8.+4sL(8-+28 


Since the comparison with experiment is made 
at constant field in the discharge, E£ is constant, 
and we may write (,+n,)=m"1, Mz,r<0) =m, 
(q.L,8,)=T, (8-/8,) =k, and solve for 


(mz,/mo) =(1+31/(0.7+0.3k) ]/[14+J(2+k)/k}. (3 


If a trace of a second gas such as argon is 
added to a self-maintaining Townsend discharge 
in neon, the potential )’7 across the discharge is 
greatly reduced, since the excited neon atoms 
ionize the argon atoms by collisions of the 
second kind.** If the conditions are otherwise 
constant and the argon concentration is very 
small, one may safely assume that V7, is a 
function of the product of the excited atom and 
argon concentrations. Penning’ has measured 
Vr in a neon-xenon mixture at 35.5 mm Hg 
pressure as a function of the _ illumination 
intensity and the xenon concentration vy; for a 
series of points at constant Vr, m_v_,=Mor, Or 
(nz/no)=(v/vz). Fig. 1 shows the curve of 
(n;,/mo) vs. L derived in this way. 

By proper adjustment of & and g,a, this curve 
may be fitted over its entire length by Eq. (3 
The appropriate values of the constants are 
k=50, g,a,=0.27. To show the sensitivity of the 
method, two other sets of points for nearby 
values of the constants are plotted. Clearly, 
qa, can be determined to a few percent, while 
k can be found to about twenty percent. Eq. (3 
may also be written in terms of g,a,; the best 
fit in this case is given by k=50, g,a,=0.01, or 
g.a,=0.25. If one sets g,=q,, instead of g,=24,, 
then the best fit results from k= 30, g,a,=0.28. 

Penning*® has estimated g, by an extension of 
the method used above to compute (4,/4); 
using the theoretical emission and absorption 
contours experimentally confirmed by Ladenburg 
and Levy,’ and Dorgelo’s measures of the efh- 
ciency of the neon positive column. His result, 


2) M. Penning, Physica 12, 66 (1932). 
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g,~0.5X10* in the units of reference 12, may 
be in error by a factor of two; combined with 
our value of g,a,, a,~6X10~5 sec. 

To complete the determinations of a, the 
writer has reduced the observations of Dorgelo 
and Washington!” by the method given in detail 
by Zemansky.'* Fig. 2a, plotted from their data, 
shows that the linear absorption vs. time relation 
predicted by Zemansky is actually fulfilled; 
Fig. 2b shows that K, the exponential decay 
constant for the metastable atoms in the absorp- 
tion cell, is described accurately by an equation 
of the form K=gD/p+ep, where p is the 
pressure. The significance is familiar—the gD 
term represents the loss of metastable atoms to 

8 A.C. G. Mitchell and M. W. Zemanskv, Resonance 


Radiation and Excited Atoms (Cambridge Univ. Press, 
1934), p. 249 ff. 
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Fic. 2. Reduction of absorption data from Dorgelo an 
Washington. 


the walls by diffusion, the ¢ term arises from 
destruction of the metastable atoms in the gas 
by collision with normal atoms. The numerical 
values are: a=1 gp=8X10-* p, diffusion con- 
stant D=220. The values derived by Zemansky 
from the observations of Meissner and Graf 
funder'* are: a=9X10~* p, D=470. (The dif- 
fusion constant of normal neon is 260.) 

Using these values of D, the loss to the walls 
in Penning’s discharge is negligible, so that a, is 
directly comparable with 1/¢, or a=210~*/p, 
in satisfactory agreement with the absorption 
values. Since practically all the errors in the 
absorption method increase the apparent decay 
constant, the most probable value of a is in the 
neighborhood of 10-*/ p sec. mm~ Hg. 

It is hoped to discuss the mechanism behind 
a in a later paper. 
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Penning’s ‘‘New Photo-Effect” in Pure Neon* 


EUGENE W. Pike, Scott Paper Company, Chester, Pennsylvania 


(Received February 1, 1936) 


The maintaining voltage of a Townsend discharge in 
pure neon was measured as a function of the illumination 
of the discharge with light from a neon arc, using a more 
sensitive circuit than that of Penning. In contrast to the 
Townsend discharge in slightly impure gas, in pure gas it is 
almost insensitive to illumination over a wide range of 
pressure, current, tube geometry and cathode material. 


With progressive cleaning of the gas, the change in 
maintaining voltage due to illumination passed from a 
small positive value through zero to a small negative value. 
The general equations of the discharge are developed, but a 
complete identification of the active atomic processes is not 
yet possible. 


: ' Abridged from a thesis presented to the Faculty of Princeton University in candicacy for the degree of Doctor of 
Philosophy. The work was done at Palmer Laboratory, Princeton, N. J., in 1934. 
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ROM experiments on the response of the 


Townsend discharge’ in neon containing 
traces of a second gas such as argon to illumina- 
tion with visible light from an intense auxiliary 
discharge in neon, Penning? has shown definitely 
that such illumination displaces the equilibrium 
between metastable atoms and resonance radia- 
tion in the discharge in favor of the latter, result- 
ing in a net decrease in the concentration of 
excited atoms in the discharge.* One would expect 
such a displacement to affect the maintaining 
voltage of the discharge in pure gas, but extrapo- 
lation of Penning’s curves to zero impurity (a 
matter of a few parts per million) indicates that 
this is not the case. The author undertook to 
measure the effects of this displacement in the 
purest possible neon, in the hope of setting upper 
efficiencies of various atomic 


limits to the 


pre cesses. 
I. EXPERIMENTAL 


A number of tubes, differing in unessential 
details, were used in the preliminary work. All 
had plane cathodes, all but one had plane parallel 
anodes. The odd tube had a grid of fine wire in a 
plane parallel to the cathode. A variety of 
cathode materials were used—iron, nickel, plat- 
inum, copper, Monel metal and platinum with a 
molecular layer of sodium. The ratio cathode 
diameter /electrode spacing varied from three to 
ten, and the pressure from 2 to 40 mm Hg. 

Fig. 1 shows tube //, the final form with which 
most of the data were taken. This tube was 
pumped on a fast, all-Pyrex system by standard 
high vacuum procedure; the residual pressure 
after sealing off was in the order of 10-* mm Hg. 
It was then reconnected to a gas-handling system 
through a fragile internal seal which was broken 
by a magnetic plunger after evacuation of the 
system, and pure neon introduced through a fine 
leak and a charcoal trap in liquid air. After again 
sealing off the neon was cleaned for about ten 
days by passing about | ampere between the 
misch metal in F as cathode and A as anode. It 


‘Used here to denote a space-charge-free, self-main- 
taining, cold-cathode discharge. 

°F. M. Penning, Phil. Mag. 11, 961 (1931); Zeits. f. 
Physik 57, 723 (1929). 
f. E. W. Pike, Phvs. Rev., this issue. 
C. C. Van Voorhis et al, Rev. Sci. 
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flecting shields; Q, seal-off; M, side tube to misch metal 
F, misch metal in iron cup; L, positive column discharge it 
neon; £, electrostatic shielding. 


was not possible to measure the final residual 
pressure of impurity, but aside from innocuous 
helium it was certainly less than 10-° mm Hg 
During observations the pressure could be ad- 
justed by withdrawing’ gas through a second in- 
ternal seal and about 30 cm of 1 mm capillary 
tubing immersed in liquid air. 

The electrical circuit finally adopted is shown 
in Fig. 2. This arrangement made it possible to 
follow major changes in ly, the potential across 
the discharge, and still measure to 0.1 volt the 
smaller changes (Al’r) illumination 
Observations of Al’r were made with the current 
in L alternately 1 amp and 0.01 amp, using (asa 
deflection instrument. 

The experimental results may be summarized 


due to 


as follows: 


(a) After normal cleaning, AV7 came to a constant value 
of one to two volts positive (i.e.,illumination increased Vr) 
This value was quite reproducible from tube to tube, and 
stable against further cleaning with misch metal. After one 
especially thorough pumping of tube H, however, removal 
of residual argon and mercury by immersion of a side tube 
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Fic. 2. Electrical circuit. B,, B., 600-volt batteries: 


P,, P»:, 6000-ohm potentiometers; 7, Townsend discharge 
tube; G, galvanometer, 10°-ohm sensitivity; Q, quadrant 
electrometer, sensitivity 4 mm volt; V, 0—750-voltmeter 
R,, 1.2X10°-ohm stabilizing resistor; Rs, 3X 10*-ohm 
damping resistor 


n liquid air reduced \ V7 to zero and eventually brought it 
ibout half a volt negative. The conditions were: oxidized 
tungsten cathode, 28 mm Hg pressure, 25°C, electrode 
spacing 9 mm, electrode diameter 30 mm, Vr: 525, 
This negative value of \l'7 was observed earlier in one 
other tube, but because it was falsely ascribed to leakage 
currents the conditions were not recorded. 

b) Aly is independent of the pressure over a wide 
range. Table | presents the best run on the variation of 
Aly with pressure, made with tube H 


TABLE I. Variation AV of the potential across the discharge 
produced by illumination vs. the pressure. 


PRESSURI 

mm He V7 Aly PRESSURI V7 Al'7 
20 430 +0.5 3.8 325 +1.0 
7.7 350 +1.0 3.0 328 +1.0 
5.0 330 +1.0 1.8 345 +0.9 


c) AV, is independent of the cathode material and the 
electrode geometry within the limits described above. 

d) Vr is independent of the current density between at 
least 3X10-° amp./cm*® and 3107-7 amp./cm*?. AV 7 is 
independent of the current density over a considerably 


wider range. 


All measurements were taken on the flat 
portion of the I’7-current curve. 


Il. THEORY 


The electron current passing through the neon 
creates by impact ions, and as well atoms excited 
to both the resonance and the metastable s states. 
Since all three of these are sufficiently energetic in 
neon to release secondary electrons from the 
cathode, the equation for the electron balance 
there must be modified to include them. Under 
the conditions of the experiment, all ions are 
drawn to the cathode. The energy represented 
by the atoms excited to the resonance state 
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diffuses by rapid alternation between radiation 
and excited atoms ; it crosses any plane in the gas 
primarily as radiation, and the fraction which is 
intercepted by the cathode releases photo- 
electrons. Some of the metastable atoms diffuse 
as such to the walls and the electrodes, but the 
major part are destroyed in the gas, their 
energy converted to radiation by impact of 
normal atoms in a process not yet fully under- 
stood; the cathode will intercept a fraction of 
both the diffusing metastables and the radiation. 

Taking these processes in order the total 
number of secondary electrons emitted, assuming 
infinite parallel plane electrodes and reducing to 
the basis of a single isolated electron leaving th« 
cathode,? is 


y(et*—1)+y/Af’ (ee47—1) +S 
+ty’[ f’(e*—1)-—S,—S,]=1, (1 


where the symbols have the following meanings: 


¢—Coefficient of secondary emission under ion impact 
corrected for back-diffusion.' 
d—Electrode separation. 
a—Townsend coefficient; the probability that an 
electron will form an ion by impact while moving a 
distance dx along the field is a dx. The bracket is 
the total number of ions formed. 
y’ Photoelectric efficiency of resonance radiation, 
also corrected for back diffusion 
f’'a—Cognate to a, for excitation to the resonance state 
\—The fraction of the energy expended in excitations 
to the resonance state which is intercepted by the 
cathode. 
¥’’—Coefficient of secondary electron emission from the 
cathode for metastable atom impact 
So, S:—The number of metastable atoms diffusing as such 
to the cathode and to the anode (and walls) 
respectively. 
f"a—Cognate to a, for excitation to the metastable 
state. 
ty’—Corresponds to Ay’, for the radiation resulting from 
destruction of metastable atoms in the gas. The 
ratio of the photoelectric efficiencies between this 
radiation and true resonance radiation is absorbed 


in the factor t. 


No account has been taken of ionization in the 
gas by ion impact, since under the conditions of 
the experiment the energy available from the 
field per ion free path is less than 0.1 volt. 





5 This reduction is possible because all processes in the 
discharge are linear in the current. (Cf. Paragraph d), 
Part I 

*°][. Langmuir, Phys. Rev. 38, 1656 (1931); C. Kenty, 
Phys. Rev. 44, 894 (1933). 
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The left-hand side of (1) must necessarily 
equal unity if the discharge is to be stable, since 
each electron leaving the cathode must on the 
average just replace itself. 

External illumination transfers some of the 
metastable atoms to the resonance state;> the 
energy intercepted by the cathode is similarly 
redistributed. In consequence, V7 changes to 
(V7+Alr), and a to (a+Aa). The equation for 
electron balance at the cathode now is 
Y(elatsad— 1) 4y/nf’ (e@tsed— 1) 4." Sy 


+[8/(8+8’) lw’ Lf’ (eets@4—1) —(So++S,*)] 


+p’ (8+ ,' Isv’C hf’ (e a+Ja tf ) 


Nm 


— (Sot +S,+) ]=1, \ 


where the new symbols are: 
Sy”, Si7 —Correspond to So, S;, corrected for the effect of 
the external illumination. 

8, 8’—These are the relative probabilities of destruc- 
tion of metastable atoms by normal atom impact 
and action of the external illumination, re- 
spectively; i.e., a metastable atom has a proba- 
bility Bdt of being transferred to the normal 
state by impact of a normal atom in time dt. 

s— The fraction absorbed by the cathode from the 
resonance radiation created by the external 
illumination from metastable atoms. 


Subtracting (1) from (2), and setting 
e*(1-+-Aa) = efetso, 
Civtw Af’ + tf’ et? JAad 
+ty’[(So+S$1) —(Sot+S1*) J 
—"(So— Sot) +L (s —0)8"/ (8 +8") ] 
x Lf" (efeta@4—1) — (Sot +$,+) ]=0. 
By comparison with (1), the first bracket is very 
nearly unity. With similar approximations, 
Aad =" (So— Sot) — tL (So+ Si) — (Sot + S17) J 


—[y/(s—2t)B", (B+ 8") ILf"%et4 — (Sot +517) ]. (3) 


III. Discussion 


The chief interest of Eq. (3) lies in the possi- 
bility which it offers of determining y’’, /, and 


\ PIKE 


(s—t),’ all at present unknown. The evaluation of 
these quantities involves accurate values of y 
and the S’s, however, which are at present un- 
obtainable. Kenty*’ has measured the photo- 
electric efficiency of an oxidized tungsten surfacy 
for neon resonance radiation, finding a value 0.03. 
but substitution of the data from Section I] (a 
into Eq. (1) gives an upper limit of 0.005 for y’, 
so that the unknown correction for back diffusion 
is evidently all-important. Similarly, the values 
of the S’s depend on 8 and on the diffusion 
constant of metastable neon atoms in normal 
neon,’ and the present uncertainty in these latte: 
allows a range of over 100 to one in the former 
The remainder of the constants are fairly accu- 
rately known.!° 

To obtain an upper limit of (s—/) from Eq. (3 
one must set a lower limit to ¥’; until further 
experimental values are obtainable this applica- 
tion of the theory must be deferred. 

One might note in passing that the experimen- 
tal work above documents an excellent method 
for detecting and measuring traces of impurities 
in neon and helium. 

In conclusion, it is a pleasure to thank Pro- 
fessor G. P. Harnwell for his generous super- 
vision of this work, and to acknowledge stimu- 
lating discussions with Professor R. Seeliger, of 
the University of Greifswald, and with Dr. C 


Kenty. 


7 A knowledge of (s—?) should throw light on the mech- 
anism which returns metastable neon atoms to the normal 
state by collision with normal atoms, since s is certainly 
based on the motion of resonance radiation in the gas 
and ¢ on the motion of whatever kind of radiation it is 
which carries off the energy of the metastable atom. The 
currently accepted theory, i.e., that the metastable atoms 
are raised to the resonance level by collisions of the first 
kind with normal atoms, is in conflict with some observa- 
tions of Found and Langmuir, and possibly also with the 
theory of molecular structure. It is hoped to discuss this 
fully in a later paper. 

*C. Kenty, Phys. Rev. 44, 891 (1934). 

°Cf. E. W. Pike, Phys. Rev. 40, 314 (1932). 

10 @, J. S. Townsend and S. P. McCallum, Phil. Mag. 8, 
657 (1928); (f’+f") J. S. Townsend and S. P. McCallum, 
Proc. Roy. Soc. Al24, 543 (1929); F. M. Penning and M 
C. Teves, Physica 9, 97 (1929); ff” cf. reference 3 
section ona. 6/8’ F. M. Penning, Physica 12, 66 (1932) 
and reference 3. 
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Capture of Slow Neutrons 


GG. Breit AND E. WIGNER, Institute for 


Advanced Study and Princeton University 


Received February 15, 1936 


Current theories of the large cross sections of slow 
neutrons are contradicted by frequent absence of strong 
scattering in good absorbers as well as the existence of 
resonance bands. These facts can be accounted for by 
supposing that in addition to the usual effect there exist 
transitions to virtual excitation states of the nucleus in 
which not only the captured neutron but, in addition to 
this, one of the particles of the original nucleus is in an 
excited state. Radiation damping due to the emission of 
y-rays broadens the resonance and reduces scattering in 
comparison with absorption by a large factor. Interaction 


1. INTRODUCTION 


ETHE,!' Fermi,? Perrin and Elsasser,*? Beck 

and Horsley* gave theories of the anoma- 
lously large cross sections of nuclei for the cap- 
ture of slow neutrons. These theories are es- 
sentially alike and explain the anomalously 
large capture cross sections as a sort of resonance 
of the s states of the incident particle. Resonance 
is usually helpful in causing a large scattering 
as well as a large probability of capture and it 
has been shown [H. B. Eq. (35) ] that large 
scattering is to be expected by nuclei showing 
anomalously large capture at thermal energies. 
This consequence of the current theories is ap- 
parently in contradiction with experiment, there 
being no evidence of a large scattering in good 
absorbers. It also follows from current theories 
that with very few exceptions the capture cross 
section should vary inversely as the velocity of 
the slow neutrons. Experiments on selective 
absorption recently performed® indicate that 
there are absorption bands characteristic of 
different nuclei and it appears from the experi- 
ments of Szilard® that these bands have fairly 


'H. A. Bethe, Phys. Rev. 47, 747 (1935). We refer to this 
paper as H. B. in the text. 

*E. Amaldi, O. d’Agostino, E. Fermi, B. Pontecorvo, 
F, Rasetti, E. Segré, Proc. Roy. Soc. A149, 522 (1935). 

* Perrin and Elsasser, Comptes rendus 200, 450 (1935). 

* Beck and Horsley, Phys. Rev. 47, 510 (1935). 

*Moon and Tillman, Nature 135, 904 (1935); Bjerge and 
Westcott, Proc. Roy. Soc. A150, 709 (1935); Arsimovitch, 
Kourtschatow, Miccovskii and Palibin, Comptes rendus 
200, 2159 (1935); Ridenour and Yost, Phys. Rev. 48, 383 
1935); Pontecorvo, Ricerca scientifica 6-7, 145 (1935). 

*L. Szilard, Nature 136, 950 (1935). 


with the nucleus is most probable through the s part of the 
incident wave. The higher the resonance region, the smallet 


will be the absorption. For a resonance region at 50 volts 


the cross section at resonance may be as high as 10°'® cm 
and 0.5X10°-°° cm? at thermal energy. The estimated 
probability of having a nuclear level in the low energ) 
region is sufficiently high to make the explanation reason 
able. Temperature effects and absorption of filtered radia 
tion point to the existence of bands which fit in with the 


present theory 


well-defined edges. It has been pointed out by 
Van Vleck’ that it is hard and probably impos 
sible to reconcile the difference in internal phase 
required by the Bethe-Fermi theory with 
reasonable pictures of the structure of the nu- 
cleus. The combined evidence of experimental 
results and theoretical expectation is thus against 
a literal acceptance of the current theories and 
it is our purpose to outline an extension which is 
capable of explaining the above facts by a 
mechanism similar to that used for the inverse 
of the Auger effect by Polanyi and Wigner.* 

It will be supposed that there exist quasi-sta- 
tionary (virtual) energy levels of the system 
nucleus+neutron which happen to fall in the 
region of thermal energies as well as somewhat 
above that region. The incident neutron will be 
supposed to pass from its incident state into 
the quasi-stationary level. The excited system 
formed by the nucleus and neutron will then 
jump into a lower level through the emission of 
y-radiation or perhaps at times in some other 
fashion. The presence of the quasi-stationary 
level, Q, will also affect scattering because the 
neutron can be returned to its free condition 
during the mean life of Q. If the probability of 
y-ray emission from Q were negligible there 
would be in fact strong scattering at the reso- 
nance, the scattering cross section being then 





7 J. H. Van Vieck, Phys. Rev. 48, 367 (1935). 

®O. K. Rice, Phys. Rev. 33, 748 (1929); 35, 1551 (1930); 
38, 1943 (1931); J. Chem. Phys. 1, 375 (1933). A similar 
process was used by M. Polanyi and E. Wigner, Zeits. { 
Physik 33, 429 (1925). 
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of the order of the square of the wave-length. 
Estimates of order of magnitude show that it is 
reasonable to assign 12 volts to the “‘half-value 
breadth” of Q due to radiation damping and that 
the ‘“‘half-value breadth’ due to passing back 
into the free state is about one-fortieth of the 
above amount. This means that when the system 
passes into the state Q it radiates practically 
immediately and the neutron has no time to be 
rescattered. It will, in fact, be seen from the 
calculations that follow that the ratio of scatter- 
ing to absorption is essentially the ratio of the 
corresponding half value breadths. The hardness 
of the emitted y-rays is of primary importance 
for the small ratio of scattering to absorption 
the y-Tay 

high. the 


makes probability of 


because it 


sufficiently Inasmuch as 


emission 
interesting phenomena occur for low energies 
we may suppose that in most cases the coupling 
of the incident state occurs through its s state, 
i.e., in virtue of head on collisions. It will be 
seen, however, that the possibility of obtaining 
observable effects by means of p states is not 
excluded even though it is less probable and leads 
to smaller cross sections. Calculation shows that 
with resonances of the type considered here one 
may obtain appreciable probability of capture at 
energies of the order of 1000 volts. It is possible 
to have at such energies cross sections of roughly 
10-* cm? with a half-value breadth of about 20 
volts. It is therefore not necessary to ascribe 
all large cross sections to neutrons of thermal 
velocities and the probability of finding a 
quasi-stationary level in a suitable region is not 
so small as to make the process improbable. 

We are presenting below the theory of capture 
on this basis in some detail not because we be- 
lieve it to be a final theory but because further 
development may be helped by having the pre- 


paratory structure well cemented. 


2. THEORY OF DAMPING 


The process of absorption from the continuum 
into a quasi-stationary level and a subsequent 
reemission of a photon is related to the phe- 
nomena of predissociation discussed by O. K. 
Rice* who made the first application of quantum 
mechanics to this type of process since Dirac’s 


AND 


E WIGNER 


first approach.® It is essential for us to consider 


two continua and in this respect the present 
problem is more general. It resembles closely the 
problem of absorption of light from a level a 
to a level c which is strongly damped by radia- 
tion in jumps to a third level 6. The absorption 
from a to ¢ corresponds to the transition of the 
neutron into the quasi-stationary level and the 
jumps from c to } correspond to the emission of 
y-rays in a transition to a more stable level of the 
nucleus. The absorption probabilities can be ob 
tained by using the principle of detailed balance 
from the solution which represents emission! 
from the level c to the levels a, } or else by a 
direct application of the theory of absorption! 
The usual theory as developed for either process 
is not accurate enough to represent the effect 
of the variation of matrix elements with velocity 
which is essential for our purpose, inasmuch 
as it is responsible for the existence of two regions 
of large absorption. The usual type of calcula- 
tion will now be generalized so as to take the 


Variation into account. 


(a) Calculation of the absorption and scattering 
process 


Let a, denote the probability amplitudes of 
states in which the neutron is free and in a state 
s. Similarly let 6, stand for the probability 
amplitude of a state in which the neutron is 
captured and there is a photon r emitted and let 
c be the probability amplitude of the quasi- 
stationary state having energy /iv. The states r, 
s are here considered to be discrete but very 
closely spaced in energy. The average spacing of 
the levels 7, s are written AE,=hAv,, \E,=hdv 
so that the number of levels s per unit energy 
range is 1 AE,. The matrix element of the inter- 
action energy responsible for transitions from 
a, toc, c to b, will be written, respectively, 


M,=hA,, M,=hB,. I 


The damping constants for c due, respectively, 
to the possibility of emitting a, or 6, are then” 


*P. A. M. Dirac, Zeits. f. Physik 44, 594 (1927). , 
WV, Weisskopf and E. Wigner, Zeits. f. Physik 63, 54 
1930). 

'"'\’, Weisskopf, Ann. d. Physik 9, 23 
2G. Breit, Rev. Mod. Phys. 5, 91, 104 


1931) 
117 (1933) 
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CAPTURE OF SI 


trr,) '=l,=[m! A,|?/Av, |,,2003 2) 
dir, l=l,=[2 B ?/Av, |, », (=I,4+FT,, 
where 7,, T» are respective mean lives of c due 


to emission of a, and b,. The quantities I’ repre- 
sent one-half of the “half-value breadth” 
measured in frequency. In discussing line emis- 
sion and absorption '°: ' the directional averages 
of A, * and 8B, * can be taken for any energy 
within the breadth of the line because the line 
can be usually considered to be sharp. In the 
present case it will be necessary to distinguish 
among directional averages of |A,'* for differ- 
ent energies. 

The states s will be thought of as plane waves 
modified by a central field due to the nucleus and 
satisfying boundary conditions at the surface 
of a fundamental cube of volume I’. The equa- 


tions satished by a,, },, c are 


d d 
( +p, Ja, =A,c; ( +yv, }b,=B,c, 
2ridt 2ridt (3) 


and 
as,=e °* - c=zA,,*e ** 
a,.=A.A.,*{ e~** ('_@ 2" 


In Eq. (6) s#5». The quantities y and »—», 
are small compared with I’; they will go to zero 
with increasing volume. From (3), one finds for 
them the equation: 

(vy, — vo tty) (v— vy—tl) =!A,, |? (7) 
so that 
Y= A,, *T'/[(v— 9)? +1? ); 


Vo=Vs,t+(ve—v)(y/T). (8) 
In obtaining Eq (7) the approximations 


: 1 —e?ti (vo ve 
>.’ |A,|? 


— 


=mt A,\*/Av, (9) 
Vs— Vo) 


are made. These correspond to replacing the 
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d 
( +y)e= EA." +> B,*b,. 
2ridt 


In these equations the influence of only one 


quasi-stationary level is taken into account and 
for this reason they are not quite accurate. They 
are sufficiently good for the present purpose 
because it will be supposed that different quasi- 
stationary levels do not fall closely together. At 
t=0 it will be supposed that 


Be = Once, b az (), c=), i=). 4 


An approximate solution of (3) satisfying this 
initial condition can be obtained by forming a 
linear combination of 


c=e- : 
a,=A,fle ies ; —@e ** (y,—v+il"’), 5) 
b,=B,Le~** Pt e-2eerrt) /(y, —y+al”), 
with 
| resor ' regior (5°) 
(y—vyo—tl 
“T/ (¥s— vo tty) (v— v9 — 17) 
(6) 
'\/(v-— votty) (y— vo —1P) 


sums by integrals and extending the range of 
integration from »v,.=—* to »v,=+2% and 
similarly for v,. In addition it is supposed that 
A,|?, |B, 
in which the integrand is large that they may be 


vary so slowly through the region 


taken outside the integral sign. These approxi- 
mations are, therefore, valid only if the contribu- 
tions to the sums (9a), (9b) are localized in a 
sharp maximum. Such a maximum exists for 
v, =v because: (1) y vanishes as the fundamental 
volume is increased and therefore one may con- 
sider yt1 and (2) for any v,—v» it is possible to 
—vo (>1. 
For such times the most important part of the 


choose ¢ sufficiently large to make |», 


integrand oscillates rapidly with v,. However for 


v.—Vvo, ~y, the values of ¢ which satisfy y/<1 
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t<1. The integrand 
the 


are always such that v,—y, 
is thus not oscillatory for v,.=vpty and 
values of A, ?, B,? on the right side of (9 
are to be understood as corresponding to v,.=y 
with an uncertainty of the order y. It can be 
verified by calculation that the contribution to 
(9) due toa finite region at a distance |v,—v9| >>y7 
contributes imaginary quantities decreasing ex- 
ponentially with 27 v,—v»9 ¢ and real quantities 
which contribute to a frequency shift® of ». For 
the present this shift will be neglected. Eqs. (6) 
are thus approximate solutions which become 
increasingly better as ¢ increases, provided 
y!<1. In our application [ is mostly due to the 
radiation damping I,. The directional averages 
of B,* vary smoothly since the energy of the 
y-ray is of the order of several million volts and 
is large compared to I. 

The quantity I’ which enters (5) is not de- 
accurately by the method 
A, * which enters in this case is some 


termined present 
because 
sort of average over the resonance width. This 
complication causes no trouble because: (a) 
for times £>1,/4rI’ 
states a,, >, are, respectively, 47yI,/T, 477°, T 
and depend only on T and not on I’; (b) the 
in comparison with I, makes 


the rates of emission of 


largeness of T 
r’—l'<«r. Thus I’ 
determining the initial transients but not the 
steady rate of absorption. This can be expected 
from the fact that the solutions (6) represent a 
the rate 
solution” 


is of importance only in 


condition in which sp» is absorbed at 
4ry. The addition of the 


(5) is only needed to enforce the condition c=0 


“emission 


at /=0; it modifies the emission of states 6,, a, 
during times comparable with the mean life of 
the nucleus but leaves them unchanged over 
longer periods very similarly to the way in which 
analogous transient conditions are of no impor- 
tance in the absorption of monochromatic 
radiation by classical vibrating systems. 

The total cross section ¢ which corresponds to 
the disappearance of the incident states Sp is 
given by 


(10) 


c= diary V/2, 
where v is the neutron velocity because the 
modified plane waves denoted by s were normal- 


8 Appendix | 


AND E 


WIGNER 


ized in the volume |’ and thus represent states 
of density 1. 
The number of possible plane waves in |" per 


unit frequency range is 


1/Av,=42 V /vA?, 11 
where A is the de Broglie wave-length. From (2), 
(10), (11) we have 
A* i 
o =yA*/Av,=—S 12 
nm (vy—vo)?+TI° 


Here the statistical factor S takes account of the 
fact that the state s) may be more or less effec- 
tive in its coupling to the quasi-stationary level 
than the average modified plane wave in the 
same energy region. If the quasi-stationary level 
has an orbital angular momentum Lf and if 
there is no spin orbit interaction then A 

=(2L+1)|A,!* because coupling to c can take 


place only through 1/(22+1) of the total 
number of states. Thus. 
S=2L+1 13 


in these special circumstances. For s terms S=1 
The total cross section 


o=0.4+0:%, 


where o, is the cross section due to scattering 
and o, is the cross section due to capture. We 


have 
A? eg A? i," 

o¢.=—S - ¢g,=—S 14 
r (vy—vo)?+I? xr (v—vo)?*+I° 


The above value of ¢, corresponds to the value 
+ a,\? and does not take into account the fact 
that there is scattering in the abscence of the 
quasi-stationary level. If this is strong one must 
correct o, for interference of the states s with 
the spherical wave present in so. In the applica- 
tions made below the scattering effect due to 
either cause will be small and the correction need 
not be considered. According to (14) the extra 
scattering can be expected to be of the order 
Def hs 
small T’,. 

It should be noted that the order of magnitude 
of o, at resonance is changed by taking into 
account the radiation damping. If this were 


times the capture and is quite small for 
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neglected and if one were to calculate simply by 


using Einstein’s emission probability for the 


stationary states of matter then one would ob- 
tain an incorrect value, 


o‘a- S ; ~* (14°) 


For resonance o,¢,'=I'7/T? and approximately 
SodE fo.'dE is T,/T,. No paradox is involved 
here because it is not legitimate to apply Ein- 
stein’s emission probability formula to levels 
separated by less than their breadth due to 
radiation damping. Eq. (14’) gives too high 
values to the cross section. If »—v» >IT there is 
no difference between ¢.’ and ¢,. For sufficiently 
large values of y—vo the discussion which led to 
Eqs. (13), (14) will break down because Dirac's 
frequency shift® is neglected in these formulas. 
A more complete formal discussion including 
the frequency shift is given in Appendix I. The 
calculation shows that one should change the 
frequency of the quasi-stationary level v by 


A,|? dv, »|B,|? dp, 
por-— -| - (15) 
Avs Vs—VY Avy v¥r—Vo 


where the integrations are extended over the 
complete range of states s, r and where the 
principal values of the integrals are to be taken. 
The last part of (15) represents the frequency 
shift due to electromagnetic radiation and can 
be incorporated in v as a constant because vp need 
be varied only in a range small in comparison 
with the frequency of the y-ray. It is dangerous 
to take this shift into account on account of the 
well known inconsistency of quantum electro- 
dynamics. The second term on the right side of 
Eq. (15) is due to interactions between free 
neutron states and the quasi-stationary state. 
It is physically correct and it is necessary in 
order to bring about agreement between (14) 
and calculations away from resonance by means 
of the Einstein emission probabilities. The shift 
is large in the applications. Nevertheless changes 
in it are small in the relatively small range of 
values which need be considered and its effect 
is therefore primarily that of displacing the 
resonance frequency by a constant amount. 


4t 
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b) Resonance of one-body systems 


The above discussion cannot be applied directly to 
cases in which resonance consists simply in a sharp increase 
of the wave function of one neutron to a maximum inside 
the nucleus because there is no intermediate state c under 
such conditions in the same sense as in the previous section 
For low velocity neutrons such resonance can be sharp 
for states with L=1. Formally one could try to appl) 
the discussion already given by starting with wave func 
tions which are solutions of the wave equation for an 
infinitely high barrier somewhat outside the nucleus. The 
difference between the actual height of the barrier and = 
can be then treated asa perturbation essentially responsible 
for the matrix elements AA,. Such a procedure leads ap 
parently to correct results which can be verified by other 
methods. It is troublesome to justify it completely because 
the region where the infinite barrier must be erected should 
be such that the wave functions within are small for all 
energies. It is preferable to use a more direct calculation for 
such a case. We consider a plane wave of neutrons incident 
on the nucleus. Resonance takes plac e to the wave functions 
of angular momentum LA. We surround the nucleus by a 
large perfectly reflecting sphere of radius R and we calcu- 
late the rate at which states of angular momentum Lh 
disappear by radiation. There is no essential restriction on 
the possibility of forming wave packets out of the plane 
waves if we admit only those states L which satisfy the 
boundary conditions on the sphere. The radius will be 
made finally infinitely large and the spacing between the 
levels infinitely small. This provides the necessary flexi- 
bility for the formation of the wave packets. 

The spacing between successive possible neutron levels is 
given by 
Av=0/2R. 16) 


The radial function will be expressed as F/r where F will be 
by definition a sine wave with unit amplitude at a large 
distance from the nucleus. The normalized wave function 
is then Y,(F/r)(2/R)* where Y_, is a spherical harmonic 
normalized so asto have f Y,|*d2=1. The wave function 


for the bound state will be written 
Vioif/r; | ftdr=1. 17) 


The damping constant which corresponds to the emission of 
radiation from the state F is obtained by using the formula 
for Einstein's emission probability and is 
ve=(C/R)| f- Ffrdr|?, 18) 
~0 
where 
7 32x%e’ Y L+ } = } . 
C= : 18’) 
she 2L+1 


the upper sign applying to jumps L~L+1 and e’~e/2 is 
the effective charge of the neutron nucleus system. As 
R—-x, both Av and yg decrease towards zero but their 
ratio remains constant. The cross section due to capture 


computed directly from the emission probability is 


oc’ =(2L4+1)A*ye/Av 19) 


If this expression approaches A? then yg¢/Av becomes com- 
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parable with unity and the levels are close enough together 
to make Eq. (19) meaningless. It is then necessary to take 
into account the mutual influence of neighboring levels. 
This can be done by means of the damping matrix.'* The 
successive states of angular momentum Lh will be denoted 
by indices j, / and their probability amplitudes by a;. 


d ’ 
( ~+ pv; Ja; =1Dy! ‘ay, 20) 
2ridt 


(20’) 


These satisfy 


where 


yt=CJ;*Ji/R; Jj=S Fifrdr. 


In our case only states with the same magnetic quantum 
number can interact so that a complete specification of the 
states is obtained through their energy. Solutions of (20) in 
which all quantities vary as exp | —27t(vo—ty)t} correspond 
as closely as possible to the notion of a stationary state 
decaying under influence of radiation damping. From (20) 
one obtains 

a votty)a; =4D7’ al. 20’’) 
These equations with the complex eigenwert »»o—?y can be 
reduced making use of the fact that y’‘ is a matrix of rank 1 


Thus eliminating the a; one finds 
sC J;\*? 
lo—3 
R vj—votty 


(21) 


for the secular equation which determines vo and y. This 
equation will be solved approximately for the case of 
sharp resonance. The resonance will be supposed to take 
place at an energy hve and to have a “half-value breadth” 
2hT'r. Close to resonance 
I'p?| Z|? 
J,|\*o— (22) 


(vj— vp)?+ lp? ° 


where |J|* is the maximum value of | J|*?. This approxima- 
tion will usually apply only in a region of a few I'v. The 


value of I'p can be estimated using '° 


4aT'p =0,/ f G'dr, (23) 


where G is F for resonance, 2; is the velocity at resonance, 
and the integration is to be carried through the range of 
large values of G. The quantity I'y is analogous to I’, of 
section (a). The state represented by G is analogous to the 
quasi-stationary state of section (a). In order to bring out 
the analogy we introduce a damping constant similar to 
the previous I, 


Tr=C|I|2/2SGedr =2xC|1\*0p/2,, (23’) 


which is the damping constant of the state represented by G 
when that state is normalized within the nucleus and its 
immediate vicinity. Substituting (23’) into (21), replacing 
the sum by an integral everywhere except in the vicinity of 
vo and performing the summation in that region on the 
assumption that the Av can be considered as equal to each 
other in that region gives 


4G. Breit, Rev. Mod. Phys. 5, 117 (1933); G. Breit and 
I. S. Lowen, Phys. Rev. 46, 590 (1934). 

16 G, Breit and F. L. Yost, Phys. Rev. 48, 203 (1935). See 
also Eq. (32’). 
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= w(vjo—votty) 
1=t\acot]- +1ta 


Av;o 
‘i Ur Jj *T'pdv; 
’ 
24) 
“0 xt; I\*T p(y vot 
v, J; "TR 
é= 7?r ’ 
v;|4/*T'r 4,, 
where the integral must be extended over all »; and the 


region around yo is integrated over the real axis. The 
quantity vjo is any one of the v; located so close to vy that 
the variation in Av in between can be neglected. In the 
approximation of Eq. (22) the integration over the reson- 
ance region I'y leads to an equation which to within a 


sufficient approximation reduces to 


— = vjo| J |2¢ 
1+Th=(Th+it)(a+12b); b=- : 25 
vr| J jo|72(1+¢ 
with 
; a Ty w(vo— 
vo—vp=ql'r; Th=tanh t=tan 25 
vj0 Vio 
By eliminating ¢ 
Th+1 Th=a +1/a+ b a, 25 
which has the approximate solution 
1/Th=a+1/a+0"/a. 26 
For values of vo which lie in the region where Eq. (22 
applies and where Avjo~Av, we have approximate 
my rr 
26 


Avy V'R?+Tp*(1+q?)’ 
where it is supposed that 'g>>Tp. If, however, I'r >T'z then 
wy Iz 
=  . 26 
Av 'p(1+g?) 
which is equivalent to using the yg of Eq. (18), (19); in this 
case one may compute using emission probabilities. If one 
is so far away from resonance that b?/a<a, 1/a Eq. (25’ 
gives 
tit 26 
provided the right side is <1. Here again the simple 
emission point of view applies. For I'y>>I'p all regions are 
approximated by 
ry v-|Jjo\*U Rl 1+ ¢°) oe 
Avjo vjo|7\*( 0x? +P r*(1+ *) ] 
The treatment of scattering by means of the damping 
matrix is somewhat involved and will not be reproduced 
here. The phase shift due to »o—vjo when added to the 


phase shift already present in Fjo gives the phase shift 
required. The scattering is diminished by Ig in much the 


same way as it was diminished by it in section (a). By 
comparing (27) with (19) 

A*v, | J jo |? eT r( 1+") . 

oc = (2L+1)- aM od len 28) 


WV jo I *( TR? T Ip" 1+") 


which is similar to Eq. (14), close to resonance. The factors 
| J|?/|Z|? and v,/v;9 take into account the deviations from 
the dependence of | /J{|* 14) this is 


on » given by (22). In 
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analogous to the dependence of I’, on v.9 combined with 
Dirac’s frequency shift. 


c) Sharpness of resonance for single-body problem 


The upper limit of integration in Eq. (23) has been left 
indefinite. By Green’s theorem 


; [ # Sa Py |+ “(E2—E:) Fi F:=0, 
dr dr dr h* 


where F;, F: correspond to energies /,, KE, and need not be 


regular at r=0. Hence * 


a Oo OF 2u 
[= ; |+ F2=0, 29) 
OrL OF Far h? 


In this section let F; be the function inside the nucleus, and 
let F stand for the regular solution of the wave equation for 
rXradial function on the absence of the nuclear field. The 
normalization of F is such as to make it a sine wave 
sin (kr+¢) of unit amplitude at «. Similarly G is defined 
as satisfying the same differential equation as F but it is to 
be 90° out of phase with F at = i.e. cos (kr+¢). The 
regular solution of the differential equation in the presence 
of the nuclear field, normalized in the same way as Fand G, 
will be called F, At the nuclear radius ro 


pao /{fo(E-F)}+Lee6-2)]}- om 


Here the accent stands for differentiation with respect to 
kr. At resonance F,'/F;=G'/G and the second term in the 
curly bracket is then 1, while the first term is zero. As E 
changes to either side of the resonance value £, the first 
term may become 1 for E=F,+AF where SE is properly 
chosen. The half-value breadth is then 2AE and AE =AT,y. 


The value of AE can be estimated by 


Of (G F, 
AE- | - ~) =I. (31) 
dE G F; ro 


Using Eq. (29) and calculation the 0/dE for E=E, one 
obtains a result which can be expressed in terms of integrals 
up to R where R is any value of r which is greater than 1. 
The function which is G for r>ro and F;(G/F;),, for r<ro 
is continuous at ro and at resonance its derivative with 
respect to r is also continuous. The function will be called G 
for 0<r< x. We have then " 


E 


"R_. G°Ed aG 
seme, Gdr +[ : 


; AE=h"y. 32 
= kaE a al \ ) 


rhe right side of this result is independent of R and is 
hnite. The term outisde the integral should be included in 
Eq. (23) changing 





Sj. A. Wheeler. We are indebted to Dr. Wheeler for 
communicating to us other applications of this relation. 

."' Cf. Eq. (22) reference 15. In calculations with Coulom- 
bian fields it is sometimes convenient to transform Eq. (32) 
of the text into 


E k (7 k (‘% EG*a ( kr 
ioe) &. 2dr — lin See 
AE [Fs J, ° Fear rJ ~ ( 7) | 


all quantities outside the integrals being taken for r=rpo. 
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Eq. (32) has a well-defined meaning only if resonance is 
sharp. Otherwise the 0/dE entering in Eq. (31) cannot be 
supposed to be sufficiently constant through the half- 
breadth 2hT'p. It cannot be expected to hold for the broad S 
resonance discussed by Bethe. 


(d) Capture by / states 


For a potential well of constant depth 


F;=sin s/s—cos 2; F=sin p/ p—Cos p; 
G=cos p/p+sin p, (33) 

where 

z=Kr; p=kr; K=wpv,/h; k=pr/h (33’) 


v;, v being, respectively, the velocities inside and 
outside the nucleus. The resonance condition is 


zsin z [sin z/2—cos z|= p COS p [cos p/p+sin p |. 


For slow neutrons p<1 the right side is <p* and 
therefore very small. The first resonance point 
is obtained for z=r—e, e~p?/z. It will suffice 
to take s=7. By substituting into Eq. (32) it 
follows that 


AE/E=2p/3=4ar,/3A. (sorry 


For E=(1/40) volt, A=1.8X10-* cm, Ar =5.8 
x10°* volt. For 3-MEV y-rays a reasonable 
value of AI'z is 5.8 volts. The cross section at 
resonance is by Eq. (28) 3A°Tr/rlr=300 
xX10-* cm*. Since scattering is of the order 
I'y/T pr times capture the scattering cross section 
is small. According to Eq. (33’’) the cross 
section at resonance for p terms with ['p>I'r 
can be expected to vary as v and Alp as vo’. 
The range in which p terms can be expected to 
give large capture cross sections and small 
scattering is therefore roughly from 1/40 volt to 
1 volt. At higher velocities Aly is likely to be 
higher than AI’,. In the absence of an apparent 
reason for nuclear p levels to fall in this narrow 
velocity range, an explanation in terms of p 
terms although possible is improbable on account 
of the small range of neutron velocities required. 


3. CAPTURE THROUGH s WAVE 


(a) 


This section will contain the calculation of the 
A, used in 2a. It is supposed that the system 
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‘“‘nucleus+neutron” can be treated in first ap- 
proximation by means of an effective central field 
acting on the neutron. The difference between the 
Hamiltonian of the system and the Hamiltonian 
corresponding to the central field will be called 


H’. On 


transitions from the s wave of the incident state 


account of this difference there exist 
to quasi-stationary excited states of the ‘‘nucleus 
+neutron” system. Normalizing the s waves 
within a sphere of radius R the wave function 


inside the nucleus is 


Csin Kr/r; C-*=[1+(U,/E) cos? Kr j27R, 


where K?/k?=(U+E)/E 
and U is the depth of the potential hole. The 
energy [’ r also 
internal coordinates x. The wave function of the 
system in the incident may be 
Cyo(x) sin Kr/r 
stationary state Wo(r, x). 
M, of Eq. (1) is then 


interaction involves besides 
state 

the 
matrix element 


whole 
and in quasi- 


The 


written 


M,= J volr, x)H’C sin Krdo(x)dv/r, (34) 


where dv is the volume element of the whole 
system. The state Q is by definition such that 
the integral of We.? through nuclear dimensions 
is unity. The order of magnitude of .\W/, is 


therefore 
M,=CHr,', (34’) 


where // is an average of H’ through the nucleus 
and may have reasonably a value of 0.5 MEV. 
It cannot be specified further without detailed 
calculation which would probably be unsatis- 
factory in the present state of nuclear theory. 
Since AE =hv,2R, 
IT’r, 
AT. = : (34’’) 
2AU cos? Kro 


According to Eq. (14) 
Aro H Hhr, 
2x Ucos? Kro h?(l,+T.)?+(E—E,)? 


is the value of E for resonance. Ac- 


where E, 
cording to this formula there are two maxima 


AND E \ 
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PABLE I. Calculated cross sections for neutron cap 
POSITION OF Avs 10249 
RESONANCE H hl, AT RESONANCE 

volts ME\ volts volts RESONANCE E> 

1/40 0.1 10 0.01 90 000 

a 1 01 900 000 
1 1 10 £05 14 000 
l 1 05 140 000 
50 1 10 37 2 000 3500 
5 10 9 13 400 80000 
1 1 0.37 11 000 350) 
a l 9 1 800 9000 
1000 B 10 1.6 320 9 
5 10 40 $20 200 
10000 1 10 5 60 0.09 
5 10 125 18 2 
for o., one for E=E, and one for E=0. The 


expected cross sections are given in Table I to 
about ten percent accuracy. The numbers corre- 
spond to A(kT)=1.8X10-§ cm; A(1 volt 
x10-° cm; ro=3X10°™" U cos? Kr, =10' 
volts. For £,=1/40, 1 volt the table shows large 
cross sections at thermal energies and above 


=2.9 


cm; 


The condition is similar to Bethe’s except for a 
relatively sharper resonance determined by AT 
For 50 volts one sees the development of two 
maxima one at resonance and one at thermal 
EF, = 1000 10,000 


maximum at thermal energies decreases as E 


energies. For and volts the 
and the maximum at resonance roughly as E 


For such high values of £, scattering has a 
chance of becoming comparable with absorption 
or even greater than the absorption at resonance. 
In the thermal energy region the 1 7 law is 
low E, the 


v law and 


obeyed for high values of £,; for 
maximum at £, interferes with the 1 
the region of its validity is displaced below 
thermal energies. 

In Table I only the effect of a quasi-stationary 
level at E 
be effects of other levels as well as radiation 
Bethe and 


Fermi which do not depend on the existence of 


is considered. In addition there may 
jumps of the kind considered by 


virtual levels. It is thus probable that in most 
cases there is a region with a 1/v dependence 
masked by a 


although it may be at times 


resonance region. 

(b) Dirac’s frequency shift 
In the above estimates the effect of Dirac’s 

frequency shift was neglected. This is given by 
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Pht, dE 
« T E-E, 


0 


HI°ro x°*dx 


ore (x? — x9?) (x? +a?) 


where x=F?, a?=Ucos* Kry and the value of 
ht, was substituted by means of Eq. (34”’). 
Here the subscript 0 refers to the neutron energy 
E, and the principal value of the integral is 
understood. Evaluating the expression 


IT’r,U cos Kro 
(hAv) p= ~ a, (36) 
2Ac( Ey + U cos? Kro) Eo! 


The shift is seen to be of the order of 3000 times 
hY, for Eg=1 volt. The shift is nearly inde- 
pendent of the velocity. In the approximation of 
Eq. (36) 


d(hAv)p Ahr, Ep} 
== - (36') 


dE, E, U?\cos Kro 


which shows that the variation in the shift ts 
small and of the order of 2K10°°(4—£E,) for 
H=0.1 MEV. 


4. Disct SSION 


(a) Absence of scattering 


According to Dunning, Pegram, Fink and 
Mitchell'® the elastic scattering of slow neutrons 
by Cd is less than one percent of the number 
captured. According to A. C. G. Mitchell and 
E. J. Murphy" scattering as detected by silver 


s about the same as absorption for Fe, Pb, Cu, 
Zn, Sn while for Hg scattering is about 1/80 of 
the absorption. In the later communication of 
Mitchell and Murphy" it is also found that Ag, 
Hg, Cd are poor scatterers of slow neutrons 
detected by silver. It is interesting that Ag 
shows small scattering in these experiments 
because the detection took place by means of 
silver and that Hg and Cd show small scattering 
because they have large absorption cross sec- 





'8 J. R. Dunning, G. B. Pegram, G. A. Fink and D. P. 
Mitchell, Phys. Rev. 48, 265 (1935). 

had 5 C. G. Mitchell and E. J. Murphy, Phys. Rev. 48, 
653 (1935). Cf. also Bull. Am. Phys. Soc. 11, paper 27, 
Feb. 4, 1936. 
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tions.'* The observation of scattering by a 
material having large absorption is difficult 
because the neutrons entering the material are 
absorbed before they can be scattered and it is 
possible that to some extent the failure to observe 
scattering in good absorbers is due to this cause. 
The absence of observed scattering in the region 
of strong absorption is therefore not a surprise, 
particularly in view of the relatively small 
numbers of neutrons available for experimenta- 
tion. It seems more significant, however, that 
strong absorbers do not show, so far, strong 
scattering in any velocity region because, ac- 
cording to the Fermi-Bethe theory, the scattering 
cross section should be large in a wide range of 
energies. The experimental evidence says little 
about the ratio of scattering to absorption near 
resonance. It indicates that this ratio is less than 
1/10 in most cases. It is impossible, therefore, 
to ascertain definitely the ratio [,/T, until more 
detailed experimental data are available. Ac- 
cording to Table I the condition [,/T,<1/10 
can be satisfied in many ways up to velocities of 
over 1000 ev. 


(b) Magnitude of interaction with internal states 
and probability of internal state in required 
region 

In Table I arbitrary assignments of values of 

r,, ', were made. It will be noted that at low 

neutron velocities the desired large capture cross 

sections are easily obtained through relatively 
wide bands having a half-value breadth 2T,. 

Keeping I’, fixed one can decrease the interaction 

energy H to 10,000 ev for hT,=1 volt, E,=1 volt 

and still have a cross section of 1000 10-* cm? 
in an energy range up to 2 volts. In some cases 
relatively weak radiative transitions will come 
into consideration leading to smaller [,. For 
such transitions /7 need not be as large as 

10,000 ev in order to have cross sections of 

1000 X10-** cm? in the resonance region. For 

the large energies involved in nuclear structure 

it is reasonable to expect interaction energies of 
the order of 10,000 volts between practically any 
pair of levels not isolated by a selection rule and 
interaction energies of the order 100,000 volts 
between a great many levels. 

There are about ten elements among 72 ob- 
served that show cross sections of more than 
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500 X 10°** cm?®. Allowing for the fact that there 
are more isotopes than elements it appears fair 
to say that the chance of such an anomalous 
cross section is about 1/20. One can try to 
account for these solely by the low velocity 
regions which exist for any resonance level, thus 
probably overestimating the necessary number 
of levels. In order that o.>50010-*4 cm? at 
1/40 volt for a nucleus having r7>=10-" cm and 
hY,=10 volts the resonance region must be not 
farther than at |E#,|=H/420 from thermal 
energies by Eq. (35). We do not wish AT, at 
thermal energies to be greater than 0.1 volt so 
as not to have too much scattering and therefore 
H should be below 210° ev at the higher £,. 
Thus £, should be kept below about 460 volts 
in order to give the large capture cross sections 
for E=1/40 volt together with small scattering. 
A level below ionization will also be effective in 
producing an increased absorption. The observed 
number of large absorptions corresponds in this 
way to one level in 900 volts for 1/20 of nuclei 
or one level every 18,000 volts for a single 
nucleus. In addition some cross sections will be 
caused by direct resonance. Just how many is 
uncertain but it is clear that such effects exist 
in Cd, Ag, Au, Rh, In. 

The average spacing between the 4 
of Th C” as given in Gamow’'s book is about 
100,000 volts and this is apparently the order of 
magnitude usual for y-ray levels of radioactive 
nuclei. There appears to be no reason why the 
energy levels found through the analysis of y-ray 
spectra should include all the nuclear levels and 
there may be as many as one level in 20,000 of 
a kind that may be responsible for coupling to 
incident neutrons. It should be remembered here 
that some of the levels may be active even though 
the coupling is weak so that more possibilities 
are likely to matter than for the y-rays of 
radioactive nuclei. 

For a complicated configuration of particles 
it seems reasonable to consider a total number 
of 100 possible levels per configuration because 
protons and neutrons can be combined separately 
to give different states. On this basis we deal 
with an average spacing between configurations 
of about 2 MEV which is not excessively small. 
It is, of course, impossible to prove anything 
definitely without calculating the levels; this 


ray levels 
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appears to be premature at present on account 
of uncertainties in nuclear theories. 


(c) Existence of two maxima 

According to the calculation given above it is 
expected that there will be in general two 
maxima one of which should be at resonance and 
another at v=0. According to the experiments of 
Rasetti, Segré, Fink, Dunning and Pegram?" the 
1/v law is not obeyed by Cd but is obeyed by 
Ag. Cadmium has therefore a resonance region 
close to thermal velocities. In the classification 
of Fermi and Amaldi* this region must be 
affected by the C group since absorption measure- 
ments by the Li ionization chamber which was 
used in these experiments agree for most elements 
with the measurements of Fermi and Amaldi on 
the C group.” The verification of the 1/v law 
for Ag by the rotating wheel indicates that in Ag 
the resonance band is located above thermal 
energies. This conclusion is in agreement with 
the smallness of the temperature effect for the A 
neutrons detected by silver which was recently 
established by Rasetti and Fink.** Since Rh 
behaves similarly to Ag in these temperature 
experiments Rh also has a resonance region 
above thermal velocities. Fermi and Amaldi 
have evidence that D neutrons, which affect Rh, 
are different from A neutrons which affect Ag 
It is very probable that both of these groups lie 
above the thermal region and they may reason- 
ably cover a range of 30 volts inasmuch as the 
B group overlaps weakly with both A and D. 

According to Szilard® In shows strong selective 
effects outside the C group and according to 
Fermi and Amaldi*' the same period of In (54 
min.) detects the D group. The number of 
neutrons in the groups is presumably in the 
ratios C/80=B/20=D/15=A/1. One could try 
to conclude that the order of increasing energies 
is C, B, D, A on the assumption that the number 
of neutrons increases towards low energies. Such 
a conclusion is dangerous because little is known 
about the velocity distribution, because within 

20 F, Rasetti, E. Segré, G. Fink, J. R. Dunning and G. B. 
Pegram, Phys. Rev. 49, 103 (1936). 

21 E, Amaldiand E. Fermi, Ricerca scientifica 2, 9 (1936); 
E. Fermi and E. Amaldi, Recerca scientifica 2, 1 (1936). 

2 Unpublished results of F. Rasetti. We are very grateful 
to Professor Rasetti for informing us of these results. 


23 F, Rasetti and George A. Fink, Bull. Am. Phys. Soc. 
11, Paper 28, Feb. 4, 1936. 
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each group there may be several bands at 
different velocities, and also because the number 
of expected neutrons in a group should depend 
on its width. Temperature effects show that 
practically all captures increase as the energy 
is lowered. The effects are strongest™* for Cu, V 
are smaller for Ag, Dy weaker for Rh and 
weakest for I. The absorption coefficient for C 
neutrons is, however, larger for Rh than for Ag 
indicating that the smaller temperature effect in 
Rh is due to a relatively greater importance in 
it of a band above thermal energies. All temper- 
ature effects agree in indicating the presence of 
a region in which the 1/v law is followed approxi- 
mately but again no definite conclusion about 
the order of bands is possible. The low tempera- 
ture effect in I would tend to indicate that its 
absorption region is high and detection-absorp- 
tion experiments on I and Br tend to indicate 
that their bands are isolated from the others 
discussed here; perhaps these isotopes have 
resonance bands at higher energies. A new band 
was recently discovered in Au _ by Frisch, 
Hevesy and McKay*® which represents strong 
absorption on a weaker background. The large 
number of selective effects observed makes the 
present explanation reasonable and the existence 
of a region of low energies in which the absorption 
decreases with energy is seen to fit in well with 
expectation. 


(d) Other possibilities 


One may consider weak long range forces as a 
possible explanation of the same phenomenon. 
Potentials of the order of neutron energies in a 
region comparable with the neutron wave-length 
would produce strong effects on absorption and 
scattering. For thermal energies the wave-length 
is of atomic dimensions and one would therefore 
expect the binding energy of deuterium com- 
pounds to be different from that of hydrogen 
compounds by an amount comparable to 1/40 
volt if such potentials were present. Such energy 
differences do not exist. It would be possible to 
devise potentials which fall off sufficiently rapidly 
with distance to make the interaction potential 
negligible for chemical binding and which would 

*P. B. Moon and R. R. Tillman, Proc. Roy. Soc. A153, 
476 (1936). 


80. R. Frisch, G. Hevesy and H. A. C. McKay, 
Nature 137, 149 (1936). 
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cover a total region appreciably larger than the 
nucleus. Such hypotheses seem improbable with- 
out additional argument. Besides special rela- 
tions between the phase integrals through the 
nuclear interior and the part of the range of force 
outside the nucleus would have to be set up in 
order to make absorption large and scattering 
small. It is improbable that the large number of 
bands could be accounted for by any single 
particle picture 

Forces between electrons and neutrons even 
though they may exist are not likely to have 
much to do with the.bands Thus it has been 
shown by Condon” that electron neutron inter- 
actions would give rise to scattering cross sections 
varying roughly as the square of the atomic 
number Z on the assumption that the electron- 
neutron forces alone are responsible for the 
scattering. Forces inside the nucleus must also 
be supposed to contribute to the phase shifts 
responsible for scattering. Since these forces also 
vary with Z one could obtain a more complicated 
dependence of the scattering cross section by 
suitably adjusting the nucleus-neutron and 
electron-neutron potentials. On such a picture 
one could try to account for sharp resonances by 
making the electron neutron interaction re- 
pulsive. However, Condon’s calculation shows 
that isotope shifts would be also produced by 
these interactions. It is improbable that the 
isotope shift is due solely to neutron-electron 
interaction because the deviation from the in- 
verse square law inside the nucleus due to 
smearing out of protons produces a considerably 
larger effect than the observed shift. But it 
would also be unreasonable to try to combine 
the proton and neutron effects in the nucleus so 
as to have each large but their difference small. 
It is therefore probable that the electron-neutron 
interaction is not much larger than that which 
corresponds to the observed isotope shift. Since 
the density of the Fermi-Thomas distribution 
varies for small 7 as r~! the effective potential 
acting on the neutron will become high for small 
r. However, calculation shows that it is not 
high in a wide enough region to account for 
sharp resonances if the limitation due to the 
isotope shift is considered. 





** E. U. Condon, in press. We are indebted to Professor 
Condon for showing us his manuscript before publication. 
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Bombardment of light nuclei with charged par- 
ticles has also shown the existence of resonances. 
Thus there are resonances*’ for the emission of 
y-rays in proton bombardment of Li, C, F and 
similarly there are the well known resonances in 
disintegrations produced by a particles. Experi- 
mental methods have not been very suitable so 
far for the detection of resonance regions on 
account of the scarcity of monochromatic sources 
and the necessity of using thin films. In Li 
protons are apparently able to produce y-rays 
in two ways; by resonance at 450 kv and by 
another process at higher energies. In fluorine 
there are several peaks. In carbon there was an 
indication of the main resonance peak being 
double. It appears possible that many more 
levels will be detected inasmuch as neutron 
experiments indicate a high density of levels. 
Calculations on the radiative capture of carbon 
under proton bombardment" lead to a higher 
vield than is observed by a factor of several 
thousand. In these calculations the capture was 
supposed to occur by a jump from the # state 
of the incident wave to an s state of the N"® 
nucleus. The calculated half-value breadth due 
to proton escape from the quasi-stationary p 
level was of the order AI'ry~10,000 ev and thus 
much larger than the width due to radiation 
damping. The yield in thick targets under these 
conditions is nearly independent of the special 
value of AI’y. It is clear from the formulas given 
here for neutron capture that one can decrease 
the theoretically expected yield either by ascrib- 
ing the capture to a transition having a small 
probability of radiation (small I’,) such as would 
forbidden 


correspond to quadruple or other 


transitions or else by using an intermediate state 
App! 


Variation of damping constant with energy and Dirac’s 
frequency shift 


Eq. (6) of the text lead to [cf. Eqs. (126’) to (129’) of 
reference in footnote (12) ] 
(vso— votty)(y—v +17) = |A,o\” 
a _il—e** é 
+(vso—vot+ty)[Z’| As |? 
vs —votty (38) 
1 —e?* "9 ¥ iy 
+ ='B, 
Ve—vottly 


27L. R. Hafstad and M. A. Tuve, Phys. Rev. 48, 306 
(1935); P. Savel, Comptes rendus 198, 1404 (1934), Ann. 
de physique 4, 88 (1935). 
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of excitation of the nucleus with a small transi- 
tion probability to the incident state of the 
proton (small I’,). In the latter case this transi- 
tion probability would have to be made so small 
as to have [,<I, and the observed width of 
resonance would have to be ascribed to experi- 
mental effects. If [,<I, the thick target vield 
depends on I, and is proportional to it for 
small I',. The apparent disagreement between 
theory and experiment previously found for 
carbon is thus not alarming from the many-body 
point of the belief that 


excitation states of the nucleus have often to do 


view and supports 
with the simultaneous excitation of more than 
one particle. 

The excitation states responsible for the neu- 
tron absorption bands make it possible for a fast 
neutron to lose energy by inelastic impact with 
the Estimates that the 
sections for such processes are likely to be small 


nucleus. show cross 
when energy losses are high. The cross section is 


estimated to be 
Ay Hr? 
4iaAy U* cos' Kry 

where A, A, are, respectively, neutron wave- 
lengths in the incident and final states. For large 
energy losses As> A, and only a small effect need 
be expected. The excitation levels responsible 
for neutron capture will give small values A; / Ao. 
Excitation levels located lower are more favor- 
able and probably the excitation of Pb to about 
1.5 MEV has to do with such a possibility.** 

We are very grateful to Professors R. Laden- 
burg and F. Rasetti for interesting discussions 
of the experimental material. 


NDIX I 


which determines [ by comparison with (7) [T!>y]. 

It is by no means natural that this equation can be 
satisfied because the right side depends on ¢. If the A, as 
well as the B, were all essentially equal and if y were great 
in comparison with the frequency differences of consecutive 
levels the sums could be transformed into integrals in the 
well-known that (9) as well as (6) would 
follow. We shall attempt here a more exact procedure. 

Consider the =,’ in the square brackets. It is natural to 
divide the range of », into two parts: one for which 
>a>y7 Since 


way'®-!? so 


Vy % and one for which |v,—vo Sa. 


% J. Chadwick and M. Goldhaber, Proc. Roy. Soc. A151, 
479 (1935). 
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A,,*, Av, changes slowly this quantity will be replaced by 
a constant in |»,—vo| a and its value may be taken to 
be that at vo for the evaluation of the contribution of this 
region. We have then to consider 

A.|*%.. J Sona Daten 39a) 
1 Vs—VoTty 

An exact evaluation of this sum is not simple because 7 and 
Av are of the same order of magnitude and the replacement 
of (39a) by an integral is somewhat objectionable. This 
point has never been completely cleared up and we have 
only qualitative arguments in favor of the correctness of 
the replacement of (39a) by an integral. For v,—vo of the 
order of a few y such a replacement is indeed meaningless 
but fortunately this region is not vital for <1/y since 
the numerator of (39a) is then small. For larger | »,—v» 
the terms of (39a) vary more smoothly and finally they 
become rapidly oscillating for | »,—vo|/t>>1 which can be 
satisfied simultaneously with t<1,/y provided a>y. The 
smallness of y is thus not as serious as might appear from 
the fact that y/Av~1. It should also be observed that the 
treatment of Rice® using real eigenwerte for a single one- 
dimensional continuum is in agreement with replacing 
(39a) by an integral. The result of doing so is given by (9), 
In addition one has the contribution of |v,—v9!| >a. 
This integral can be treated neglecting y because it is of 
interest to‘evaluate y only to quantities of order y/(v—v») 
and because the discussion is supposed to apply only to 

yt<1. This integration gives 
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| - A, |? we 39b) 


which means that the principal value of the / is under- 
stood. Similarly one obtains a contribution due to |B, 
These two integrals give the Dirac frequency shift which 
is included in Eq. (15). 

As stated in the text the difference bet ween I" and T' does 
not affect the absorption for ¢>>1/Il'. Thus for the initial 
condition given by Eq. (4) 


v— vo) +12] n vo)? + 
1+enttl t¢— 2¢-28lt cos 2x(: vol 
+ . > +cross product term 7° 
vp—v)?+r" | 


Only the first fraction in the curly brackets contributes to 
the steady increase of =| 6, 2 in times >>1/T.. Its contribu 


tion is 


[ y—vy)?+I? JA», 


The last factor is for practical purposes 4ry/. rhe second 
and third terms in the curly bracket give terms exp( —47I"'t), 
exp(—2zIt) and constants. The first two kinds die off and 
the last kind represents the effect of transients which do not 
matter in the long run, so that for times not too large as 
compared with 1/y and yet great as compared with 1,1 


one may consider the rates of change of 2/4, /? and of 
D'\a,\2 to be 4ryI,,/ Tl and 4ryI,/T. These are the results 


used in the text. 
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The Theoretical Cross Section for K Electron Ionization by Electron Impact 
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Born’s approximation is applied to the general problem of single K shell ionization without 


excitation. Slater’s wave functions are used and a formula is derived which expresses the cross 


section as a function of the atomic number, the ionization potential of the A shell and the 


energy of the colliding electron. By comparison with experimental measurements the theory is 
shown to give good qualitative results for collision energies above about three times the 


ionization energy. 


HE Born approximation leads, in the case 

of single ionization by electron impact, to 
a value of the collision cross section given by 
the relationship 


o(k, c, J) =(m?, 2xh'*k) | k’ | UPR dk" dada’, (1) 


‘ad 
where U= | W,e'% fo VY *e- ik’ tog" tidy, (2) 


Here m is the mass of the electron, &, k’ and k”’ 
are equal respectively to 2rm/h times the 
velocity of the colliding, deflected and ejected 
electrons. VY, is the wave function of the normal 
atom and W; that for the ionized atom. 7 and 7, 
are, respectively, the radial coordinates of the 
colliding and ejected electrons. dQ’ =sin 6'dé’d¢’ 
and dQ” =sin 6’d6’'"do"" where the angular co- 
ordinates are those of the vectors k’ and k”. 
The origin of the coordinate systems is chosen 





. ws 


Jt 
we 
to 


at the atomic nucleus leaving the polar axis to 
be designated later. The integral in Eq. (2) is 
extended over the space coordinates of all the 
electrons. 

The perturbation energy is 


where 7,0= |f2—Tfo 

Let us build the functions representing the 
normal and the ionized atom from which a 1s 
electron has been removed by a product of 
functions representing individual electrons. 


/ 


vY.= Iy, and 
l=] n 


Y; = 


II Ym’. 


The y's will be selected from the appropriate 
functions of Slater.'! The primed functions will 
be identical with the unprimed in the angular 
factors but their radial factors will differ from 
those of the unprimed by a suitable increase in 
the effective nuclear charge. 

Following Bethe’s procedure® we may integrate 
Eq. (2) over the coordinates of the colliding 


electron. Thus 


a eat 4r 
| —divy) = —e'a'™, 
3 


Tnod q~ 
where q= k—k’ gives 


U=(47e? ¢) | VV te © (s—De'a™) dv’. (3) 


e n=1 


Since q is not a function of the électron coordi- 
nates we may select it as their polar axis. Then 
q:t,=9r, cos 6, and we may expand the expo- 


nential as 


e* arn cos On = 1 +197, cos 6,—-***. (4) 

The important contribution to the matrix 
element occurs when the scattering angle is 
small, i.e., q near its minimum value, and when 
r, is not much greater than the radius of the 
Bohr orbit occupied by the mth electron. Under 
these conditions gr, is less than unity and we 


1]. C. Slater, Phys. Rev. 42, 33 (1932). 
2H. Bethe, Ann. d. Physik 5, 325 (1930) 
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may as a first approximation neglect terms of 
the expansion beyond the second. This approxi- 
mation becomes more exact with increase of 
collision energy. 
Combining Eqs. (3) and (4) gives 
trie? zs 
U=- ¥ | V.Vie*' r,, cos 6,dv’. 5 
q n=le 


Because Slater’s functions are of the hydro- 


genic type and because ionization does not 
affect the angular factors there will arise in the 
integrals of Eq. (5) for all values of m other than 


unity a factor of the form 
° 


| cos 6,|P,'(cos 6,)|* sin 0,dé@,. 


e 


It is possible to show by repeated application of 
a reduction formula in a process too lengthy to 
be of interest here that the above integral is 
equal to zero for all values of s and ¢. 

Eq. (5) is therefore reduced to one term of the 
sum, that is where n=1. 

4rie* 7 
U= ——— | WV ,*e7 ©" tir, cos 0,d0". 
gq e 

It may be shown that changing the effective 
nuclear charge by unity does not affect greatly 
the normality of the ¥-functions we are using. 
That is 


| vey2"'dv,~ 1. 


Making this assumption and integrating over 


all coordinates except those of the ejected 
electron gives 
4rie* 
U=-— “j yie®"’-" cos 8; 
q 
Xsin 6,7;°dr,d6\d¢y, 6) 
where ¥i=(a*/r) te, a=[(z—0.30)/a0], 


do is the radius of the first Bohr orbit in hydrogen 
and z is the atomic number. 

Let « be the angle between q the polar axis 
and k’’. Then, neglecting subscripts since the 
remaining coordinates are all those of the 
ejected electron, we have 
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r-k’’=rk'’(cos @ cos e+sin 6 cos ¢ sin e), 


where ¢ is measured from the plane containing 
q and k’’. Eq. (6) may now be integrated over ¢* 


o-7 
e k’’r sin € sin 8 cos dd 
. u ) ” . . 

=2rJo(—k’’r sine sin 6). 
The @ integration becomes 
ot 

2x | e ik’’r cos €cos 8 J)( — p!’y sine sin A) 

al 


Xcos 6 sin 6dé, 


which is a form of Gegenbauer’s integral having 
the value 


(27)! J3,0( —k’’r). 


The Bessel function may be expressed in trigo- 
nometric form and the final integration is 


woo sin k’’r 
. ; , 9 r 
— (427 cos €/k’ 7 ree -~—cos k’’r }dr. 
7 
e k"'r 


0 
This may be integrated by parts to give 
— 32rk"ia cose/(a?+k’”’)3, 
We have finally 
Rk” COS € 
U=A , where A = —2?7 99/2 q5/2¢@2, 


q (a?+k” )* 


Substituting this value of U in Eq. (1) and 
expressing sin 6’d6’=qdq/kk’ from the relation- 
ship g@=k?+k’ —2kk’ cos 6’ we obtain 


m*A? pk’ cos? «¢ 
o(k, 2, J)= =| -- 
2rhtk. q (a?+hk’’*)6 
<sin 6'd0"dg''dg'dqdk”’. 
Since k and k’ remain fixed during integration 
over 6” we are at liberty to choose q as the polar 


axis from which 6” is measured. Thus 6” =«. 
Integrating over 6’’, ¢’ and $” gives 





*H. S. W. Massey and C. B. O. Mohr, Proc. Rov. Soc. 
A132, 612 (1931). 
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4arm*A? ok” max. > Imax. 
o(k, 2, [)= | 
3htk? Jo - 
dq as 
x dk". 


q (a?+k”’ . 


Performing the integration over g and making 
the substitution of V and V” for k and k” from 
the equations 


k?=8r?meV/300h?, 


Rp” - Sr°me ] hdd 300h?, 


we have 


a( V, 2, V;)= | log 
2). 


(a? +k” i 


where V; is the ionization potential in volts of 
the 1s electron, V the potential through which 
the colliding electron has dropped before en- 
countering the atom and V” is the potential 
drop necessary to stop the ejected electron. The 
energy relationship are V=V'+V"+V;, and 
Rak? +k’?+J2. [2=822meV;/300h. 

Eq. (7) gives the ionization cross section for 
collisions in which the 1s electron is ejected from 
the atom. The cross section is seen to depend 
only upon the energy of the colliding electron 
and two constants of the atom, the ionization 
potential and the atomic number. The final 
integration must be performed numerically. 


COMPARISON WITH EXPERIMENT 


Experimental cross sections for the ionization 
of K shell electrons have been measured for only 
the two atoms helium and silver. Smith‘ and 
Liska® have done the last work on helium and 
their absolute values are used in this comparison. 

Fig. 1 shows the measured‘ and calculated 
values (from Eq. (7)) for low energy collisions in 
helium. For values of V>4V; the qualitative 
agreement is good. It will be noticed particularly 
that the maxima of the two curves agree well. 





4P. T. Smith, Phys. Rev. 36, 1293 (1930). 
5 John W. Liska, Phys. Rev. 46, 169 (1934). 
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Fig. 1. The cross section for ionization of helium by 


electron impact. The circles are calculated values and the 


solid curve represents the experimental values determined 
by P. T. Smith 
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Fic. 2. The cross section for ionization of helium by 
electron impact. The circles are calculated values and the 
solid curve represents experimental values determined by 


John W. Liska. 


Quantitatively, however, the error is rather large 
and agreement at these low energies is not as 
good as that obtained in an earlier calculation® 
where the expansion of Eq. (4) was not used. 
This is to be expected from the nature of the 
assumption made in neglecting terms of Eq. (4). 

Fig. 2 shows a comparison for large collision 
energies in helium. Again good qualitative agree- 
ment is obtained but the calculated results are 


®W. W. Wetzel, Phys. Rev. 44, 25 (1933). 











Fic. 3. The cross section for K shell ionization of silver 
by electron impact. The circles are calculated values, the 
solid curve represents a prediction of what the absolute 
value of Webster’s measurements should be and the crosses 
represent the absolute value as fixed experimentally by 


Clark 


too high by varying percentages of the total 
cross section. 

Fig. 3 shows the values for the cross section 
for 1s ionization of silver. Again the calculated 
and experimental maxima check and the theory 
seems to be approaching the experimental values 
at high collision energies. 

The predicted curve was determined by cor- 
recting the relative experimental values given by 
Webster and his associates.’ Since this correction 
involved an unjustified assumption, namely, 
that the percent error between the calculated 
and observed curves at a point near their 
maxima should be the same in the case of helium 
and silver, the publication of this paper was 
delayed. Recently Clark* has made a determi- 
nation of the absolute values of the cross section. 
The results of adjusting the experimental points 
according to his data are shown as crosses on 
Fig. 3. The predicted curve agrees with Clark's 
measurements within the limits of his probable 
error. This gives some justification for using 
Eq. (7) corrected in this fashion to determine 
absolute values of 1s ionization for atoms 
intermediate between helium and silver. 


7D. L. Webster, W. W. Hansen and F. B. Duveneck, 
Phys. Rev. 43, 839 (1933). 

8]. C. Clark, Phys. Rev. 48, 30 (1935). References are 
given in this paper to other theoretical treatments of the 
problem. 
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Vibrations of Tetrahedral Pentatomic Molecules 


Part V. General Criteria for Potential Functions 


* 


Jenny E. ROSENTHAL, 


Chemistry Department, Columbia University 


Received February 17, 1936 


The requirement that the forces along the bonds of the 
YX, molecule be fairly large as compared with “interaction 
forces’ is shown to be sufficient to remove the ambiguity in 
the values of the constants in the potential energy function; 
no further special assumptions as to the types of forces 


involved are necessary. The ratio of the masses of the Y 


[' has been pointed out in the earlier papers' 
f this series that the four fundamental 
frequencies of the tetrahedral molecule YX, 
enable us to calculate directly two of the five 
constants in the potential energy function and 
to express two of the others as functions of the 
third one considered as a variable parameter. 
While the values of some of the constants are 
thus undetermined, certain limits are imposed 
on them by the requirement that they should 
all be real. If the rotational structure is known, 
the absolute value of the parameter D may be 
found from the vibration-rotation interaction.’ 
However, even in the latter case there remains 
an ambiguity, for corresponding to every value 
of the parameter there are two possible sets of 
values for the other constants. This is due to 
the fact that the constants A and B, or rather 
A,/4u and 8B are obtained as the roots of a 
quadratic equation. (u=1/,(4m+M) and M 
and m are, respectively, the masses of the Y and 
X atoms.) There are no mathematical grounds 
for deciding which of the two quantities corre- 
sponds to the larger root, so that both possi- 
bilities have to be taken into account. Since 
only the square of the parameter D may be 
obtained, the one linear and two quadratic 
empirically known functions of the three con- 
stants A, B, and D give four different sets of 
solutions. These again lead to four different 
sets of vibration frequencies for the isotopic 





*Sarah Berliner Research Fellow of the American 

Association of University Women. 
wet II, Phys. Rev. 45, 538 (1934) and Parts III, 

, Phys. Rev. 46, 730 | 1934). The notation in the present 
paper conforms to the one introduced there. 

*M. Johnston and D. M. Dennison, Phys. Rev. 48, 868 
(1935). Cf. also D. M. Dennison and M. Johnston, Phys. 
Rev. 47, 93 (1935). 


and X atoms determines in practically every case which set 
of force constants is the chemically acceptable one. 
Limiting values are given for the constant of the YX- bond 
for CH,, CCly, TiCly, SO4", and other molecules. Wherever 
possible comparison is made between the limiting and the 


actual values. 


molecules. If we introduce some simplifying 
assumption, for example, the Urey-Bradley 
field,* we can discard all but one set of constants 
for chemical reasons. However, the validity of 
such an assumption may be questioned, in 
particular since it involves a smaller number of 
constants than are permitted by symmetry 
considerations.‘ A method will now be discussed 
which makes it possible to set up a very general 
criterion for discarding certain sets of constants. 

If we write the potential function in the same 
form as has been done for the XYZ; type 
molecules,> we can ascribe a physical meaning 
to the constants involved and make definite 
predictions as to the relative orders of magnitude. 
Since this particular form of the potential energy 
has been given before, only the terms essential 
for the further discussion will be repeated here. 
Let 67; as before denote the displacement 
between the Y and the ith X atom and let 6g;, 
be a similar quantity for two X atoms. The 
potential energy is' 


3H. C. Urey and C. A. Bradley, Jr., Phys. Rev. 38, 1969 
1931). 

*The only way to check this assumption would be to 
compare the values of the constants A, B, and D calculated 
by means of it with those obtained with the use of addi- 
tional experimental data, e.g., isotopic shifts, vibration- 
rotation interaction, or the spectra of certain derivatives of 
YX, like XYZ; and ZYX;. T his comparison has so far been 
possible for two molecules. For CH, the agreement is very 
good, as a matter of fact the potential function of CH, is of 
an even more special type, namely the valence force type. 
For CCl, the force field does not appear to be of the 
“7. Bradley type. 

H. Voge and J. 
137 (1936). 

®It should be pointed out that while the constants Gy, 
etc., are susceptible of a simpler physical interpretation 
than A, B, etc., the latter correspond to geometrical 
symmetry coordinates. As written in Eq. (1) the potential 
energy function would prove very cumbersome for any 
actual calculations. 


Rosenthal, J. Chem. Phys. 4, 
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2 V =G,,(6r1 + 6r2+ 673)? +Go26r4" 
+ 2G126r4( 571+ d5r2+ drs) 
+ 2G23674(dgi2+ dqist+ ges) +-**, (1) 


where Gi, Goo, Giz, Go, are related to the usual 
constants A, B, D and E by: 


Gi = At+4B+4E—3D, 
i229 = fg A +4B+4E+D, 
Gi.=4B+4E— 3A —D, 


(14 -4B41D)y/3. 


V“_c; 
“2 
II 


Solving for A, B, and D we get: 
A =G,,;+Ge—2G 2, 
B=} (3G22+§Gii— Giz) — 3Gesy 3 (3) 


D=3G22.—3Gi1— }Gie. 

Gx is the constant of the YX bond, and Gy is a 
quantity of a similar nature, while Gjz and Go; 
are interaction constants. Simple ideas con- 
cerning chemical structure lead us to expect Gi 
and Gz: to be positive and of the same order of 
magnitude with Gi. and G23 comparatively small 
in absolute value, particularly Gee. This assump- 
tion which has been verified for the two molecules 
CH, and CCI, where the calculation was possible’ 
seems to be a very reasonable one. We cannot, 
of course, predict, how close the values of Gi 
and Gs: should be for any given molecule, or 
how small the interaction constants. The last 
question could not be answered anyway, for 
while G2: is quite insensitive to small changes in 
the frequencies, Gi and G2; are very much 
affected. Thus for CH, the constants A, B, and 
D obtained from the observed CH;D frequency 
2205 lead to Gy2=0.04310° dynes/cm and 
Go=0.002 as compared with G».=0.442 and 
Ge3= —0.250 obtained with the use of the 
vibration-rotation interaction data. 

On the basis of this assumption concerning the 
relative orders of magnitude of the various 
constants G we may conclude from Eqs. (3) that 
D>0O. To find out, whether A/4u is smaller or 
larger than 8B, we compare their values ex- 
pressed in terms of the G’s. If A/4u>8B, then: 


7 See Table II in reference 5. 
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G),(1+ M /8m)+Geo(1—3M, 8m) 
— 2Gi2(1 — M/8m) > (M/m)Go/6 4 


This inequality should certainly be satisfied jf 
3M/8m<1, but not if 3M/8m>2. No conclu- 
sions can be made at present for the intermediate 
region, i.e., for 1<3M,8m<2. However, with 
the exception of SnCl, for which 3.8m = 1.255, 
there are no molecules with known vibrations 
for which the mass ratio has this intermediate 
value. Thus we conclude that for practically all 
molecules except the hydrogen compounds, the 
larger root of the quadratic equation represents 
A/4y. This result differs markedly from the one 
obtained on the basis of the Urey-Bradley field 
where the smaller root has always to be chosen 
for A/4u. It might be interesting to point out in 
this connection that in those cases where the 
valency assumption holds, we have to select the 
roots so that A/16>B, i.e., A/4u>8B if 2u<1, 
or 3M/8m<3/2. We shall now apply these 
results to find the limiting value for the bond 
constant Gop. 

Consider first the case 3M/8m <1. If we make 
use of the formulas giving A and B in terms of 
the parameter D* and substitute in Eq. (2) for 
Goo, we get: 


Goo m =%(2+4)(As+As) —$(2—w)[ (As —Ag)? 
—8D?/um? }}+5+D/m. (5 


Here \=4r*w*. Since D can only vary from 0 to 
sm(A3—Ay)\/(u/2) (the last condition being im- 
posed by the reality requirement), we have: 


tri tamArs+ dAy < Gee m<iry 
+4f1+y (u 2) Pas +h[1—y (u 2) PAs. (6 


Eq. (3) does not lead us to expect D to be 
actually zero, so probably Gy is not too close to 
the lower limit. 


TABLE I. Limiting values of the bond constant of various 
YX .-type molecules. 


Ger MINIMUM Gee MAXIMUM AcTUAL VALUE 
CCl, 2.640 « 10° 6.238 x 10° 3.611 x10 
SiCl, 1.745 4.237 
TicCl, 1.292 3.085 
SnBr, .996 2.311 
CH, 4.127 4.944 4.916 
SO, 5.023 8.328 
ClO, 4.986 7.081 


All values given in dynes/cm. 


8 Cf. Part III of reference 1. 
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In the case of the hydrogen compounds, the 
formulas for Gee are: 


Goo,'m = § (2+) (As+As) +4 (2 —w)[ (As — Aa)? 
—8D?/um?}i}+5+D/m (7) 


and 


na + 3X; + tury <Gr ms pAit+ 3As+ 4 HMAs, (8) 
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the minimum occurring for D=0 and the maxi- 
mum for D=my(dA:—),4)/(2+u). Here again we 
would not expect Ge to be very close to the 
lower limit. 

The numerical values of the bond constant 
calculated for various molecules are given in 


Table I. 
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The Ferromagnetism of Nickel 


J. C. SLATER, Massachusetts Institute of Technology 
(Received February 11, 1936) 


By using metallic energy levels extrapolated from copper 
to nickel, the energy difference between a nonmagnetic and 
a ferromagnetic state with permanent magnetic moment is 
calculated for nickel, and it is shown that the ferromagnetic 
state is the stable one. Both saturation magnetic moment 
and Curie point are calculated, in agreement with experi- 
ment within the limits of error of the calculation. Extra- 
polation further into the iron group, though less justified 
than to nickel, indicates that ferromagnetism should 


WO different approaches have been made 

to the theory of ferromagnetism, that of 
Heisenberg! and that of Bloch.? The former 
operates with wave functions of the Heitler and 
London type, localized around the atoms of the 
crystal, while the latter uses plane waves travel- 
ing through the whole crystal. Neither theory 
has met with great success, except in the funda- 
mental point of explaining the large energy of 
orientation of the elementary magnets, or 
magnetic spins, within a ferromagnetic crystal 
in terms of exchange energy rather than magnetic 
energy. This fundamental principle is surely 
correct, but Heisenberg’s theory postulates the 
existence of a positive sign for certain exchange 
integrals which, in every case where they have 
been calculated, have negative sign, and Bloch’s 
theory does not permit ferromagnetism except 
for lattice spacings much larger than those ob- 
served in the ferromagnetic metals. Neither 
theory has been applied specifically to the metals 





_' W. Heisenberg, Zeits. f. Physik 49, 619 (1928). See also 

san, Vieck, Electric and Magnetic Susceptibilities, Oxford, 
_?F. Bloch, Zeits. f. Physik 57, 545 (1929). For discus- 
sion of both theories, see A. Sommerfeld and H. Bethe, 
Handbuch der Physik, Vol. 24, second edition. 


persist in that group down approximately to iron. The 
criterion for ferromagnetism previously suggested by the 
author, the existence of inner unfilled electron shells (the 
3d), small in proportion to their distance apart, is justified 
The calculation is not made according to Heisenberg's 
method, which is considered to be unsuitable for application 
to ferromagnetism, except in its general principle of 
explaining the energy of orientation of elementary magnets 
in terms of exchange energy. 


which are known to be ferromagnetic, and neither 
one gives any suggestion of why just iron, cobalt, 
and nickel should have these properties. 

Several years ago the writer® pointed out that 
a necessary condition for ferromagnetism ap- 
peared to be the possession by the atoms in 
question of an incompleted inner shell of elec- 
trons, capable therefore of taking on a magnetic 
moment by proper orientation of the spins, and 
small enough so that the shells of neighboring 
atoms overlapped very little, so that these 
shells did not take part in the cohesion. The 
reason is that the theory of valence and cohesion 
demands, in agreement with experiment, that 
the stable state be that in which the spins are 
neutralized, so as to produce no net magnetic 
moment. The electrons responsible for ferro- 
magnetism cannot then be those engaged in 
cohesion, but must be much more like those in 
free atoms, in which the condition of stability 
is known to be just the opposite: the stable 
state is in general that of highest multiplicity, 
or that with the maximum magnetic moment. 
It was at that time pointed out that the incom- 


3J. C. Slater, Phys. Rev. 36, 57 (1930). 
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pleted groups of d electrons in the transition 
groups were appropriate for ferromagnetism, 
and it was found that of all the elements of these 
groups, those with the greatest distance be- 
tween the d shells in proportion to their size, 
and therefore the least overlapping, were just 
iron, cobalt and nickel. The recent discovery of 
the ferromagnetism of gadolinium, in which the 
4f group is presumably effective, fits in with this 
rule, as Bethe* has pointed out. 

Recent advances in the theory of electronic 
energy levels within a metallic crystal’ have 
made it possible to carry Bloch’s type of theory 
much further than before. In the reference just 
mentioned, it was shown to be possible to 
understand ferromagnetism in terms of this 
type of theory, and to see why those elements 
with non-overlapping unfilled shells should be 
likely to show it. Until the present time, how- 
ever, energy levels for a ferromagnetic crystal 
have not been calculated, and consequently no 
detailed comparison of the theory with experi- 
ment could be made. Such computations have 
still not been carried out, but the calculations of 
Krutter® on energy bands in copper makes pos- 
sible an extrapolation to the next element, 
nickel, which has the same face-centered cubic 
crystal structure, with a good deal of assurance. 
Similar computations are now being made by 
Dr. Manning for iron, and when they are com- 
pleted an interpolation for the whole iron group 
will be possible. In the meantime, however, it 
was thought worth while in the present paper 
to develop the method of calculation, and apply 
it to the ferromagnetism of nickel, using the 
extrapolated energy bands. The result is to 
show definitely that nickel should be ferro- 
magnetic, and to calculate saturation magnetic 
moment and Curie point in good agreement with 
experiment. Thus for the first time the theory of 
ferromagnetism has advanced far enough to 
make one fairly sure of its correctness. At the 
same time, with much less justification, the 
calculations have been extrapolated all the way 
to iron. The extrapolation indicates definitely 
that only a few metals preceding nickel in the 
periodic table should be ferromagnetic; as a 


* See article in Handbuch der Physik, Vol. 24, see p. 596. 
5 J. C. Slater, Rev. Mod. Phys. 6, 209 (1934). 
®°H. Krutter, Phys. Rev. 48, 664 (1935). 
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matter of fact, it would predict the limit of the 
ferromagnetic elements to lie about at iron 
rather than between iron and manganese, but 
the calculation depends on the small difference 
between large quantities, and a small error in the 
extrapolation could shift the limit to its correct 
position. It is hoped to remedy this when the 
calculations for iron are completed. It is worth 
noticing that while the present calculations are 
made by the general method of electrons travers- 
ing the whole crystal, there seem to be good 
reasons to doubt whether any calculations of 
Heisenberg’s type, depending on wave functions 
surrounding the individual atoms, could lead to 
the correct result in any simple way. 

The writer is indebted to Mr. R. H. Fowler, 
Professor J. H. Van Vleck and Professor Francis 
Bitter for illuminating discussions regarding the 
question of ferromagnetism. 


UNFILLED ENERGY BANDS AND 
FERROMAGNETISM 


The wave functions of an electron in a 
periodic potential field representing a crystal 
are modulated waves, behaving around each 
nucleus like an atomic wave function, but 
changing phase from one atom to another. In 
the case of rather tightly bound electrons like 
the 3d electrons concerned in ferromagnetism, 
the wave function is very much like an atomic 
one, falling to a rather low value in the region 
midway between atoms. The energy level is no 
longer a sharp one, as in an atom, but is spread 
out into a band, the spreading being greater 
and greater as the overlapping gets greater. 
Thus in the 3d functions the spreading is not 
very great. In copper, it amounts to about 0.4 
atomic unit, or about 5.5 volts. Since there are 
five types of d electrons (corresponding in the 
atomic case to the five degenerate states given 
by m=2, 1,0, —1, —2), there will be five bands 
in the crystal, but these overlap in a complicated 
way. For our present purposes we need only the 
distribution in energy of the levels. This is given 
in Fig. 1, where f(£), the number of states 
between E and E+dE, divided by dE, is plotted 
as a function of E, for the five 3d bands com- 
bined. The total area under the curve then gives 
the total number of electrons which can be ac- 
commodated in the band ; this is 10, where NV is 
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Fic. 1. Distribution of states in energy, for copper 
f(E), number of levels per unit energy range, plotted 
against F, in atomic units, lowest energies to the left. 3d 
and 4s bands shown separately. Vertical lines indicate the 
portion of bands filled by 1, 2, ... 12 electrons, re- 
spectively. 
the number of atoms in the crystal, since each 
band can hold one of each spin per atom. The 
values of Fig. 1 were found by starting with 
Krutter’s calculations, extending them from 
the particular -directions of propagation for 
which he found them to all directions of propaga- 
tion by a method of interpolation, and tabulating 
the number of states for each range of energy 
of 0.01 atomic unit. The method will be de- 
scribed in a forthcoming paper. The fluctuations 
noted in the curve of Fig. 1 are partly of signifi- 
cance, partly a result of the approximate method 
of calculation used. 

Overlapping the band of 3d electrons is an- 
other band, also shown in Fig. 1, arising from 
the 4s level of the atom, with interaction with 
the 4p and other higher levels. This band, capa- 
ble of holding 2 electrons per atom, has wave 
functions of a very different character, large in 
the regions between atoms, and acting there 
very much like plane waves. Electrons in this 
band act like free electrons, and possess high 
electrical conductivity, in contrast to the 3d’s. 
At the same time they contribute strongly to the 
cohesion, as in Wigner and Seitz’s? calculation of 
cohesion in sodium, and in Fuchs’s* calculation of 
the cohesion in copper, where wave functions 
very similar to the present ones were used. 

Asa first approximation to the structure of the 
metal, we may suppose that levels are filled, 
beginning with the lowest energies, high enough 





"E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933). 
*K. Fuchs, Proc. Roy. Soc. A151, 585 (1935). 
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up so that all the electrons of the crystal are 
accommodated. There will be two electrons per 
state up to a certain energy, none above. In 
Fig. 1, lines are drawn indicating the maximum 
energy if there are 1, 2, 3 --- 12 electrons in the 
3d and 4s bands combined. The case of 11 cor- 
responds to copper, whose atom has 10 3d and 1 
4s electrons. It is seen that the metallic 3d band 
is filled, and that the 4s contains one electron 
per atom. If now the energy levels remained 
the same for the preceding elements, the case of 
10 electrons would represent nickel, 9 cobalt, and 
8 iron, respectively. We shall not try to estimate 
the changes in the energy levels, but shall wait 
until the case of iron is actually calculated. 
We see that the ferromagnetic elements contain 
less than one 4s electron per atom (explaining 
partly the poorer conductivity of these metals 
compared with copper, as Mott® has pointed 
out), but that nickel contains only about 9.5 
electrons in the 3d band, cobalt about 8.5 and 
iron about 7.5. These missing electrons in the 3d 
band are necessary for ferromagnetism. But as 
long as the electron arrangement is as we have 
described, with each state either occupied by 
two electrons of opposite spins or empty, there 
will be no net magnetic moment, and the metal 
cannot show ferromagnetism. 

We must, then, go to a higher approximation 
to find any net magnetic moment. A similar 
situation arises in the structure of atoms. As a 
first approximation, we assign to each electron 
a set of quantum numbers, , /, m, and a spin 
quantum number, and assume that, as in a 
central field of force, the energy of each depends 
on » and /. We then assume that the energy 
levels in the central field are filled far enough 
up to accommodate all the electrons, deducing in 
this way which of the quantum states are oc- 
cupied in the normal state of the atom. As a next 
step, however, we consider the vectorial inter- 
action of electrons, resulting in the multiplet 
structure. The electrostatic interaction is al- 
ready taken into account if we set our central 
held up properly, and the same is true in the 
metal if the periodic field is chosen properly. 
But we must consider the changes of energy 
resulting from orientations of the orbits (change 
in m) and from orientations of spins. The energy 


*N. F. Mott, Proc. Phys. Soc. 47, 571 (1935). 
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of a state proves to be, to the first order of 
perturbation theory, given by the energy neglect- 
ing orientations, minus the sum of certain ex- 
change integrals over all pairs of electrons having 
parallel spins.'° The integral is positive; with 
the result that the energy depends on the ar- 
rangement of spins, and that the minimum 
energy corresponds to the state where as many 
electrons as possible have spins parallel to each 
other. Now suppose there are (10—) electrons 
per atom in the 3d (10—mn)N in the 
whole crystal. If these are divided, as we have 


shells, 


assumed, into equal groups of each spin, there 
will be (5—/2)N of each spin. With each spin 
there will then be 3{(5—2/2)N}{(5—n/2)N—-1} 
pairs, or, neglecting 1 compared with N, the 
total number of exchange integrals will be 
{(5—n/2)N}*. Suppose on the other hand that 
there are only (5—”/2—y/2)N electrons with 
negative spin, but (5—,/2+ y/2)N of positive 
spin, leaving the same number of electrons as 
before, but with a net positive spin of (u/2+ 4/2) 
= electrons per atom. Then there will be ap- 
proximately 


IL +h 


? 


-_ 


2 


\NP+[(5—n/2—p/2)N Ff} 
[(5—n/2)NP+[pN/2} 


exchange integrals, an increase of (uN/2)* in- 
tegrals. Since these terms are to be subtracted, 
there will be a decrease of energy, and a tightness 
of binding, if the spins set themselves parallel, re- 
sulting in a magnetic moment. The maximum 
possible moment is determined by the condition 
that the number of electrons of either spin in the 
3d bands must be between zero and five. The 
second limiting condition is the significant one in 
the ferromagnetic case, and as a result we see that 
the maximum value of u equals n. In this case, all 
states of positive spin are filled, but there are 
nN electrons missing from the states of negative 
spin. The exchange energy in this case is less 
than that in the case of zero spin by (nN/2)? 
times the average exchange integral between 
two electrons. 

We have just seen that the total energy of the 
crystal is decreased, on account of exchange 
effects, if as many spins as possible set themselves 
in permanent magnetic 


parallel, resulting a 


moment. There is, however, a compensating 


10 J. C. Slater, Phys. Rev. 34, 1293 (1929). 
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effect, which may be less or greater than the ex- 
change effect depending on circumstances. If we 
go from the state of balanced spins to that of 
unbalanced spins, we remove uN/2 electrons 
from states of minus spin, and add them to the 
states of plus spin. In the balanced state, the 
2)N levels of each spin were oc- 


lowest (5—m” 
cupied, up to a certain definite energy level 
In the process of unbalancing, the u.N/2 electrons 
of negative spin are removed from energies lowe: 
than this maximum value, and are raised in 
energy so as to occupy levels of positive spin 
above this maximum value. Thus there is a net 
increase of energy. This increase is likewise ap- 
proximately proportional to u*; for the number 
of shifted electrons is proportional to u, and the 
energy difference through which they must be 
raised is, at any rate for small u, approximately 
proportional to u. The net effect of this energy 
increase and the energy decrease due to exchange 
will then be positive or negative depending on 
the relative magnitude of the effects, but in any 
case it will generally remain of the same ,sign 
that 


sO 


independent of the magnitude of ux, 
either the state of balanced spins will be the 
most stable, or it will be the least stable and the 
unbalancing will go the whole way, resulting in 
an unbalanced spin of » electrons per atom. This 
latter case is that for ferromagnetism. 

The magnitude of the exchange integral is not 
very different for different atoms, and for most 
atoms the increase energy on account of 
shifting the electrons to higher levels, increasing 
the binding-energy, is far greater than the de- 
crease of exchange energy when going into a 
magnetic state. (Algebraic rather than absolute 
values of the energies are here referred to.) Thus 
most metals are nonmagnetic. For the ferro- 


in 


magnetic elements we must seek those whose un- 
filled energy bands are as little split up as pos- 
sible, so that as little energy change as possible 
will result from a rearrangement of electrons 
within the band. Thus we are led back to our 
former criterion of nonoverlapping unfilled 
shells, for it is just here that the splitting will be 
small. Furthermore, the whole calculation, 
based on the top part of the band being unfilled, 
the bottom part filled, is one for which the 
method of periodic wave functions is particu- 
larly suited, but for which Heisenberg’s method 
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of atomic wave functions could be adapted only 
with great difficulty. 


CALCULATION OF EXCHANGE INTEGRAL 


To get the exchange effect, we need only find 
the average value of the exchange integral 
between two 3d electrons in the lattice. Such 
an exchange integral between two wave func- 
tions a, and a; is by definition 


Sa *(1)aj*(2)(€/ne)a (1) (2)dr 2, (1) 


where 1 and 2 represent the coordinates of two 
electrons. We divide the volume into N poly- 
hedral cells, each surrounding an atom, and 
integrate separately over each pair of cells. 
Thus the integral is the sum, as p and g range 
from 1 to N, of the integral in which dy, ranges 
over the interior of the pth cell, dv, over the gth 


N NV 
I=>> © exp (i(k;—k,)-(R,—R,)) 
q=1 


p=1 


A single term is the potential of a charge density 
eu ;*(1)u;(1) in the pth cell, on a charge of density 
eu ;*(2)u (2) in the gth cell, times the appropriate 
exponential. First we shall consider the terms 
p=q, which prove to be the leading terms. For 
this case, R, —R,=0, so that the exponential re- 
duces to unity, and the integral becomes sim- 
ply Nfu*(1)u;*(2)(e/re)u (1)u(2)dr2, inte- 
grated through the central cell only. Now within 
this cell, the functions wu; and u; are very much 
like atomic d functions, on account of the rela- 
tively small perturbations of the d electrons of 
different atoms on each other. The normaliza- 
tion, however, is different, each function being 
approximately 1/N! times the corresponding 
atomic function, so that the square of the func- 
tion, integrated over the N cells, can equal unity. 
Thus the exchange integral finally is seen to be 
approximately 1/N times the corresponding 
atomic integral, which we shall call J. 

Next we consider the case p#g. Two different 
cases arise, depending on whether u; and u; 
refer to the same value of the quantum number 
corresponding to the atomic quantum number 
m, or not. If they do, u; and u; are approxi- 
mately equal to each other within the cell, and 
their product integrates to 1/N over the cell, 


cell. Within the pth cell, according to the 
modulation of the wave function described by 
Bloch, the function a; is equal to its value at the 
corresponding point of the central cell, multi- 
plied by the factor exp (7k,;-R,), where k; is the 
vector wave number describing this stationary 
state, R,, is the vector from the center of the 
central cell to the corresponding point of the 
pth cell. Similar relations hold for the other 
functions. Thus our integral becomes 
\ Vv 


> > exp (i(k;—k,)-(R,—R,)) 


— 


p=1 q=1 


times the integral in which dz, ranges over the 
pth cell, dv, over the gth, and where we replace 
a; and a; by their values at corresponding 
points of the central cell. Let these values in the 
central cell be u;, u;; then the exchange integral 


becomes 


a7 


ui*(1)u;*(2)(e?/ry2)uj(1)u;(2) doy. (2) 


so that the charge density eu,*(1)u,;(2) is a 
charge of approximate amount e/N, distributed 
through the pth cell. Its potential at the gth cell 
will then be approximately (e/N)(1/|R,—R,!), 
plus terms in higher inverse powers of |R,—R, |, 
arising from the fact that the charge distribution 
really has multipole components as well as a 
total charge. We shall see, however, that the 
whole effect is very small, and the multipole 
terms can consequently be neglected. On the 
other hand, if u; and uw; refer to different m 
values, the functions will be approximately 
orthogonal to each other, their product will 
integrate almost to zero, over the cell, the charge 
distribution will consist practically only of 
multipoles (beginning with a quadripole), and 
we may neglect its potential completely. We are 
left, then, only with the case of equal m values, 
in which the contribution to the integral is the 
potential of a charge e’N on an equal one at 
distance |R,—R,!. We may then convert the 
sum over p and g in (2) into a single sum in which 
the pth cell is at the origin, obtaining 


q (9) 


= J/N+ N(e®,N*)> (exp (i(k;—k,)-R,)) /R 
“ 


if uw; and u; refer to the same m value, while 
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I=J,/N otherwise. It is to be noted,that J is 
different for different pairs of m values of the two 
wave functions. 

We shall not attempt to evaluate J from (3) 
for each separate combination of functions 4; 
and u;, but only to find the average over all 
possible combinations, recognizing that this 
may involve us in certain errors in the final re- 
sult. In reference 10 it was shown that the ex- 
change integrals between two d functions can be 
written in terms of two parameters, F* and F*, 
which can be evaluated from atomic wave func- 
tions, or from observed atomic spectra. In the 
table of coefficients b*(/m,; l’m,') on p. 1312 of 
that reference, the necessary coefficients are 
given for different combinations of the m’s. 
When we average over the 25 combinations of 
the five m values with each other, we find from 


that table that 
J =2/35(F?+ F‘). (4) 


For nickel, we can estimate the quantities F° 
and F* from the observed spectrum. The 
(3d)*(4s)? configuration of Ni I involves energy 
separations determined by just these quantities. 
Taking the observed multiplet separations for 
this configuration, we can fit F values to it. 
The fit is not very good, but to a rough approxi- 
mation we have F*?=80,000 cm~', F*=50,000 
cm™~'. To indicate the order of accuracy, we give 
in Table I the energy differences between the 
multiplets, as obtained by taking the center of 
gravity of the observed multiplets, and we 
compare with the values determined from these 
parameters. Using these values of the F's, we 
have J =7430 cm approximately. 

Next we wish to average the second term of 
(3) over all pairs of wave functions. We recall 
that this term is present only when u; and u, 
have the same m values, which occurs only 1, 5 of 
the time, so that we wish 1 5 of the average 
over 7 and j of this term. Strictly, the quantities 
k; and k; should be allowed to range over a 


TABLE I. 

OBSERVED COMPUTED 
F-—'D 12,550 cm 13,250 
F—3p 14,750 15,950 
F—'G 21,150 20,750 
IF—1S not obs. 51,600 


SLATER 


polyhedral cell in the & space, but it will make 
only a small error if instead we allow them to 
range over the interior of a sphere of the same 
volume. It is easy to show that the average of 
exp (7k-R) over the interior of a sphere of radius 
K in k space is 3 (sin KR—KR cos KR)/(KR 


Thus the average of the second term of (3 


Ss 


1 3sin KR,—KR, cos KR, . 
: | /R 5 


wont? 
5 q (KR,)? 

Here K is to be so chosen that (4,3)7K® is the 
volume of one cell in & space. If R is the distance 
from an atom to its nearest neighbors in a face 
centered. lattice, it is not hard to show that this 
leads to KR=6!2'x!. Inserting this value, the 
bracket above becomes 0.000352 for the case 
where g is one of the nearest neighbors, and it is 
much less for more distant neighbors. The 
nearest ones are then the only ones which need 
be considered, and since there are 12 of them, the 
expression (5) becomes approximately 


(12/5)(0.000352)(e?/N)/R. (6 


The distance R from an atom to its nearest 
neighbor is about 4.67 atomic units in Ni 
Inserting this value, and converting to cm~, the 
expression (6) becomes 36 cm™'/N. It is thus 
evident that the terms for p#q are almost 
negligible in finding the exchange integral, which 
arises almost entirely from the atomic integral 
Adding our two expressions, we have / =7466 
cm~!, N, and 


Exchange energy of magnetic state —exchange 
energy of nonmagnetic state 
= —(n*/4)(7466)N 
= — 1866 n* cm“ per atom. (7 


On account of the discrepancies between the ob- 
served multiplets and those calculated from our 
F values, we may expect this value to be uncer- 
tain to between five and ten percent. We must 
recall, furthermore, that it represents only an 
average value, and that the exchange energy for 
particular pairs of wave functions can differ 
widely from the average. 


CALCULATION OF CHANGE IN BINDING ENERGY 


We. must next find the change of binding 
energy involved in shifting nN/2 electrons from 
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the states of negative spin to those of positive 
spin, remembering that they must go into states 
previously unoccupied. Since the binding energy 
is known only in the form of numerical calcula- 
tions, as summarized in Fig. 1, this calculation 
must be made numerically, not analytically. In 
terms of the function f(#£) plotted in Fig. 1, 
giving the number of states per unit energy range, 
the total energy of the electrons between two 
energy values, assuming them to fill all energy 
states between the limits, is /Ef(E)dE between 
these limits, or rather a corresponding summa- 
tion, since f(£) has been computed for discrete 
values of E. We must find the energy in two 
cases: (1) the state of no net spin, in which all 
states are filled from the bottom of the band to 
an energy such that there are m2 states per 
atom of each spin vacant above this value; and 
(2) the state of maximum spin, in which all 
states of positive spin are filled, but in which 
states of negative spin are filled only to such an 
energy that there are m states per atom of this 
spin vacant above this value. These energies 
have been computed, as a function of mn, by 
obvious numerical methods, and their differ- 
ence E,— EF, has been found, giving the increase 
in binding energy in going to the ferromagnetic 
state. 


Net ENERGY CHANGE FROM NONMAGNETIC TO 
MAGNETIC STATE 


We have now found that on account of the 
exchange effect, as given in (7), the energy of the 
crystal decreases by 1866 n? cm™! per atom when 
the electrons change from the state of no net 
spin to the state of maximum spin. On the other 
hand, the energy arising from binding increases 
by the amount E.—£,, as discussed in the last 
paragraph, when the same change is made. The 
absolute values of these two quantities, as 
functions of m, are plotted in Fig. 2. It is seen 
that, for small n's, the increase in binding energy 
is less than the decrease in exchange energy, so 
that the net result is a decrease of energy, the 
ferromagnetic state is more stable than the 
nonmagnetic state, and the metal should be 
ferromagnetic. On the other hand, as n increases, 
the two curves are seen to cross, so that for too 
large a value of m (too small a total number of 
electrons in the d shell, or in the atom), the in- 
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Fig. 2. Magnitudes of exchange and binding energy 
changes between ferromagnetic and nonferromagnetic 
states (in cm™! per atom), as function of » (number of net 
electrons spins per atom). The condition for ferromagnetism 
is that the exchange energy term be greater numerically) 
than the binding energy term. 


crease in binding energy becomes greater than 
the compensating decrease in exchange energy, 
so that the state of zero spin is the most stable 
one. The question now is, in which part of the 
curve is nickel located? We have already pointed 
out that, according to our extrapolated copper 
energy bands, nickel should contain about 95. 
electrons per atom in the 3d bands, so that 
n=0.5. Mott® has shown that according to ex- 
periment x is nearer 0.6. This is near enough our 
theoretical results so that the uncertainty in the 
extrapolation could easily account for the differ- 
ence. In either case, Fig. 2 very definitely shows 
that nickel should lie in the ferromagnetic part 
of the curve. It is tempting, as we have remarked 
above, to extend the extrapolation further down 
the iron group, using these same curves for 
cobalt and iron, though of course they are really 
not accurate, since both the exchange energy (as 
computed from the spectrum) and the binding 
energy (as computed from the metallic energy 


bands) will really change as we go from atom to 
atom, and furthermore cobalt and iron have 
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different crystal structures from copper and 
nickel. We have seen that approximately =1.5 
for cobalt, 2.5 for iron. Mott® has shown that 
experimentally »=1.7 for cobalt, 2.2 for iron, 
approximately. These values of » are all obtained 
from the observed saturation magnetic moments 
of the metals. From Fig. 2, on which the values 
1.7, 2.2 are indicated, it is seen that the curves 
cross at approximately the m value for iron, 
instead of between iron and manganese (that is, 
for m>2.2) as they should from the observed 
ferromagnetism of iron. But the error is one 
which we can easily attribute to the errors in our 
extrapolation. Both curves rise rapidly with n, 
and comparatively small changes in them would 
shift the intersection by the required amount. 
We may reasonably say, then, that our theory 
not only predicts the ferromagnetism of nickel, 
but also predicts that the phenomenon of ferro- 
magnetism should persist for some distance down 
into the iron group of elements, approximately 
as far as it is observed experimentally, and then 
should cease. 


MAGNETIC ENERGY AND CURIE POINT OF 
NICKEL 

Let us adopt for nickel the value 0.6 for n, 
as indicated by experiment. Then we find the 
exchange energy to be 1866 (0.6)? cm™!=672 
cm~'. The difference in binding energy between 
the two states is found to be 401 cm~. The net 
decrease in energy on going from the nonmag- 
netic to the ferromagnetic state is thus 271 
cm~! per atom. This of course is at the absolute 
zero of temperature, for with both nonmagnetic 
and ferromagnetic states we have considered 
only the lowest possible energy level, all elec- 
trons being in the lowest state consistent with 
the given total spin. This quantity is not sus- 
ceptible of experimental measurement, but it is 
closely related to the Curie point. For nickel, 
the Curie temperature @ is about 630° absolute, 
and the corresponding energy k@, where & is 
Boltzmann's constant, is 440 cm~', when con- 
verted to spectroscopic units. This is obviously 
of the same order of magnitude as the quantity 
found theoretically. We can, however, make a 
more definite comparison, on the basis of avail- 
able theories of ferromagnetism. For the present, 
we shall not try to replace these theories of the 
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temperature variation of the magnetic moment 
by a more suitable theory based on the present 
model. 

In the Weiss theory of ferromagnetism* in the 
absence of an external magnetic field, one writes 
the internal magnetic field in terms of the 
magnetic moment per unit volume as H= NJ, 
where N is an empirical constant, large com- 
pared with unity, and I is the magnetic moment 
per cc. The magnetic energy per unit volume is 
then —3}H-I=—3NF in the magnetic case, as 
compared with a nonmagnetized state. The 
formula for the Curie point on this theory can 
be written @=o¢no’Np/3mR, where ono is the 
magnetic moment per gram molecular weight 
(equal to our m, times the Bohr magneton, times 
Avogadro’s number), p is the density, m the 
molecular weight, R the gas constant per gram 
molecular weight. In an obvious way oo can 
be written in terms of J, the magnetic moment 
per unit volume, so that @ is plainly proportional 
to the same quantity NJ* appearing in the 
energy per unit volume. In fact, if we multiply 
this magnetic energy by the volume of an atom, 
to get the energy per atom, we have very simply 


k@=? magnetic energy per atom (Weiss). (8) 


We can proceed similarly in the Heisenberg 
theory.t The magnetic energy of the whole 
crystal is then written in terms of certain ex- 
change integrals (not to be confused with our 
exchange integrals) denoted by J, a number of 
nearest neighbors z, and a total spin S’ of the 
crystal, equal to NS, where S is the spin of one 
atom, N the number of atoms in the crystal. 
It is then, as given by Van Vleck’s formula (22 
equal to —zJS”/N. The Curie point, on the 
other hand, is given by Van Vleck’s formula (37 
as 6=2zJS(S+1)/3k. Thus we are led to the 
relation 


k@=2S(S+1)/S times magnetic energy 


per atom (Heisenberg). (9 


It is observed that Heisenberg’s result ap- 
proaches Weiss’s in the limit S=~, as it 
does in other respects. However, as Van Vleck 


* See for example E. C. Stoner, Magnetism and Atomic 
Structure, p. 75. 
+See for example Van Vleck, Electric and Magnetic 
Susceptibilitres, p. 328. 
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points out, the experimental results for nickel 
are best represented by S=1/2, leading to 
25(S+1)/S*=2. Thus we may expect the Curie 
energy to be somewhere between 3 and twice 
the energy of magnetization at the absolute 
zero, with the probabilities favoring the latter 
value. The experimental Curie energy 440 cm~! 
is 1.62 times the magnetization energy 271 cm~! 
which we have found. The agreement seems en- 
tirely satisfactory, considering the nature of the 
approximations which we have made. 


ELECTRONIC SPECIFIC HEAT OF NICKEL AT Low 
TEMPERATURE 

It has been recently observed by Keesom" 
that the electronic contribution to the specific 
heat, which, being proportional to the tempera- 
ture, outweighs the contribution of the lattice 
vibrations at low temperatures, is many times 
the normal amount in nickel. Both Keesom and 
Mott have suggested that this arises from the 
large concentration of energy levels in the neigh- 
borhood of the top of the Fermi distribution, on 
account of the 3d levels. We can easily test this 
hypothesis. Sommerfeld and Bethe (reference 2, 
p. 430) show that the specific heat per gram 
molecular weight on account of free electrons in 


"Keesom and Clark, Physica 2, 513 (1935); W. H. 
Keesom, Proc. Roy. Soc. Al52, 12 (1935); N. F. Mott, 
Proc. Roy. Soc. A152, 42 (1935). | am much indebted to 
Mr. R. H. Fowler for calling my attention to these 
references. 
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a normal metal is C,=(272*,/3)N(¢o)k?7, where 
N(E)dE represents the number of energy levels 
of one spin in dE, if we are dealing with a gram 
molecular weight, and {9 is the maximum energy 
occupied at the absolute zero. Thus his N is half 
our distribution function f. For a ferromagnetic 
metal, however, at low temperatures all levels of 
positive spin are filled, so that they do not con- 
tribute to the specific heat, and we count only 
electrons of negative spin. This removes the 
factor 2 in the formula above. Furthermore, ¢» 
must be taken to be the maximum energy of 
electrons with negative spin, when the electrons 
of positive spin occupy the whole band. In this 
case, from the numerical data from which Fig. 1 
is drawn, we find that N(E£) is 26.15 times 
Avogadro’s number, if energy is expressed in 
atomic units. Substituting this value, and appro- 
priate constants, we find C,=0.00117 calories 
per mole per degree. Keesom’s experimental 
value is 0.0017447. The agreement is not very 
good, but still it is correct in order of magnitude. 
To get better agreement, the peak in Fig. 1 at 
energy —0.45 atomic unit would have to be 
about half again as high as it is, presumably 
being correspondingly more narrow. This does 
not seem impossible, though it is unlikely. In 
any case, the high specific heat is definitely 
connected with the high peak in the 3d distribu- 
tion curve of Fig. 1, verifying again this general 
feature of our curve. 
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A Note on Positron Theory and Proper Energies 


ROBERT SERBER,* University of California, Berkeley 


(Received February 7, 1936) 


The origin of the infinite light-quantum proper energy 
which follows, according to Heisenberg, from the creation 
of matter by the field of the photon, is examined in some 
detail. We are led to investigate the inconsistencies which 
appear on the incorporation of the Dirac positron theory 
into the formalism of the quantum theory. These incon- 
sistencies make it impossible to regard with confidence any 
predictions of the theory for which a consideration of the 


Na field free vacuum, the infinite distribution 
of electrons in negative energy states envis- 
bef t. ; 


5 National Research Fellow. 


singularities of the density matrix is essential. it is shown 
that when the off-diagonal distance in the density matrix 
is taken different from zero, the proper energy of a light 
quantum and the electromagnetic energy of an electron are, 
in order e*, finite, and that the off-diagonal distance plays 
in these calculations the part of a generalized electron 
radius. 


aged in Dirac’s theory of the positron clearly 
should contribute nothing to the expectation 
values of observables of the system (charge and 








546 ROBERT 


current density, energy, electromagnetic fields). 
The contribution of this field free vacuum distri- 
bution has then to be subtracted from the expecta- 
tion value of every observable. Since, however, 
the difference between the contributions of the 
vacuum distribution in the presence and absence 
of an electromagnetic field does not turn out to 
be finite,' the theory in this form cannot give 
sensible and determinate answers to questions 
involving the reaction of the electron distribution 
in varying fields. Dirac himself has now tried to 
modify the formalism in such a way that all 
expectation values would necessarily be finite. 

Dirac’s treatment of the problem? is based on a 
study of the density matrix 


(x! pb” R x’ k’ ‘= dy Ash by n . (x’k’ ). 


Here y,(x’k’’) is a one-electron wave function at a 
“and k” is the spin 


, r 


point xo" =ct’’, x1", x2"", Xs 
variable. The summation is to be extended over 
all occupied states. Dirac has shown, by integra- 
tion of the equation of motion for ®, taking the 
y, as solutions of the ordinary Dirac equation, 
that ® is finite when the length of the four- 
vector x=x’—x” is different from zero, i.e., the 
density matrix is singular only when x,x* =0. The 
worst singularities in ® can be eliminated by 
between 


symmetricizing the density matrix 


positron and electron, that is, by replacing ® by 


(xk | Ri xh’) =D al xR Wnt (a k’) 
— ¥ Walx’R’ Wnt (x’h’) ]. 


Explicit expressions for the singular terms, S, 
which remain in R have been given by Heisen- 
berg.’ These singular terms may be deleted by 
replacing R by r=R-—S. It is of course possible 
also to subtract finite terms. Heisenberg has 
included such a term in S; it is chosen to re- 
normalize the polarization of vacuum for slowly 
varying and weak fields to zero. The expectation 
value of an operator F at a point X= }(x’+x’’) is 


'P. A. M. Dirac, Solvay Congress, 1933; W. H. Furry 
and J. R. Oppenheimer, Phys. Rev. 45, 245 (1934); R. 
Peierls, Proc. Roy. Soc. A146, 420 (1934). 

?P. A. M. Dirac, Proc. Camb. Phil. Soc. 30, 150 (1934). 

3 W. Heisenberg, Zeits. f. Physik 90, 209 (1934). 
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now given by 


F(X) =lim (X —}x, k' Fr X+5x,k 


ro) 


and is guaranteed to be finite and determinate 


The electromagnetic field produced by the 


electron distribution is taken into account by 
including in the total current which appears in 


the Maxwell field equations the expectation value 


of the current calculated from the density matrix, 


jh X)=—elim DO anene(X — bx, k’ |r| X+4y, b). 


z30 kk’ 


Here jo=p, ao=1. Heisenberg has shown that, 
despite the subtraction of the singular terms, the 
conservation law for the current still holds in the 
limit x0, as do also conservation laws for the 
energy and momentum of the entire system, 
material field plus electromagnetic field. 

The passage from the density matrix treat- 
ment, based on the correspondence principle, toa 
quantized field theory can be made by regarding 
the equations of motion (the Dirac equation and 
the Maxwell equations) as matrix equations in 
the variables y, ¥*+, E;:, Ai, which satisfy the 
usual commutation laws. The consistency of such 
a procedure depends upon the compatibility of 
commutation laws and equations of motion 
these equations must, namely, be such that the 
values of the commutators remain unchanged in 
the course of time. This condition is not satisfied 
by the system of equations with subtracta- 
tive terms. Thus, with A;=—E; and E=curl 
H+47ea(R—S), (d/dt)LE(X), A(X’) ]=—4re 
[a:S(X), A(X’) ]#0. In fact the time rate of 
change of this commutator is singular. 

A second scheme has also been proposed by 
Heisenberg. This consists in taking as a Hamil- 
tonian function the energy expression of the 
density-matrix treatment, now regarded as 4 
function of the g-numbers y, ~*+, F;, A;. One takes 
as canonical conjugates y and y*, A; and &,, asin 
the usual quantum electrodynamics. One sets the 
scalar potential Ay=0, and x»=0. The Hamil- 
tonian function is then‘ 


‘We shall employ rational units, measuring length 1 
terms of the Compton wave-length h/mc, time in terms 0! 
h/mc?, and mass in terms of the electronic mass m. In these 
units e? =a. 





Th 


inate. 
\ the 
nt by 
ars in 
value 
latrix, 


x, R). 


that, 
ns, the 
in the 
or the 
ystem, 


treat- 
le, toa 
arding 
on and 
ions in 
fy the 
of such 
ility of 
notion: 
nat the 
nged in 
atisfied 
ytracta- 
E =curl 
= —4re 


rate of 


sed by 
Hamil- 
of the 
d as a 
1e takes 
E, asin 
sets the 
Hamil- 


length 1 
1 terms ol 
. In these 





POSITRON 


} 10 Ox, a Se[ A j| B 4 + 5x) +A 


H=} {aX ¥ Ca; 


THEORY 547 


X —43x)]} +Be-uJWt(X + 4x, k’)W(X — Be, b”) 


—¥(X — bx, kyr (X + 4x, k’) J+ { dX>° 10 / Axo(X —}x, k’'| S|\X+4x,k') +H ew, (1) 


where J7eyw=(1/8r) f'dX(F°+J2). This Hamil- 
tonian function can be formally treated as 
Hermitian, provided one understands that in 
taking the transpose of a matrix one is also to 
reverse the direction of x. The interaction be- 
tween electromagnetic field and matter is to be 
supposed a small perturbation and the usual 
perturbation theory, expanding in powers of the 
electronic charge e, is to be applied. In the result 
of any calculation one is to pass to the limit 
x0. 

However, the Hamiltonian 
with it its own equations of motion, and these are 


function carries 


not the same as the equations of motion of the 


simple density matrix treatment. Thus the 


Hamiltonian equation for y is 
il(X j= [ a; —10/0X ;+ bel} A(X ) 
+A ,(X—x)})+B]y(X—x), (2) 


which differs from the ordinary Dirac equation 
through the appearance of x on the right-hand 
side. Since the Heisenberg subtractative terms 
were determined on the supposition that y obeys 
the Dirac equation, there is no warrant for be- 
lieving them correct when, instead, y obeys (2). 
The Hamiltonian 
magnetic field differ from the Maxwell equations 


equations for the electro- 
even more radically, for the subtractative terms 
introduced into (1) are explicit functions of EF 
and A, and hence give rise to new terms in the 
equations of motion. Thus instead of the ordinary 
equation —A;=E,;, one finds, aside from terms of 
order x, 


~A,= E,—(e* 3m) [xix jE; /x x 

— E; log (|x,x*|/a)). (3) 
The equations for E; can be written 
~E;+curl, H =42a[_ji(x) —e*S,(j) ] 


+e? (singular terms). (4) 


Here j(x) represents the contributions of the 
terms in (1) which are linear in e, and —e*S(j) 
contains the subtractative terms necessary to 
cancel the singularities which appear in the 
expectation value of j(x) in order e®. The appear- 
ance of singular terms in (3) and (4) is of course 
fatal to the theory, since £; and A; will evidently 
be singular as x0 and no provision has been 
made to remove such singularities. Moreover, in 
the derivation of the subtractative terms ex- 
plicit use is made of the relation —A; = E;. One 
would be just as badly off if one attempted to 
regard the subtractative terms in (1) as not being 
explicit functions of dynamical variables, since 
then the subtractative terms for the current 
would be lacking in the Hamiltonian equations 
for the field. 

Heisenberg has applied the Hamiltonian func- 
tion (1) to the calculation of the proper energy of 
a photon in order e* but he did not obtain a finite 
result. The reason for this is easily seen. Since 
is a function of x the transformation matrix 
which diagonalizes it will also depend on x. 
The transformed Hamiltonian function is H(x) 
= S~'(x)/7(x)S(x). However Heisenberg set x =0 
taking (x) 
= S~'(0)/7(x)S(0), presumably in order to avoid 


in the transformation matrix, 
any possible consequences of the change from the 
Dirac equation to (2). (2) of course reduces to the 
x=0. 
would be unnecessary if S(O) led to a finite result, 


Dirac equation when This precaution 


since in that circumstance it would be quite im- 
material whether one used S(0) or S(x). Actually, 
if S(O) is used, one term appears in //, namely, 
S~'(0)HewS(0), which is completely independent 


of x, and this term diverges.’ Since the off- 


* It would be of no aid to introduce x into E?+H? in a 
way analogous to that in which it is introduced into y*y. If 
the propagation vector of the photon is k, such a modifica- 
tion would introduce into Hgy a factor depending on k-x 
Such a factor does not help, since by the conservation laws 
k is equal to the sum of the momenta of the pair virtually) 
created by the perturbation, whereas the singularities arise 
from integration over the difference of the momenta. 
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diagonal matrix elements of H, as well as the 
diagonal elements, approach no limit as x0, the 
artifice of setting x=0 in the transformation 
matrix is clearly not permissible, nor can any 
conclusions regarding the characteristic values 
of JJ be drawn from such a calculation. 

If one carries out the calculation using S(x) no 
difficulties are encountered as long as x~0, since 
now every term in // involves x. The photon 
proper energy in order e’ is found to be 


W = —(ahv/3r)(z/R?—2 log }CR+O(R)), 
(log C=0.577 ---) (5) 


where z is the component of x in the direction of 
the electric vector of the photon, and R is the 
length of x. But in the limit x—0, the photon 
proper energy becomes infinite. 

The entire contribution to (5) comes from 
IT, in fact (5) is just the expectation value of 
Il zy —hv. This seems paradoxical, since the un- 
wanted singular terms in (3) and (4) do not 
affect the expectation value of /7g¢y in order e?, 
while from the remaining terms all singularities 
have presumably been removed. That the singu- 
lar terms in (3) and (4) play no part may be seen 
most easily by inspection of the ordinary per- 
turbation theory formula for Hzgy. If we write 
H=H®+eH”™+eH™ we see that, since Wry is 
diagonal in the original representation, its ex- 
pectation value in order e? is completely inde- 
pendent of H®. However, the singular terms in 
(3) and (4) can be regarded as coming entirely 
from JJ, 

Heisenberg has ascribed the ‘photon proper 
energy difficulties to the point nature of the 
photon, that is, to the possibility of large fluctua- 
tions in the quantized electromagnetic field. That 
this is not the true explanation can be seen by 
thinking of the calculation as carried out by 
direct integration of the equations of motion in 
powers of e. One then sees that in order e? the 
energy contains no terms higher in degree than 
quadratic in the unperturbed value of the electro- 
magnetic field variables, ),, and },,*. For the 
photon proper energy, we require the terms in the 
energy proportional to My, the number of 
photons. Fluctuations in the field consequently 
contribute nothing : one must have terms at least 
quartic in }, and },,* before fluctuation terms 
proportional to \/,, can appear. 
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The real source of the trouble lies in the fact 
that a Hamiltonian perturbation calculation jn 
powers of e is equivalent to integration of the 
equations of motion under the supposition that e 
is a slowly varying function of the time, which 
increases from an initial value zero to its final 
value, e. But in the calculation of the subtracta- 
tive terms it is explicitly assumed that e is inde- 
pendent of time: all time dependences are to be 
included in the vector potential A. If it were not 
for the subtractative terms it would be quite 
immaterial, in computing the current, which one 
supposed to vary, e or the electromagnetic fields, 
since only the produce eA appears in the Dirac 
equation. However the subtractative terms in the 
current also involve eA, and hence will by no 
means be the same in the two cases. Suppose that 
the perturbation is applied over a time 7. During 
this time transient terms of order 1/7 will appear 
in the expression for the induced current, before 
inclusion of the subtractative terms. If one 
supposes A to vary with time, transient fields 
will also be present, which will contribute 
additional subtractative terms which at every 
instant just cancel the transient terms in the 
uncorrected current. However if e is varied 
these additional subtractative terms will not 
appear. The transient terms which in conse- 
quence remain in the current vanish in the limit 
T— ~«, but since they also act for a time 7, their 
contribution to the electromagnetic fields will be 
of order unity as 7. If one calculates the 
contribution of these transient terms to //gy, one 
obtains just the result (5). The photon proper 
energy calculation illustrates a further difficulty 
involved in the passage to a Hamiltonian theory: 
even though it is possible to give a rule by which 
non-singular solutions of the equations of motion 
may be obtained, these solutions will not, in 
general, correspond to stationary states of the 
Hamiltonian system. 

The efficacy of the subtractative terms is in 
any event severely limited by the fact that they 
have been determined only in the approximation 
in which the density-matrix treatment is appli- 
cable, and so cannot be expected to give, in 
general, finite answers for problems in which 
fluctuations are involved. 

The possibility of constructing a Hamiltonian 
theory, along the lines laid down by Heisenberg, 
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which would avoid the difficulties discussed 
above is very doubtful. Certainly no simple modi- 
fication of the original scheme, such as the 
introduction of new subtractative terms in the 
Hamiltonian function, will give equations of 
motion and expectation values for energy which 
are both free of singularities. On the other hand 
formal arguments alone hardly provide an 
adequate basis for a more radical revision. 

In the absence of any consistent theory, the 
simplest procedure which suggests itself is to 
dispense entirely with subtractative terms, to 
carry out all calculations with x different from 
zero, and in the result to pick out the term inde- 
pendent of x as the only one of significance. This 
scheme has the flaw that a classification of terms 
as depending on x or not depending on x is 
essentially ambiguous. It is of course possible to 
pick out terms which are singular or which vanish 
as x—0, but unfortunately, (x*+? 
+2*)/R?=1, an unambiguous separation into 


since, 


“determinate” and ‘‘indeterminate”’ terms is not 
possible. In this connection, it should be re- 
membered that after the introduction of x, an 
operator F is in general no longer either Lorentz 
covariant or gauge invariant. This could lead to 
no inconsistencies provided the lack of the proper 
transformation properties manifested itself only 
in terms which were recognizably x-dependent. 
However, we have not succeeded in formulating a 
criterion for the seperation of determinate and 
indeterminate terms which, in any given problem, 
can be relied upon to lead to Lorentz covariant 
and gauge invariant results. Calculations made in 
this manner accordingly cannot lead to un- 
equivocal results. But it may be hoped that they 
will give some indication of the sort of corrections 
which must be applied to present theory, and of 
the order of magnitude of the correction terms. 

We shall introduce x into the theory in a 
somewhat different way than did Dirac and 
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Heisenberg. The operators of the usual field 
theory will be left altogether unaltered, except 
for their symmetrization in electron and positron 
(WWt(X)W(X) 3 Lt (X)W(X) —W(X)yt(X) ]). The 
commutation laws however, will be taken to 
involve x: 


[wt(X), W(X’) }+ =8(X —X’+2). 


This formalism would correspond to the im- 
possibility of making field measurements without 
mutual disturbance at two arbitrarily near 
points. It thus suggests a not unreasonable 
motive for the introduction of the off-diagonal 
distance x. 

As in the usual quantum electrodynamics, the 
wave equation associated with the Hamiltonian 
function is to be augmented by the auxiliary 
relation div E—4rp=0. This relation is, in fact, 
not compatible with the equations of motion, 
since the conservation law for the current holds 
only when x=0. The inconsistency of the auxili- 
ary relation with the equations of motion is only 
to be expected, for, with x #0, the Hamiltonian 
function is no longer gauge invariant. The same 
inconsistency appears if x is introduced into y 
and y* instead of into the commutation relations 
(Even in the density matrix treatment, it will be 
remembered, the conservation laws fail for x #0). 
Since the covariance of an inconsistent set of 
equations is in any event illusory, we have been 
content to introduce x into the theory in a 
patently noncovariant manner. As we shall see, 
this lack of covariance is much the same as that 
which appears in classical electrodynamics when 
one attempts to introduce the radius of the 
electron. 

The formalism just described finds an interest- 
ing application in the calculation of proper 
energies of photon and electron in order e*®. The 
photon proper energy is found to be (in c.g.s 
units) 


9.9 


3g 3? 1 =x? 
+0O(R) 


where x is the component of x in the direction of propogation of the photon, z is the component in 


the direction of the photon’s electric vector, and R 


is the length of x. This result may reasonably be 


interpreted as indicating a null photon proper energy in order e°. 
If we carry out the elimination of the electrostatic field variables in the usual way, overlooking, of 
course, the inconsistency of the auxiliary relation, we obtain in the energy a term 
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1 > (nt dm — Oman*) (As+a,—Aias* )Onm: sts 6 
where 

Unt (X )tm(X )ust( X’)o(X’ 
Oam: 0 = 50" fax [ax exp — }i(K,+k,,+k,+k,)-x, 
. . X—X 

and y and ¥* have been expanded in electronic plane wave functions v,, of momentum k,,: Y= Ya,,, 
exp — 3tk,-x, += doa, *v,,* exp — }7k,,-x. The factors exp — }7k,-x are inserted so that the a, and a 


n 


will obey the usual commutation relations. The diagonal matrix elements of (6) give the ordinary 

Coulomb and exchange terms, and in addition, terms which represent the electrostatic proper energies 

of electrons and positrons. The proper energy of an electron in a state r is 
Wes = LQr'; r'r— LOre: or: 

Here r and r’ refer to positive energy states, p to negative energy states. Evaluating (7) in a Lorentz 


frame in which the electron is at rest, we find 
Wes=(e? 16x?) fdr r+x fa 1+?)-* exp (tk: r) = (e?/47? fark rr+x 


where K,(r) is the Bessel function defined in Whittaker and Watson, Modern Analysts, §17.71. The 
angular integration gives 


R _ 
Wes=(e? ac ‘| rK\(r)dr- K,(r)d) 
‘ Jp 


0 


which can be reduced, by use of the recurrence relation K,(r)=(d/dr)Ko(r), to 


»f 
Wis =(e? rR) | K(r)dr= — (2e?/r)[ log }CR—1+0(R) }. 


| 


The electromagnetic proper energy in order e?can Our method of calculation, too, is strongly 
also be shown to be finite when x#0. Thus the _ reminiscént of that employed in classical electro- 
dynamics, which can likewise be described as the 
introduction of a parameter, the electron radius, 
and the final abandonment of all terms except 


those independent of this parameter. The aim, of 


artifice introduced to eliminate the singularities 
due to the sea of negative energy electrons has a 
very natural extension to the treatment of the 
electronic proper energy difficulties. 

It will be observed that the filling of the nega- course, is to separate, insofar as possible, those 
tive energy states results in the summations over — terms which obviously depend on the particular 
positive and negative states appearing in the 
electrostatic proper energy with opposite sign. If 
the negative energy states were left empty, independent of these factors, and which may be 
the proper energy would instead be 2/Q,,; +, hoped therefore to have some significance, 
+50,» »-. This gives Wes=e!/2R, which is fant despite the unquestionable inadequacy of the 


model employed or on the way in which calcula- 
tions are carried out, from those which seem to be 


theory. 

I wish to express my deepest thanks to Pro- 
fessor J]. R. Oppenheimer, to whom I am indebted 
for many illuminating discussions of this subject. 


p 


the electrostatic energy of an electron of radius 
R. The off-diagonal distance x thus appears to 
find a classical analogue in the electron radius. 
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Note on Neutron-Proton Exchange Interaction 


[he matrix elements of the interaction between a proton 


with coordinates x; and a neutron with coordinates X, as 


proposed by Majorana' may be written as 
x xo) V)x,'xo"’) = —J (|x, —Xe2| )6(x, —X2")5(X2—xX;,’). 1 
his interaction assumption corresponds to an interchange 
of the positional coordinates (x,, X2) of the two particles; 
clearly this interchange will entail a shift in the center of 
mass if the mass of the proton, m1, is not the same as the 
mass of the neutron, m,. A similar effect on the center of 
mass is also characteristic of the Heisenberg interaction 
assumption? in which both spin and positional coordinates 
‘‘mixed”” Heisenberg- 


and of any 


While the difference between the 


are interchanged, 
Majorana interaction 
masses of the proton and neutron according to present 
values is small, the difficulty introduced by the shift in 
the center of mass has perhaps a theoretical interest and 
has been discussed by Breit and Wigner.’ 
This center of mass difficulty is made at once evident by 
considering the matrix elements of the exchange interaction 
in terms of momentum variables. Thus, the transformation 


of (1) gives 


Pip: |" pips’) =f( p2’— pi )6(pi’—pit+p2’—pz), (2 


in which p; isa momentum of the proton, and p» a neutron 


momentum and the form of f depends on the particular 
for J in (1). If the matrix elements (2 
written as 


total 


form assumed 


are now functions of the momentum 
P=p,+p: and the relative momentum p = u(p2/m2—p,/m 


where «=m m2/(m, 2) becomes 


pP I’ p'P 
=f( p+ p’ ip Me m ml 


+ 91 


+ Me 6(P—P’), (3 
in which form the flaw in the exchange interaction assump- 
tion is obvious: the matrix elements of the interaction are 
explicit functions of the momentum P of the center of mass 
so that it is not possible to separate the energy or the wave 
function characteristic, respectively, of the 


into parts 


internal motion and the motion as a whole. (3) is to be 


regarded as a consequent interaction law only insofar as 


one may put m,=m», for which case 
pP | V p’P’)=f('p+p’|)s(P—P’), 3’ 

and the difficulty mentioned above disappears. 
There are many possible ways of modifying the inter- 
action assumption (1) or (2) without losing the essential 
feature of the exchange, and without introducing at the 


same time any error due to a difference in masses. One 


on 


tn 


simple possibility is to take instead of (2) the following 
pip: V'| pi'ps') =f( ——| mi(pe' +P 

ms(pi' +p )o(p\’-pitp Pp 4 
Phen, pP V p’P’)=f( p+p’! )4(P—P’). 5 


and (3 
that the wave functions ¥(p) and the energy values E(P 


rhe identity of (5 proves the interesting result 


found by using the Majorana interaction under the as 
sumption m, =m, are exactly the same as the corresponding 
wave functions and energy values deduced from the 
modified interaction which is valid regardless of the values 
of m, and mz. 

[he same conclusions as given above should of course 
follow from consideration of the interaction in a coordinate 


representation. In fact, the coordinate representation of 


4) is 
X\X2| V | x,'x.’) = —J(|x,—x, 
. ” 
o| x Xe + x x 
” aL 
m ml 
0] Xe—X; + x x 6 
nl + Ms 


which is the same interaction as has been proposed by Breit 
and Wigner. If 


proton is written as y =u(X)g(é 


the wave tunction tor a neutron and a 


where X = (\x,+ ox 


mi+my) and &=x;—X.2, and if V is the operator defined 
by (6), V-u(X')g(E’) = —J(\E! )u(X)g(—E), while the 
operator defined by (1) is such that V-u(X’)g(& 


I(\E| (K+ (m2—m,)E/(me+m,))g(—-E 

It is clear that it is possible to extend the exchange 
interaction law so as to include the case of two particles of 
different masses without adding any essential complication. 
Che greater complexity of (6) as compared with (1) is 
only apparent and should be regarded as characteristic of 
the particular coordinate system used. Although 
that 


there 


may be other reasons for supposing the exchange 


interaction is accurate only insolar as (m, Me 


i M+ Me 
can be neglected in comparison with unity, it does not 
appear that such a limitation® can be found by examining 
the exchange hypothesis itself for consistency and sim- 
plicity. 
MiLTon S. PLESSET 
University of Rochester, 
Rochester, N. Y 
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March 11, 
E. Majorana, Zeits. f. Physik 82, 137 (1933 
?W. Heisenberg, Zeits, f. Physik 77, 1 (1932); 78, 150 (1932). 
G. Breit and E. Wigner, Phys. Rev. 48, 918 (1935). 
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The F Terms of Ce IV 


Following the publication by the author of a partial 
analysis! of the spectrum of Ce IV considerable interest 
appeared concerning the position of the F terms. These 
terms have now been located and this opportunity is 
taken to give a brief note concerning them. A rather 
complete analysis of the spectrum is now nearing com 
pletion and details will be published elsewhere. 

Ihe 5F terms are quite regular as compared with neigh- 
boring ions of similar electronic structure and are given 


by the following lines: 


$Dy—5Fy 741.79 (40) 134,809 
5Dy—5Fy 754.60 (30) 132,520 
§Dy4—5Fy 755.75 (12) 132,319 


[he question whether the 4F or the 5D terms are the 
deepest in the spectrum is a very interesting one and seems 
now to be definitely settled in favor of the former possi- 
bility. The lines involved are: 

4Fy4—5Dy 2009.94 (100) 49,737 
4Fy—5D., 2000.42 (100) 49,973 
4Fy—5Dy 1914.75 (35) 52,226. 

rhe following term values are based upon three S terms 

using a Ritz formula. 


4Fy 296,197 
4Fy 293,944 


5Fy 111,652 
SFy 111,451. 


The reality of these terms is confirmed by the 5G 
combinations and also by the fact that the 4F interval 
corresponds sufficiently well with the proper intervals in 
the low terms of Ce III. 

R. J. LANG 

University of Alberta, 


Edmonton, Alberta, Canada, 
March 10, 1936. 


R. J. Lang, Can. J. Research Al3, 1 (1935). 


Magnetic Moment of the Li’ Nucleus 


The value of the magnetic moment of the Li? nucleus has 
been determined in two different ways. Breit and Doer- 
mann! calculated it to be 3.28 nuclear magnetons, on the 
basis of hyperfine structure measurements of the *P)—*S;, 
group of Lit. Fox and Rabi,? by the molecular beam 
method, found the hfs separation of the ground state (*.S;) 
of the normal Li atom to be A(s) =0.0267+0.0003 cm™. 
Che resulting value of the nuclear magnetic moment, 
estimated from the Goudsmit-Fermi-Segré formulae, was 
found to be 3.20 nuclear magnetons. However, even apart 
from other objections’ to the GFS formulae, it is not to be 
expected that a statistical theory will yield correct results 
for an atom as light as lithium. Therefore, a more accurate 
calculation is desirable. 

Recently, Fock and Petrashen‘ have published values of 
wave functions (self-consistent with exchange) for lithium. 
For the 2s electron, the value of x=(foo/r) at the origin 
is 1.44. The value of the nuclear magnetic moment, in 
terms of nuclear magnetons, is’ g=949 A(s)[J/(27+1)] 


THE EDITOR 


x [1/x2(0) =4.58. This is in disagreement with the above 


value of 3.28 as found by Breit and Doermann. 

Fock and Petrashen have also calculated the transition 
probabilities for some members of the principal series of 
lithium. The agreement with observation® is mostly 
qualitative, and not quantitative. This work indicates 
then, that the Fock functions are inadequate for transition 
probability calculations. Our work shows that the same is 
true for nuclear moment calculations, if the theory be 
supposed to be correct. 

It may be possible, with more accurate wave functions, 
to resolve the above discrepancy between 4.58 and 3.28 
rhis is being investigated. 

James H. BartLert, Jr 
J. J. GrBBons, Jr. 
Department of Physics, 


University of Illinois, 
March 6, 1936. 


Breit and Doermann, Phys. Rev. 36, 1262 (1930 
* Fox and Rabi, Phys. Rev. 48, 746 (1935). 
’ Gibbons and Bartlett, Phys. Rev. 47, 692 (1935 
‘Fock and Petrashen, Physik. Zeits. Sowjetunion 8, 547 (1935 
>A. Filippov, Zeits. f. Physik 69, 526 (1931). The disagreement 
between these observations and Fock's theory seems to be quits tsicke 
the experimental error. 


A New Form of Crystalline Quartz at —183.5°C 


When X, Y and specially oriented Y cuts,' encased in 
an air-tight brass crystal holder (Fig. 1), were subjected 
to low temperatures, they ceased abruptly to oscillate 
piezoelectrically near the boiling point of liquid air. In one 
set of experiments these cuts of quartz were driven by a 
Hartley oscillator. Resonance responses of the plates were 
noted by listening for the corresponding characteristic 
disturbances produced in the horn of an amplifier which 
was connected across a portion of the Hartley plate 
circuit. This method is perhaps the most sensitive test? 
for the piezoelectric property. In another set of experiments 
specially oriented Y cuts were placed in the usual nonre- 
generative crystal oscillator circuits. These circuits ceased 
abruptly to oscillate near the boiling point of liquid air. 
Piezoelectric activity increased steadily from room tem- 
perature to approximately —175°C. By immersing the 
crystal holder in liquid oxygen boiled under reduced 
pressure, the temperature at which alpha quartz is trans- 
formed to the nonpiezoelectric form, called 6 quartz, was 
found to be —183.5°C. The temperature of the ja 
transition occurs above the boiling point of liquid oxygen 


at atmospheric pressure. The range permitted by our 
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Fic. 1. Crystal holder. B, one of six brass bolts; H, Hartley oscillat 
/, insulators; Q, quartz plate 














Polarimeter crystal er. Q, quartz plate; / ypt ixis 
e; 7, insulator for mounting one of the brass electrodes 
pools; F, iron tu G 
H, Hartley oscillator 





apparatus showed that this transition must occur above 
—182.4°C. The existence of two separate transition tem- 
peratures is perhaps due to the tendency of either state to 
persist. 

In the most satisfactory arrangement for studying the 
discontinuity in the optical rotatory power at the a+~é 
transition the quartz plate was mounted in a holder which 
is incompletely described in Fig. 2. This holder was 
mounted in a Lippich three-prism polarimeter.’ For the 
5461A line of Hg the optical rotatory power of both right 
and left-handed quartz decreased from 25.59°/mm at 
room temperature to a minimum of 24.89°/mm at the 
as) transition. Here, optically right or left-handed 
quartz remained optically right or left handed, respectively. 
Below the transition point the rotatory power of 5 quartz 
increased to 24.95°/mm at the lowest temperature (perhaps 
—190°C). In calculating these rotatory powers the z 
dimensions of the plates were assumed to remain constant. 
Changes in the rotatory power were more abrupt at the 
a—§ transition than at the é—a transition. This was also 
true as regards the piezoelectric property, but the changes 
in rotatory power were sluggish as compared to the rapid 
piezoelectric changes. Both of these changes occurred 
together. Corresponding changes in volume must be small 


since quartz can be plunged into liquid air without the 


slightest damage. 
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THE EDITOR 


Upon the assumption that 6 quartz is a member of the 
hexagonal system, the loss of the piezoelectric property is 
interpreted to mean that it falls among classes 9, 12, 21 or 
24 (Voigt) or 24, 25, 30 or 32 (Schoenflies). Of the four 
types (a, 8, y‘ and 54) of crystalline quartz y and 6 quartz 
are nonpiezoelectric. 

lime was lacking for determining the 5—a transition 
temperature and measuring the optical rotatory power ot 
5 quartz as a function of temperature. 

My thanks are due to Robert E. Adams for his assistance 
in this work and to Professor J. R. Roebuck for helpful 
discussions. 


HAROLD OSTERBERG 
University of Wisconsin, 
March 10, 1936. 
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Errata: Continuation of Work on Scattering of Slow 
Neutrons 


ALLAN C. G. MITCHELL, EDGAR J. MuRPHY AND 
LAWRENCE M. LANGER, New York University 
Phys. Rev. 49, 400 (1936 


AND 
Selective Scattering of Slow Neutrons 


ALLAN C. G. MITCHELL, EDGAR J. MURPHY AND 
MARTIN D. WHITAKER, New York University 
Phys. Rev. 49, 401 (1936 
Through an inadvertence Fig. 1 of the first letter appears 
with the caption and text of the second letter on page 401 
and Fig. 1 of the second letter appears with the caption and 
text of the first letter on page 400. An interchange of these 


two figures should be made to rectify the error 








